FERMILAB-PUB-08/028-E 



Review of Top Quark Physics Results 



R. Kehoe 

Physics Department, Southern Methodist University, Dallas, Texas 
Dallas, TX 75275, USA 
kehoe@physics.smu.edu 

M. Narain 

Physics Department, Brown University 
Providence, HI, 02912 USA 
narain@hep. brown, edu 

A. Kumar 

Physics Department, State University of New York at Buffalo 
Buffalo, New York 14260, USA 
ashishk@fnal.gov 

Received 14 December 2007 
Published 12 February 2008 

As the heaviest known fundamental particle, the top quark has taken a central role 
in the study of fundamental interactions. Production of top quarks in pairs provides 
an important probe of strong interactions. The top quark mass is a key fundamental 
parameter which places a valuable constraint on the Higgs boson mass and electroweak 
symmetry breaking. Observations of the relative rates and kinematics of top quark final 
states constrain potential new physics. In many cases, the tests available with study 
of the top quark are both critical and unique. Large increases in data samples from 
the Fermilab Tevatron have been coupled with major improvements in experimental 
techniques to produce many new precision measurements of the top quark. The first 
direct evidence for electroweak production of top quarks has been obtained, with a 
resulting direct determination of Vjj,. Several of the properties of the top quark have 
been measured. Progress has also been made in obtaining improved limits on potential 
anomalous production and decay mechanisms. This review presents an overview of recent 
theoretical and experimental developments in this field. We also provide a brief discussion 
of the implications for further efforts. 
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1. Introduction 

The top qu ark has played a key role in particle physics for well over a decade. 
Its discoverjlH2J marks a triumph of the modern theoretical framework of SU (3) x 
SU(2) x U(l) interactions. It also hallmarked a new generation of experimental ap- 
proaches that will be central to the intensifying search for the Higgs boson and new 
physics. Subsequent analyses have measured the mass of the top quark with remark- 
able precision. What is more, this particle is uniquely placed to probe strong inter- 
actions, electroweak interactions, and the Higgs mechanism itself. Many challenges 
confront the experimentalist, particularly because of the diversity and complexity 
of the final state. 

We review this exciting new field with the underlying goal of illustrating not only 
the measurements themselves, but also the techniques being developed to extract 
them. In the remainder of this section, we synopsizc the role of the top quark within 
our understanding of fundamental interactions and particles. In subsequent sections 
we discuss production and decay expectations, as well as the experiments that study 
top quark events. Measurements of strong and electroweak production are reviewed 
in Section |4] and El respectively. Section [6] treats precision measurement of the top 
quark mass. This is followed in Section [7] with coverage of other properties of this 
particle. Lastly, we finish in Section [8] with searches for non-standard physics within 
the top quark sector and we conclude with some considerations for future efforts. 
There are excellent reviews of this subj ect which primarily cover results obtained 
at y/s = 1.8 TeV from 1992-1996 ISM termed Run I of the Fermilab Tevatron. 
Informative reviews with a more phenomenological emphasis are also availabl e ! ^ !. 
We concentrate in this review on results obtained at y/s — 1.96 TeV and much 
higher integrated luminosity in the Tevatron's Run II which started in 2001. We 
also make a point of including the later Run I results that may not have been 
completely discussed before. To provide completeness and also a better indication 
of the direction of the field, we include several preliminary results from CDF and 
D0 when published results are not yet available. As these results may change before 
publication, we indicate them in the text. 

1.1. Electroweak and strong interactions 

We now view the weak and electromagnetic interactions as different manifestations 
of the same underlying electroweak interaction. The crucial progress in this realiza- 
tion came with the development of a viable SU(2) xU(l) modeW This interaction is 
mediated by four gauge bosons. The W ± bosons mediate the charge changing weak 
currents for left-handed fermions. The Z boson and the photon propagate the weak 
and electromagnetic neutral currents. The symmetry of the theory's Lagrangian is 
broken in the ground state of the system by introducing the Higgs fieloP which has 
a non-zero vacuum expectation value. This then provides mass to the W ± and Z 
bosons. This model, originally designed around first and second generation leptons, 
was extended first to include two generations of quarks, and then a full third gen- 
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eration of fundamental fermions. The first generation consists of the electron (e) 
and its neutrino {y e ) and the up (u) and down (d) quarks. The second generation 
includes the muon (fx), muon neutrino (l/p), charm (c) quark, and strange (s) quark 
and the third generation consists of the tau (r), tau neutrino (v T ), top (i) quark, and 
bottom (b) quark. In each generation, the left-handed charged lepton and neutrino 
form a doublet under SU(2) and so do the left-handed quarks. The quantum num- 
ber associated with the SU(2) symmetry is called weak isospin and the left-handed 
fermions have a third component of weak isospin, t^L = Right-handed fermions 
are isosinglets. The Higgs couples to these fundamental fermions / with Yukawa 
couplings, yf. Each fermion propagates through the Higgs field and, by virtue of 
this coupling, acquires mass. The value of yf is therefore related to the mass of the 
fermion, m/, by yf — y/2mf/v, where v is the vacuum expectation value of the 
Higgs field. The value of yf is different for each fermion and is not predicted by the 
model. We can only ascertain its value by a measurement of each fermion's mass. 

Strong interactions were given a coherent theoretical description by the non- 
abelian SU (3) theory of quantum chromodynamics in conjunction with 

the quark model of hadronJ^. Eight massless gauge bosons (gluons) propagate 
the interaction and carry a color charge. Quarks are the underlying constituents of 
hadrons. Their quantum numbers explain the isospin and baryon number proper- 
ties of hadrons in nature. The strength of the strong interaction, governed by a s , 
is large, making perturbative calculations difficult. The attraction between quarks 
exhibits asymptotic freedom, meaning that at large energies or very small distances, 
the interaction strength declines. The converse means that isolated quarks are pro- 
hibited from being extracted from hadrons. This 'confinement' means that bare 
colored particles are not observable in nature. Two colored partons receding from 
each other will produce a series of colorless hadrons in the final state. This process is 
called hadronization. These hadrons are collimated into jets moving roughly in the 
directions of the original par tons. The process is understood phenomenologically by 
several models, eg. Ref. ^4 In the string fragmentation model, a narrow tube or 
'string' connects two colored partons. As they recede, particles are created along 
the string which give rise to jets and energy flow along lines of color in an event. 
The physics of hadronization provides great challenges for precise measurements 
that involve jets, like those of the top quark. 

1.2. Top quark and the flavor spectrum 

The existence of the top quark is firmly placed in our picture of fundamental inter- 
actions. In our current understanding, all matter is made of fundamental fermion 
fields in two categories. Quarks are sensitive to the strong interaction; leptons are 
not. The electroweak eigenstates of quarks are a mix of the mass eigenstates, de- 
scribed by the CKM matrix-121. As a result, generations that contain only a single 
quark lead to the expectation of flavor changing neutral current interactions. The 
absence of such interactions in experimental observations motivated the successful 
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prediction of the charm quarlffl 

The discovery of the third generation tau lepton (rjpl anc j bottom quark (6^21 
led to the search for the top quark (t) as the SU(2) partner of the bottom quark. 
Experimentally, several measurements showed that the bottom quark was a mem- 
ber of an SU(2) doublet. Several early modelfPHSni allowed for the absence of 
a sixth quark. If the bottom quark were an SU{2) singlet, then flavor changing 
neutral current interactions would result in the B meson system. The limit on the 
branching fraction for B —> < 10~ 3 , h owe ver, put a stringent limit on th ese 

processes and ruled out the singlet hypothesis^! Measurements of B°-B° mixing^l 
are sensitive to the magnitude of the CKM matrix elements Vtd and Vt s . These were 

observed to be non-zero, indicating the existence of the top quark, although this 
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could also be explained without the presence of a top quark^. Measurements of the 
third component of the weak isospin of the b— quark, t\ L: provided strong additional 
evidence. The b— quark coupling to the Z boson is dependent on t\ L + ^sin 2 ^, 
where 6 W is the weak mixing angle. Measur ement s of the Z — > ra te, as well as 
the bb charge asymmetry in e + e~ collisioniP^S] both probed this coupling. Both 



measurements indicated t\ L = — necessitating a t^L = +i partner to the b 
quark. A model was develop eoP2l to explain the Z width measurement without the 
top quark, but this could not account for the charge asymmetry measurement. 

Theoretical arguments also favored the existence of a top quark once the r lepton 
was identified. In the SU(2) x U(l) electroweak model, triangle diagrams with V 
and A couplings give rise to anomalies unless the total charge in a given generation 
sums to zero. With three color degrees of freedom per quark, anomalies are avoided 
in the first two generations. The existence of the r (and its surmised neutrino) 
indicated the presence of a third generation quark doublet. Once the b quark was 
found, another quark of charge = +| was required. 

As a result, a search for the top quark was pursued from the late 1970's onward 
(e.g. Refs. I27l- l3i"jl . As accelerator energies continued to rise, the lower limits on 
the top quark mass increased. A fuller description of this period of research can be 
found in Ref. [3l By the mid-1990's, it was established that the top quark was very 
heavySH. This had the result of focusing the search strategies of the CDF and D0 
collaborations and led to the discovery of top-antitop quark pair production in 1995 
at the TevatroiJH21. 



1.3. Top quark properties 

A massive top quark has several properties which make it quite interesting as a probe 
of known strong and electroweak physics, as well as a sensitive window to potential 
new physics. Perturbative QCD calculations of pair production can be carried out 
with significant precision. Electroweak production of single top quarks gives direct 
sensitivity to \Vtb\- The value of m t places an important constraint on the mass 
of the Higgs boson. It has been speculated that the resulting top quark Yukawa 
coupling, could point to new dynamics beyond the standard modeP^I The 
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exceedingly short predicted lifetime of the top quark of O(10 _25 s) is an order of 
magnitude smaller than hadronization timescales and this permits a glimpse of the 
properties of a bare quark. The decay modes of the top quark may harbor evidence of 
new physics if alternative models are correct. Measuring the charge of the top quark 
is important to establish that the top quark is as expected. A direct determination 
of the \Vtb\ matrix element probes for potential new physics. 

2. Production and Decay of the Top Quark 

For the next several years, the top quark will continue to be produced solely in 
hadron collisions, specifically involving protons. In such collisions, the top quark 
can be produced strongly in pairs or electroweakly alone. In general, the proton 
can be considered to harbor three 'valence' quarks (uud) which dictate its quantum 
numbers. These valence quarks typically carry much of the momentum of the proton. 
There are also virtual or 'sea' quarks and gluons in the proton which carry less 
momentum individually. When a proton and antiproton collide, a hard interaction 
occurs between one of the constituents ('partons') of the proton with a parton of 
the antiproton. Soft interactions involving the remainder of the hadron constituents 
produce many low energy particles which are largely uncorrelated with the hard 
collision. Because of the large mass of the top quark, production usually involves 
the higher momentum valence quarks at the Tevatron. 

2.1. Strong pair production 

Production of a top-antitop quark pair (ft) occurs dominantly via strong processes. 
At leading order (LO), valence quarks supply the primary production probability 
at the Tevatron via quark-antiquark annihilation (qq). Approximately 15% of the 
cross section comes from gluon (gg) fusion. Fig. [T] shows these diagrams. 

In the center-of-mass frame in which the proton and antiproton are rapidly mov- 
ing, the hard interactions between constituent partons are fast relative to the time 
for partons to interact. As a result, the hadronic collision can be factorized^l into 
a parton collision weighted by 'parton distribution functions' (pdf's), Fj(xj) which 
express the probability for parton i to carry momentum fraction, Xi, of its parent 
hadron. These pdf's are properties of specific hadrons and are independent of the 
specific hard scatter interaction at parton level. They encompass non-perturbative 
soft processes. As a result, they are extracted from the examination of inelastic 
interactions involving hadrons. The pair production cross section is then calculated 
as 

a{pp—f tt + X) = Sjj J dxidxj x F[(x it Uf)Fj{xj, (j,f)aij(xi,Xj,m%, a}) (1) 

where the sum runs over gluons and light quarks in the colliding proton and antipro- 
ton, and <xy is the perturbative cross section for collisions of partons i and j. The 
factorization scale, fiF, defines the splitting of perturbative and non-perturbative 
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Fig. 1. LO Feynman diagrams for quark-antiquark annihilation (qq) and gluon fusion (gq). 

elements. It is common to treat this scale in common with the scale appropriate 
to the renormalization of the perturbative cross section, /x^j, since both parameters 
are arbitrary. This common scale, /i, is usually taken to be equal to the top quark 
mass (mt)- An exact calculation would not depend on these scales. Finite order 
calculations have a sensitivity that must be assessed as an uncertainty in the theo- 
retical calculation, usually by bounding the predictions with /i = m t /2 and /i = 2m t 
calculations. 

Initial LO cross secti on cal culations were performed in Ref. 34, Next-to-leading 
order (NLO) calculations^^! accounted for associated quark production and gluon 
bremstrahlung, and virtual contributions to the LO processes. Because top quark 
production at the Tevatron occurs near threshold, a large uncertainty in these cal- 
culations results from initial state soft gluon radiation. In order to fully account for 
this, several groups included resummation of dominant soft logarithms to all orders 

in perturbation theor y for hadronic cross sections. Results from resummed Drell- 
1071001 

Yan productio n 11 were generalized to handle the color elements in the initial 
heavy quark calculations^' - '^. While the corrections are small for ti production, 
the dependence of the cross section calculation on the choice of scale is reduced. 
This work has been improved with more recent pdf's and m ore accurate uncertain- 
ties in the determination of the pair production cross sectiorl^. The results of this 
calculation for Tevatron energies are given in Table[T] A total uncertainty of 10% to 
15% is obtained for the NLO cross sectiorl^. Parton distribution functions, partic- 
ularly for gluons at high x, provide the main source of theoretical uncertainty. The 
modest increase in collision energy from 1.8 TeV to 1.96 TeV causes a 30% increase 
in the production cross section. This is because of the substantial gain in valence 
quark phase space from the parton distribution functions. 

Resummed calculations have also been performed including additional higher 
order termJ22Il2] Both Refs. I4T1 and l43l are based on [35l [36] with resummation at 
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Table 1. Calculated it cross sections from Refs. [41. 43 assuming rat = 175GeV/c 2 . 
The former is a complete NLO calculation where the uncertainty is primarily from the 
choice in pdf. The latter provides an NLO calculation with additional higher order terms 
and an uncertainty taken from the choice in kinematic scheme. The rightmost col umn 
provides the change in calculated cross section (Act) appropriate to a change in rat 



(TeV) 


Kidonakis, et al.HS 


Cacciari, et al.ED 


Act (pb) 
{ra t = 170, 180 GevjSH 


1.8 

1.96 


5.24 ± 0.31 pb 
6.77 ± 0.42 pb 


5.i9tg;g pb 
6.7D±°;S pb 


+0.91, -0.76 
+1.13, -0.95 



least to next-to-leading logarithms. A difference between them is that the former 
resums soft gluons to all orders, while the latter ignores small terms beyond next- 
to-next-to leading order. Because top quark production is at threshold and both 
calculations are not complete beyond NLO, there is an ambiguity in the soft-gluon 
resummation calculation. We do not attempt a detailed review of this topic here, 
but in short the situation is the following. These calculations are performed in terms 
of a choice of kinematic parameters, and different choices will give slightly different 
parton level calculations^. To address this, Ref. [43] has expanded the ersummed 
cross section to next-to-next-to leading order. The resulting corrections are small 
for choices of scale = m t , but the observed dependence of the cross section on scale 
and kinematics is further reduced. The results of this calculation for y/s = 1.8 and 
1.96 TeV are given in Table [TJ The remaining difference from choice of kinematics 
provides the uncertainty estimate Ss3. 

Because of its short lifetime, the tt pair is not expected to form a bound meson. 
However, non-standard production mechanisms have been proposed which do form 
a tt resonance, eg. Ref. 1441 Testing such production mechanisms will help determine 
if the top quark fills a special role among fundamental fermions. Such analyses will 
be discussed in Section l8Tl 

2.2. Electroweak production of single top quarks 

In addition to the pair production of top quarks via the strong interaction at the 
Tcvatron, they can also be produced singly in electroweak interactions. There are 
three modes of single top quark production, which differ in the virtuality, Q 2 , of the 
participating W boson, where Q 2 is the negative square of the W boson four mo- 
mentum q. Two of the three modes are labeled by the corresponding Mandelstam 
variables t and s involved in the transition matrix elements. The three dominant 
modes of single top quark production are listed below in descending order of their 
expected production cross section: 

• i-channel: In the process pp — > tqb + X, the W boson is spacelike (— Q 2 = 
q 2 = t < 0). The predicted cross section for this channel, for the Tevatron 
center of mass energy y/s = 1.96 TeV, computed at NLO is 1.98±o'il P*^- 
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Inclusion of NNLO and NNNLO threshold soft-gluon corrections leads to 
a cross section of 2.30±0.14 ptP^H Two Feynman diagrams, as shown in 
Fig-El contribute to this channel. In the leading order diagram, the b quark 
is from the sea of quarks in the proton or anti-proton which couples with 
the virtual W to produce a top quark. In the next-to-leading-order diagram, 
the anti-bottom (6) quark comes from the splitting of the gluon into a bb 
pair. The b quark couples to the virtual W boson and produces the top 
quark. Thus this channel is also known as VF-gluon fusion. The b quark in 
the final state has low transverse momentum (pr) and is at high rj. The 
i-channel denotes the processes with the following quarks in the final state: 
tqb,tq,tqtq. We denote this channel as tqb. 

• s-channel: In the process pp — > tb+X, the W boson is timelike (— Q 2 = q 2 = 
s > (m t + nib) 2 > 0). For this channel the predicted cross section at NLO 
for y/s = 1.96 TeV is 0.88±g;^ prM Corrections at NNLO and NNNLO 
leads to a predicted cross section of 1.08±0.08 The two initial-state 
quarks annihilate into a virtual W boson which decays into a top quark 
and a bottom quark, see Fig. O The s-channel includes both tb and tb and 
is also referred to as W* production. We will refer to this channel as tb. 

• associated single top quark production: In the process pp — > tW, Q 2 = m^, 
an on-shell W boson is produced together with a top quark. The production 
cross section for this channel is predicted to be very small, about 0.28±0.06 
pb at NLO plus NNLO and NNNLO threshold soft-gluon corrections^. 

There are also other diagrams for single top quark production which involve Wts or 
Wtd vertices, but they are highly suppressed due to the small CKM matrix elements. 
Their contribution to the cross section is expected to be less than a percent. 



(a) (4) 




Fig. 2. Feynman diagram for leading order s-channel single top quark production (left) and t- 
channcl single top quark production (middle and right). The middle diagram is the leading order 
and right one is the 0(a s ) W-gluon fusion diagram . 

One of the primary interests in establishing the cross section for single top 
quark production is that its production cross section is directly sensitive to the 
transition width of t — * Wb, and consequently to \Vtb\- The partial decay width 
of the top quark, T(t — > Wb), can be measured, and hence it will be possible to 
measure the top quark lifetime. At the Tevatron, single top quark events have the 
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same final state signature as the production of a Higgs boson in association with 
the W boson, thus understanding single top quark productions is an important 
prerequisite for Higgs searches at the Tevatron. Single top quarks are also one of 
the irreducible backgrounds to associated Higgs production. Studies of single top 
quark events could lead to insights into the existence of non standard couplings (e.g. 
right-handed couplings and polarization) in the top quark decays to W boson and 
6-quark. Single top quark production is very sensitive to new physics. The s and t- 
channels are sensitive to different new physics and hence can be used to distinguish 
between various exotic models (Top Flavor, 4th generation or those which produce 
flavor changing neutral currents etc). 

2.3. Top quark decays and final states 

With three generations of quarks, the unitarity of the CKM matrix^Sl i n conjunction 
with the measured values of V u b and V c b means the decay probability of t — > Wb is 
virtually 100%. The W decay is divided approximately equally between each of the 
three lepton-neutrino (If) pairs and each of the three colored doublets of up-down 
quarks and of charm-strange quarks. The quarks will hadronize to produce jets and 
heavy flavor mesons from b quarks will decay. Typically, the & — jet and W decay 
products will have large momenta in the plane transverse to the beam ('transverse 
momentum', or because of the high value of ra ( . 

Events from top quark pair production consist of two W bosons, and two b- 
quarks: it — ► W + bW~b. Both W's may decay leptonically, one may decay lepton- 
ically and the other hadronically, or both may decay to hadrons. The first case, 
termed 'dilepton', gives a final state of £ + v£~Vbb and manifests itself in the detec- 
tor as two high px leptons, significant unobserved px from the fs ('$r'), plus two 
6-quark jets. The second '^+jets' scenario provides £ ± fqq'bb and is observed as one 
high pt lepton, large $x, an d four jets including two 6-jets. The last decay chain 
is the 'all-jets' channel and it results in a qiq2q 3 q i bb final state appearing as six 
jets including two from 6-quarks. If the leptons are restricted to e and a, the usual 
case, then the total branching fractions are approximately 5%, 30% and 45% for 
dilepton, £+jets and all-jets channels, respectively. Quark color degrees of freedom 
are primarily responsible for the substantially larger branching fraction of ti to the 
£+jets and all-jets final states. 

Final states arising from single top quark production consist of a top quark, a 
bottom quark and a light quark (i-channel or tqb), or a top quark and a b quark (s- 
channel or tb) . To successfully identify single top quark production in the presence 
of backgrounds, the W — > ef(uf) decays are required. Thus single top quark events 
contain an isolated hight px lepton, significant Ifix due to the neutrino, two 6-jets 
(and a light quark jet in the t-channel production). 

In addition, for both the ti and single top quark production events, the heavy 
flavor hadrons from b quark fragmentation will have a long lifetime. This means that 
final state particles will originate from a location which is some distance away from 
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the event primary vertex. Sometimes these 6-jets will exhibit a soft, non-isolated e 
or fi from semileptonic b decay. 



2.4. Modeling top quark signal and background 

Measurements of the top quark require a firmly understood simulation of signal and 
background processes. This understanding includes a need for reliable estimates of 
event selection efficiency and the kinematic and angular distributions in the final 
state. Event generators are used to simulate the physics of production, decay and 
hadronization. Two kinds are employed. Exact leading-order calculations such as in 



Alpgen 1.2p2l result in the partonic final states at the origin of observed events. 
Alpgen in particular also includes the spins of the particles. Showering Monte Car- 
los such as Pythia 6.2^ or Herwig 6.4^ simulate the hard scatter with leading 
order elements and employ a phenomenological showering mechanism to account for 
hadronization. An approximate perturbative QCD description of gluon radiation is 
implemented in these simulations. Pythia also adds the underlying event from the 
proton antiproton interaction. The two types of generator are often used sequen- 
tially so that the hard and soft physics are dealt with as completely as possible. In 
more recent experimental analyses, a significant effort is being made to correctly 
account for potential double-counting of soft QCD processes when this is done. A 
prop osed m echanism to match observed jets with partons is now employed in these 
case d^^ . An alternative approach is also available^!. Specialized simulators are 
also employed to perform important calculations for r decays (Tauola^I) or b and 
c quark decays (QQ 9.ll^, EvtgeiJSSI) . These generators are reviewed in detail in 
Ref. 1661 The specifics of their use are documented in some of the papers reviewed 
below (e.g. Ref. EH). 

Various Monte Carlo generators were used to model single top quark events. D0 
uses the SINGLETOP package^, based on the CompHep Monte Carlo generator^ 
to produce the parton four-vectors of the single top quark signal events. This pack- 
age uses a leading order simulation for the s-channel and a next-to-leading-order 
simulation for the ^-channel processes. Spin information in both the production and 
decay is included. The matrix element event generator MadGraph /MadEveniEOI 
is used by the CDF experiment to simulate the signal processes. 

Most of the Run II CDF analyses described in this paper use these genera- 
tors in the following way. In general, the CTEQ5L structure functions are used. 
Monte Carlo samples for tt to calculate acceptances and efficiencies, and to un- 
derstand kinematic shapes, are generated with both Pythia and Herwig. Alp- 
gen plus Herwig is used with the parton matching algorithm to model the W+ 
heavy flavor backgrounds. Diboson backgrounds are also generally modeled with 
Alpgen+Herwig. Single top quark backgrounds for £+jets channels, and Z/j* 
backgrounds for dilepton channels are studied using Pythia. In general, b and c 



quark decays are handled by the QQ generator^. 

The D0 Run II simulation generally uses Alpgen fed into Pythia for signal 
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estimation. VF+jets and diboson production are also simulated with Alpgen and 
Pythia. Z — > tt backgrounds are studied using Pythia. In general, TAUOLA and 
Evtgen are used for r and b/c hadron decays, respectively. The Q 2 scale used by 
D0 single top quark s-channel and tt samples is mf , while the t-channel is generated 
at (m t /2) 2 . CTEQ6M (CTEQ6.1M) parton distribution functions are used for single 
top quark (tt) events. The components of the IF+jets process: heavy flavor process 
qq' — > Wg with g — > bb or g — > cc, and t/g — ► VFc) are included in their relative 
proportions estimated using Alpgen and normalized to the data. 

3. Experimental Facilities and Reconstruction 
3.1. Experiments 

The Tevatron provides sufficient energy to create the extremely massive top quark. 
The two experiments at the Tevatron, D0 and CDF, are therefore unique in their 
ability to study the top quark directly. The accelerator complex at Fermilab accel- 
erates protons through a linear accelerator and three synchrotrons, the Booster, the 
Main Injector, and the superconducting Tevatron ring^l The Main Injector was 
built for Run II to increase the number of antiprotons that can be injected into the 
Tevatron. Antiprotons are produced by extracting protons at intermediate energy 
onto a nickel target. They are selected from the secondaries that are produced and 
then cooled and accelerated before they are injected into the proton acceleration 
chain. Proton and antiproton beams are accelerated to a final energy of 980 GeV 
in the Tevatron to provide proton antiproton collisions at 1.96 TeV, up from 1.8 
TeV in Run I. The use of antiprotons is important because it provides the ability 
to use valence quark annihilation, with available higher parton center of mass ener- 
gies, for top quark production. Instantaneous luminosities are high, so far reaching 
2.6 x 10 32 /cra 2 /s. Typically around three proton-antiproton interactions occur with 
each beam crossing. 

Collisions occur in the two locations on the Tevatron Ring where the D0 and 
CDF experiments are located. Both detectors (see Fig. [3]) have magnetic central 
tracking regions which include silicon microstrip trackers for precise vertex mea- 
surements. Field strengths of 1.4 T and 2.0 T are achieved via superconducting 
solenoids for CDF and D0, respectively. The higher value for D0 partially compen- 
sates for a smaller tracking volume. Primary charged particle tracking is performed 
for CDF by a drift chamber which obtains 96 position measurements for each par- 
ticle. Acceptance ranges to \rj\ < 2.0 for silicon and \rj\ < 1.0 for the drift chamber. 
D0 utilizes a scintillating fiber tracker composed of sixteen doublet fiber layers. In 
tandem with the silicon microstrip d etector an acceptance of |ry| < 3.0 is achieved. 
All trackers replace Run I equivalents^^!, resulting in better 6-tagging efficiency 
for CDlffl and a new 6-tagging and track momentum measurement capability for 
Dgffl 

Sampling calorimeters with large acceptance provide energy measurements for 
electrons, photons, and hadrons. D0 utilizes uranium absorber bathed in a liquid 
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Fig. 3. The Run II D0 detector (left) and CDF detector(right). Each detector has an inner 
tracking region surrounded by calorimeters and outer muon spectrometers. 



argon sampling medium to achieve near-compensating performance, good energy 
resolution and low backgrounds for the surrounding muon spectrometer^. Coverage 
extends to \rj\ < 4.2. Scintillator-absorber layers provide these measurements for 



CDF in the region \rj\ < 3.8^21. New endplugs for CDF^S, and electronics and 
hardware trigger upgrades for D(zf^ deliver improved performance for the high 
instantaneous luminosities achieved in Run II. 

Drift tubes and scintillator arranged in l ayers o utside of the calorimeter are used 
for muon identification for both experiment J 58 | 59 J j n rj0 casej these layers are 
distributed in front of and behind thick background-suppressing magnetized iron 
toroids which provide a 1.8 T field for the full coverage to |^| < 2. The CDF 
detector has coverage to \rj\ < 1.0. The prim ary en hancements from Run I are 
improved forward coverage for both experiment J60 | 6 ^j 

Data are taken using a three-tiered trigger system composed of one hardware and 
two software levels. Luminosity is measured via gas Cerenkov or plastic scintillation 
counters mounted in the forward regions of CDF and D0, respectively. Analyses 
utilize data taken from 2002 until Spring 2006. By Spring of 2007, each experiment 
had accumulated about 2.5 fb _1 of collision data. 



3.2. Reconstruction 

The effective identification of e, fi, v, light quark jets, b jets and ts is crucial to the 
analyses reviewed in this paper. One cannot properly demonstrate the requisite un- 
derstanding of signal and background efficiency without a thorough understanding 
of the performance of the detector to measure their energies or momenta, their ef- 
ficiencies, and their rates to be incorrectly reconstructed. This subsection describes 
the methods of identification for these primary objects as well as the ways in which 
their performance is determined and propagated into experimental models of the 
processes being considered. 
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Electrons and Muons: Electrons are reconstructed as grouped or 'clustered' elec- 
tromagnetic calorimeter (EM) elements which exhibit significant energy. Cells are 
the fundamental elements used by D0. Three adjacent towers are used by CDF. 
These clusters are associated with a well-reconstructed charged particle track. A 
range of parameters are used to identify good electrons, including the ratio of EM 
to hadronic energy, isolation from more energy in a wider cone around the centroid, 
cluster energy divided by track momentum {E/p), and shower shape. CDF applies 
these cuts for 'tight' electrons, and omits the isolation requirement for 'loose' elec- 
trons. D0 applies basic cuts on electron candidates and then assembles a likelihood 
discriminant which combines these and other calorimeter and tracking parameters. 
For D0, only isolation and EM-fraction cuts are applied for loose electrons. Muons 
are identified from track signatures in the outer muon spectrometers of each ex- 
periment, which are then matched to charged particle tracks. In the case of D0, 
isolation cuts are applied based on tracks and calorimeter energy. Loose muons may 
pass a looser isolation selection. 

The efficiency, momentum resolution and energy scale performance of e and /j, 
in the detector are typically measured by examination of Z — > ee, fifi events. By 
requiring one 'tag' lepton to be well-identified, and employing other cuts to suppress 
non-Z backgrounds, the other 'probe' lepton can be measured in an unbiased way. 
The reconstructed mass of the Z boson is a direct reflection of the electron energy 
scale. Isolated leptons can be misidentified typically from jet activity. This can occur 
through fragmentation to 7r° — > 77+ track where the track is from a conversion or 
charged hadron overlap. Semileptonic decay of hadrons in the jet can also rarely 
provide isolated leptons with significant px- The probability at which jets mimic 
leptons is termed the 'fake rate'. It is often extracted from multijet samples where 
W and Z/j* contributions are suppressed. For instance, D0 uses such samples 
with leptons passing very loose requirements. It is important that signal trigger 
conditions on the faking lepton be reproduced. 

These selections are applied to find leptons from the W bosons in top quark 
decay. It is also important to be able to identify leptons from semileptonic heavy 
flavor meson decay. Such leptons can be used to flag a reconstructed jet as having 
come from the fragmentation of a 6— quark. These 'soft leptons' do not have to 
satisfy the isolation requirements. Estimation of the rate in data to tag b— jets with 
a soft lepton is usually ascertained from a dijet sample where one jet is tagged and 
a non-isolated lepton is sought in the other. The rate for typical jets, primarily light 
flavor jets, to produce a 'mistag' is measured from a multijet sample, such as 7+j'et 
events. 

Jets, rs and $t' An iterative, fixed-cone algorithm" is applied to electromag- 
netic and hadronic calorimeter elements to reconstruct jets. Because of the high jet 
multiplicity in top quark events, cone sizes used are generally small: for D0 (CDF) 
most analyses use AR = y/ Ac/) 2 + Arj 2 = 0.5(0.4). Jet energies are corrected for var- 
ious effects such as non-linearities in single particle response, unregistered energy 
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at detector element boundaries and cracks, particle showers, and underlying event 
and multiple interactions. The jet energy scale is ultimately measured using 7+ jets 
events where the photon in the electromagnetic calorimeter section is calibrated via 
the Z — > ee resonance. The general approach compares the momentum of the jet 
in the transverse plane ('transverse momentum', or 'pt') to the reference provided 
by the photon px, although in practice several effects mean that a direct px com- 
parison is not optimal^. A subset of jets can be further identified as originating 
from the fragmentation of heavy quarks, particularly b— quarks, if they encompass 
an associated, non-isolated low px electron or muon from semileptonic decay or if 
they match a secondary vertex. Hadronic r lepton decays are also identified as jets. 
In this case selections are applied based on shower profile in the calorimeter, as well 
as requirements of one or three tracks in the jet cone. 

Neutrinos escape the detectors undetected. Their presence can be inferred from 
a significant imbalance in the px of all observed particles in an event. Recall that, 
since most final state particles are essentially massless compared to the magnitudes 
of their momenta, p ~ E. So the event-wide pt imbalance is generally termed 
missing Ex, or $t- The $x is calculated as the negative of the vector sum of the 
Ex's of individual calorimeter elements. This quantity is then corrected for the pr's 
of any reconstructed /x's in the event, as well as the energy scale corrections for jets 
and electrons. 

Vertices: The pp collision point is called the event primary vertex. It is recon- 
structed with well-measured tracks that are consistent with origination from the 
same point. In many cases, more than one primary vertex are reconstructed. These 
may be due to extra pp interactions resulting from the high operating instantaneous 
luminosity. There will generally be only one vertex corresponding to the hard scatter 
giving rise to a top quark. 

The secondaries from the decay of long-lived B hadrons produce charged par- 
ticle tracks that originate from a vertex that is displaced from the primary vertex. 
Both CDF and D0 employ algorithms that identify ("tag") jets as originating from 
b— quarks by reconstructing such secondary vertices within a particular jet ('sec- 
ondary vertex tag'). Such secondary vertices are reconstructed from at least two 
well-reconstructed tracks with px > 1.0 GeV that match a jet. These tracks must 
have an impact parameter in the transverse (x, y) plane with a significance of > 3a 
with respect to the primary vertex for CDF. D0 applies a 3.5<r cut. CDF adds a 
second, three-track vertex with relaxed track px and impact parameter cuts. The 
vertex reconstructed from these tracks must be significantly displaced from the pri- 
mary vertex in the transverse plane, given the known tracking resolutions and it 
must be on the same side of the primary vertex as the jet, to tag a jet as a 6-jet. 

Other tagging algorithms exist. For instance, an 'impact parameter tag' algo- 
rithm has been used by D0 which is based on counting the number of tracks with 
impact parameter significance above a certain value. Both experiments have also 
implemented 'jet probability taggers'. These use the knowledge of track resolutions 
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to calculate a probability that tracks in jets originate from an event primary vertex 
and compute a jet probability based on the combined track probabilities. The jets 
are tagged if the jet probability has a low value, indicating an inconsistency with 
the hypothesis for that jet to originate from the primary vertex. D0 and CDF also 
use neural networks to combine the various parameters associated with a 6-tag such 
as impact parameter, momentum and invariant mass of all tracks associated with 
the vertex, vertex displacement, etc. into a more powerful discriminant. 

It is important to understand the efficiency for b— jets to be tagged, and light 
quark and gluon jets to be mistagged. Jets from c— quarks and their relevant tagging 
rates must also be accounted for. The efficiency for a particular b— tagging algorithm 
is generally measured in heavy-flavor enriched inclusive jet samples in data. Such an 
enrichment can be achieved by lepton-tagging a jet for which the b— tag efficiency is 
to be measured. Monte Carlos are used to obtain an estimate of the residual c— quark 
content and other selection biases. The mistag rate to falsely identify light quark or 
gluon jets as b— jets is estimated from samples of inclusive jet events. The primary 
backgrounds produce displaced vertices that are equally probable to be 'behind' the 
primary vertex (negative tags) from the point of view of the jet they correspond 
to, and 'in front of the primary vertex (positive tags). However, heavy flavor will 
produce an asymmetric distribution with more positive vertex displacements. The 
negative tag rate is a useful first approximation to the mistag rate. It must be 
corrected for heavy flavor in the jet samples, as well as for long-lived particles in 
actual light flavor jets that are not reflected in the negative tag rate. Final tag and 
mistag rates are parametrized as functions of px and rj. 



Performance and parametrization of simulation: The final element tying the 
simulation of signal and background sam ples to the data involves the simulation of 
the detector. For CDF and D0 , GEANT3^is used with a full detector simulation. 
At this stage minimum bias Monte Carlo events were added to the hard scatter such 
that the Poisson distribution mean matches the average instantaneous luminosity 
of the data sample. 

In general, the simulation does not replicate exactly the detailed performance 
parameters of the detector. The parameters of interest are: e and [i identification 
efficiency, momentum resolution and energy /momentum scale; jet efficiency, energy 
resolution and energy scale; $t resolution; and b— tag algorithm efficiency. As a 
result, it is necessary to quantify each element of performance in the simulation 
and in the actual detector using isolatable control samples. For leptons, these are 
generally the Z — > ee, fifi samples. For jets, an important sample is the 7+jet 
sample, b— tagging requires simulated inclusive jet samples. The simulation is then 
corrected so that the performance for real data is replicated. A few performance 
measures must be measured only for use in data. These are generally the lepton 
fake rates and the 6— tagging mistag rates. 
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4. Top Quark Pair Production Cross Section Measurements 

Measurement of the top quark pair (tt) production cross section permits a unique 
test of perturbative QCD predictions. Unveiling discrepancies between measured 
rates and expectations for various final state channels provides a potential indi- 
cation of new physics. The selection of clean, well-understood analysis channels 
facilitates the measurement of the top quark mass and other properties. In this sec- 
tion, we discuss the primary channels that have been employed for the cross section 
measurement. These consider variants of the dilepton, £+jets and all-jets channels. 
Analyses of r channels, while they have been pursued, have not yet provided signif- 
icant constraints on a t i and are omitted from this section. 

Two primary aspects of these analyses are general in their impact. The event 
selection that is used for all top quark measurements has been developed based 
on extensive studies of the expected properties of tt events as well as those from 
known backgrounds. The lepton, jet and $t cuts are derived from these studies, 
as are more sophisticated parameters which consider angular correlations, scalar 
Et sums, or event topology. Given a particular event selection, the behavior of 
signal and background events in terms of these variables must be understood well 
enough to keep uncertainties in the measurements under control. It is particularly 
important that the background models be validated wherever possible to optimize 
the analysis sensitivity. Achievement of an effective selection whose performance can 
be confirmed not only ensures an optimal cross section measurement, but it also 
serves as the foundation for measurements of the properties of the top quark. We 
discuss these two issues in this section and then review the individual measurements 
from D0 and CDF. 
4.1. Selection variables 

Events containing a tt pair have several distinctive characteristics that arc not 
reflected by the anticipated backgrounds. In general, the total scale of energy in 
the event, particular in jets, is quite large. Decay products are roughly isotropically 
distributed in top quark events. These two qualities follow largely from the high mass 
of the tt pair. A third category comes from the observation that some backgrounds 
come from mismeasurement and this reflects itself in angular correlation between 
objects. The discriminating variables used fall into the following categories: 
Energy scale: Due to the large mass of the top quark, the characteristic energy 
scale of the tt event is significantly larger than that of the average QCD background 
event. This means that tt events generally have more energetic jets and larger multi- 
jet invariant masses. This is especially true of the leading two jets from the b quarks 
which come directly from top quark decay. This total energy scale is most often 
expressed as a scalar sum of object Et(pt)s, termed Ht- Sometimes the discrimi- 
nating power can be increased by also adding observed lepton pxs and perhaps also 
the Et- Invariant mass parameters, such as for the leading two jets (M^ n ), the 
second leading two jets (M^ n ) or all of the objects (£, jets, Et), also indicate the 
production of a high mass state. 
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Event topology: In general, backgrounds with jets originating from gluon radia- 
tion provide steeply falling px spectra for jets. The manner in which color flows in 
such events tends to produce particles that congregate in a planar geometry. Addi- 
tionally QCD events usually have a more back-to-back spectrum because of their 
hard scatter origin. The jets from ti decay on the other hand, are almost isotropi- 
cally distributed. These differences can be quantified using event-shape parameters, 
such as sp her icity, iS, and aplanarity, A, calculated from the normalized momen- 
tum tensoJ^. Top quark events tend to have higher A and S than the background 
events. 

Rapidity spectrum: ti events on average are expected to have more jets of higher 
energy and with less boost in the beam direction, resulting in events with many cen- 
tral jets. The QCD background tends to have jets that are more forward-backward 
in rapidity. One variable considered is the centrality (C), defined as C = ^Ji T , 

where is the total energy in the event (or it may represent the scalar sum of the 
momenta (energies) of jets Y,Ei ets ). Another variable is < rj 2 >, the pr-weighted 
mean square of the 77 of the jets. The QCD multijet background is expected to have 
a broader distribution in < if > than the ti signal. 

Angular Correlations: Several of the backgrounds arise because the Et(pt) of 
one or more objects in candidate events have been significantly mismeasured. Often, 
this mismeasurement will create the appearence of $t either in the azimuthal direc- 
tion of the offending object, or opposite it. Jet energy and \x momentum resolutions 
are most often the cause. As a result, cuts are placed on the angular correlation of 
objects: Acj)(x,yi) is the azimuthal angle between the transverse momentum vector 
of objects x and yi. Here, the subscript i denotes leading, second leading, etc. of 
object y. Common uses are A4>(J£t, Mi) or A(j)(^T,j) where /ii is the leading muon 
or j is a jet. 

4.2. Background extraction from data 

Production of top quark pairs presents a unique set of final states. However, the high 
mass of the top quark also means that the production cross section is quite small. 
Use of the variables discussed in Section ^. ll cntail a heavy reliance on measurements 
of jets and $t- These are complex observables defined by sophisticated algorithms 
and subtle physics at a hadron collider. Simulation of background processes must be 
complemented by study of those processes in the data. Often, the attempt to validate 
the background model will result in corrections to the simulation that will improve 
agreement. In other cases, the background will be taken from the data alone or in 
combination with the simulation. The simulation can still provide a valuable cross- 
check of the model, since it represents the best attempt to estimate the background 
from first principles. This in turn gives confidence in the signal model which can only 
be obtained from Monte Carlo. Common schemes to be elaborated in the following 
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sections are: 



Modeling Jet Production: The physics behind jets involves soft processes and 
the algorithms used to reconstruct jets have complex behavior. For all analyses, it 
is therefore important to demonstrate an understanding of the background (and 
signal) at all jet multiplicities. In low jet multiplicities, background will dominate. 
Agreement between models and data provides a confirmation of the background 
model in the signal rich bins. Typically this will be reflected as an estimate of the 
sum of signal and background ('5 + £?') vs. jet multiplicity, with the observed data 
in comparison. Another approach is to consider the agreement of the modeled and 
observed jet ptS- This has occasionally been used, for instance, to determine the 
background level from a fit to the leading jet pt- 

Lepton Backgrounds: Jets also impact analyses because they can mimic charged 
leptons or evidence of ^s. In the former case, it is difficult to precisely model those 
jets that shower in the detector to look like electrons, or that produce a lepton 
from semileptonic decay which can pass lepton isolation and px cuts. The rate for 
jets to fake leptons is obtained in data as mentioned in Section [3. 21 Two things are 
crucial in the use of these rates. First, samples must be identified in data which 
have negligible signal contribution. Second, it is important that the faking, which 
is a result of hadronization or decay, be uncorrelated with the kinematics of the 
event. Typically, samples are selected with all kinematic selections and one lepton 
identification relaxed. In the case of $t, the rate to satisfy certain $t cuts in a 
f-less process (termed the '$t fake rate') is studied in Z — > ££(£ = e,/x). 

Jet Flavor Modeling: Identification of b— jets from evidence of displaced sec- 
ondary vertices or non-isolated soft leptons plays an important role in selecting a 
sample enriched with top quark events. In all 6— tag analyses, the rate for light 
quark or gluon jets to accidentally satisfy the soft lepton or displaced vertex crite- 
ria is called the 'mistag rate' as described in Section 13.21 The application of these 
mistag rates to an untagged background sample can be somewhat more compli- 
cated than for the leptons. The reason is that jets arising from the hadronization 
of c— quarks can also harbor high impact parameter tracks. Also, all backgrounds 
contain some relative fraction of light quark, c and b jets. In VF+jets production, 
most of the associated jets arise from light quarks or gluons. Occasionally, gluon 
splitting will generate a bb or cc pair which give rise to heavy flavor jets. More 
rarely, sea s— quarks radiate a W boson giving a single c— quark and subsequent jet. 
It is common in this instance to use Monte Carlo to provide some estimate of the 
relative flavor content, and then to apply the mistag rates or tagging efficiencies to 
the relevant background fractions. 
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4.3. Dilepton final states 

The signature of two isolated, high px leptons and $t in association with two 
high Et jets is a striking consequence of the tt quark decay chain where both 
W bosons decay leptonically. The only processes exhibiting the leptonic signature 
are diboson production, particularly WW where both W& decay leptonically as 
in the ti case, or Z — ► tt production where both r's decay leptonically. Diboson 
cross sections are of the same order as that expected for tt production. Z — ► tt 
production is much greater but the dilepton branching fraction is small and the 
charged lepton and v pr's are soft. In both cases, additional jets are only produced 
at a rate of approximately a s each and so jets become a primary discriminator of 
these backgrounds. 

Beyond these processes, it is possible to mimic the top quark signature through 
instrumental mismeasurement. The chief background arising this way comes from 
Z/7* — ► ee, pp+ jets production. Here the reconstructed $t is the result of tails 
in the jet or lepton energy or momentum resolutions. An important background 
also comes from W — ► e, p+ jets production where one of the jets fakes an iso- 
lated lepton. QCD multijet production is a background when both leptonic and 1£t 
mismeasurement occur. 

In order to extract a measurement of the top quark production cross section 
in the dilepton channels, different techniques have been developed. One of the 
most widely used has been a selection where loose and tight versions of the lep- 
ton identification are employed. For 'explicit' dileptons, these distinctions still refer 
to selections requiring convincing signals in the subdetectors responsible for lep- 
ton identification, particularly the electromagnetic layers of the calorimeter and the 
outer muon spectrometers. However, these detectors have significant gaps and holes 
in their coverage, and lepton identification and isolation requirements are ineffi- 
cient. The resulting loss in acceptance seriously reduces the signal event yield in the 
dilepton channel. In the most extreme 'implicit' dilepton strategies, identification 
requirements are relaxed to mere requirements of an isolated, high pt track or even 
to no evidence of a second lepton at all. Such an approach can accept dilepton events 
where one of the leptons is a r which decays hadronically. Event-wide kinematic 
variables and tagging of jets from b— quarks are valuable ways to isolate the dilepton 
final state, particularly when the lepton requirements are relaxed and backgrounds 
are higher. 

4.3.1. Explicit channels 

Explicit dilepton analyses are defined as those in which both leptons are fully recon- 
structed as either eor/x. W — > r — > e, p decay modes cannot be distinguished from, 
and so are grouped with, the direct W — > e or W — > p channels. Three channels 
have been examined by both Tevatron experiments: ee, ep and mi. 

The D0 experiment has used explicit channels to p roduce a cross section mea- 
surement at yfs = 1.96 TeV with 243 pb _1 of data Events were triggered with 
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two leptons in the first level hardware trigger and one or two leptons in the higher 
level software triggers. In general, the Run II measurements derive many techniques 
from the Run I analyses-^ (see final combined Run I cross section in Table [5]) . Al- 
though the differences are less than in Run I, momentum resolutions are significantly 
different for e's and /i's. This means event selections across the three channels are 
somewhat varied. Also, because 6-jets carry away a large momentum directly from 
top quark decay, D0 uses relatively soft lepton pt cuts and stiff jet Et cuts. The 
primary background which arises after these cuts is mainly Z/7*+jets production. 
A substantial cut on $t reduces this background considerably, but it is still the 
primary background in the like- flavored channels. Suppression of this background 
involves rejecting events with dielectron or dimuon invariant mass (M ee or M^) 
consistent with that measured for the Z boson mJ^. In the ee channel, the ex- 
cellent electron energy resolution allows a narrow window to be drawn around the 
reconstructed M ee . The cut on ^t is elevated for events within this window. For the 
/i/j, channel, the muon momentum resolution from the tracker degrades at higher 
Pt and has significant non-gaussian tails. Instead of a window cut, a kinematic fit 
is employed which quantifies a \ 2 test with respect to Mz, Xai z - I n Run I, this test 
used the muon and calorimeter measurements and knowledge of their resolutions. 
The improved Run II detector has permitted a simpler test based on the muon 
momentum measurements alone. 

The selections were optimized to produce the smallest fractional statistical un- 
certainty on the cross section, taken to be proportional to S/ yS + B. These opti- 
mizations used several variables, including event shape variables, dilepton mass and 
lf>T cuts. The most important of these variables have been the jet Et and Ht- The 
final selection criteria are given in Table [2 

Instrumental backgrounds consist of Z — > ee, background events where the 
Et may be mismeasured, and VF+jets and multijet events where the jets fake lep- 
tons. These are generally determined from data. Lepton backgrounds are estimated 
by folding fake rates into a signal-like sample with one tight lepton. For the Z — > ee 
background, the rate for $t to pass the D0 cuts is assessed in a sample of 7 + 2 
jets events that are kinematically similar to the Z + 2 jets sample. Agreement in 
the IfiT spectrum is observed for the Z+jets data sample confined to the region 80 
GeV< M ee < 100 GeV and a high statistics 7+jets sample. The 7 + 2 jets behavior 
also agrees well with a high statistics Z{— * ee) + 2 jets Monte Carlo sample where 
all known effects on the $t resolution in data have been accounted for. The ratio of 
events to pass divided by those that fail the $t cut is extracted from the data and 
multiplied by the number of ee events failing the $t cut in data. The Z — > 
jets background is estimated by extracting the $t fake rate from fully simulated 
Alpgen Monte Carlo and multiplying that by the number of observed /i/i events 
in data. 

Physics backgrounds are those in which the full signature from the top quark 
pair is mimicked at the particle level. These are Z — > tt — * 11 + 2i/s, WW — > 2^ + 2^ 
and WZ — > 3£+v. The former is the dominant background for the efj, channel. These 
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backgrounds were estimated from Monte Carlo corrected for efficiencies measured 
in data. 




Jet multiplicity Leading Lepton pj (GeV) 

Fig. 4. Kinematic characteristics of dilepton events in 243 pb — 1 of D0 dataES Event jet multi- 
plicity (left) and leading lepton pj- (right) are shown. Data are indicated by points and error bars, 
and the sums of signal and background distributions are overlaid as solid histograms. 

The cross section for y/s = 1 .96 TeV is taken by maximizing the likelihood among 
the three channels that, given a value of the assumed cross section, the estimated 
backgrounds and top quark efficiencies can produce the number of observed events. 
The cross section is 8. 6+2 7 (siat) ± l.l(sys) ± 0.6(/wm) pb. Fig. 2] indicates the jet 
multiplicity of the dilepton sample with all other selections applied. The excess of 
events over background in the 2-jet bin is consistent with the expected top quark 
contribution. Fig. U also provides the px distribution for the leading lepton in the 
final event sample. Kinematic distributions of dilepton events have been of some 
interest. The CDF collaboration has observed some unexpected properties of their 
dilepton samples in Run I and Run ii^. The D0 distributions are consistent with 
standard model effects. The dominant uncertainties for the D0 measurement are 
roughly equal from lepton identification and triggering, and from jet reconstruction 
and energy scale. 

This basic analysis has been performed in 425 pb -1 of data with somewhat modi- 
fied selection on the e/i and channels. In the latter case, substantial improvement 
was made in the rejection of Z bosons. Figure [5] shows the Ht distribution for the 
e[i data, and for tt and background expectations. Good agreement is observed and 
the statistical sample is starting to permit a fairly distinct tt component to become 
more evident. In combination with the £+track channels described in the next sec- 
tion, a combined cross section measurement of 7.4 ± lA(stat) ± l.O(sys) pb was 
obtained ^1 A preliminary result in I fb _1 from the dilepton and ^+track channels 
has given 6.2 ± 0.9(stat)±Q%(sys) ± 0A(l urn) pb These results are given in Table 

El 

For Run II, CDF has analyzed an explicit dilepton data sample from 197 ± 12 
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Fig- 5. Kinematic characteristics of dilcpton and £+track events in D0 analysis of 425 pb -1 
sample uJA, Event Ht distributions are shown for the efj, channel (left) and the combined £+track 
channel (right). Data are indicated by points with error bars, the estimated signal and background 
contributions are shown via overlaid histograms. 



pb-iLL^ Triggering involved one high px lepton, e or /x. In general, the CDF analyses 
are characterized by a common approach across channels which aids in background 
understanding and justification of event selection. The selection in Run II is similar 
to that used in Run I E01 ( see combined measurement in Table [5]). Top quark events 
were selected by requiring two leptons with px > 20 GeV, plus two energetic jets 
with pt > 15 GeV and large |?t(>25 GeV). One of the two leptons must satisfy 
a tight selection, while the other may pass looser cuts. For both the ee and fj.fi 
channels, events are rejected if the dilepton invariant mass, My, is near the Z 
boson mass, Mz, to suppress Z backgrounds. The basic kinematic selection is given 
in Table [21 An opposite sign requirement for the two leptons is applied. Events 
were removed when A0($r,j) is small and when A(f>($T,£) < 20deg. CDF did not 
require that explicit and implicit channels have orthogonal event selection. 

Remaining backgrounds from Z/j* were estimated by comparing the number of 
dilepton events in each jet multiplicity in the data and using that to normalize a 
fully simulated Pythia sample. Lepton instrumental backgrounds were estimated 
by folding fake rates measured in data with a sample selected to be kinematically 
signal-like, but without lepton identification selection on one of the leptons. Diboson 
backgrounds were estimated from events simulated with Alpgen plus Herwig or 
Pythia. Fig. [6] shows the $t and Ht distributions of the dilepton events after the 
application of all the selections. 

The Run II explicit dilepton channels were incorporated into a ti cross section 
measurement from all dilepton channels, as shown in Table El In this channel, 13 
events were observed with expected signal of 8.2 events and expected background of 
2.7 events. Since the selections overlap, explicit channels have significant correlation 
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Fig. 6. Kinematic characteristics of dilepton events for CDF analysis in 197 pb -1 . Th e T^t for 
explicit dilepton events (left) and the event Ht for £+track events (right) are shown UQi. Data are 
indicated by points with error bars, and estimated signal and background contributions are shown 
via overlaid histograms. 

with the implicit channels described in Section f4. 3. 21 The combined dilepton cross 
section for Run II is provided in Table [5l The cross section measurement was re- 
peated for subsamples that have at least one b— tag, and that have two-tight leptons, 
and were found to b e consistent with these results. Preliminary analyses of 360 pb -1 
ISHand 1.2 fb _1 ^21 have been pursued by CDF. The former combined both dilepton 
and £+track measurements to obtain a cross-section of S.Sj^'^siai^Q'^sys) pb- 
The latter resulted in a measurement of 6.2 ± l.l(stat) ± 0.7 (sys) ± 0.4 pb. 

4.3.2. Implicit channels 

Analysis of a dilepton channel where one lepton is not identified began with 
an initial effort by D0 in Run which looked for one isolated, high pt electron 
plus exceptionally high $t- For a concise review of this analysis, see Ref. |3J A cross 
section for this channel alone was obtained to be 9.1 ± 7.2 pb, which is incorporated 
into the final D0 Run I result given in Table [5] 

In Run II, CDF has performed a c ross section measurement in 197 ± 12 pb -1 
using an implicit dilepton selection^! A tight lepton (e or fx) was required with 
strict calorimeter, muon system or tracking requirements, and an isolation cut. 
Another high px track isolated from significant momentum in nearby tracks was 
also required. This lepton plus track (£+track) approach provides measurement of 
all final state particles, aside from vs. Thus, top quark properties, such as m t , can be 
measured with techniques appropriate to an explicit dilepton analysis. The cost of 
the looser lepton selection can be borne by tighter requirements on other quantities, 
such as jet Et or jet b— tagging. 

As with the explicit channels, Z/~/* backgrounds were estimated by comparing 
the number of dilepton events in each jet multiplicity in the data and using that to 
normalize a Pythia plus full detector simulation sample. Instrumental backgrounds 
were extracted from a kinematically signal-like sample in data without final identi- 
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Table 2. Selection cuts for D0 and CDF Run 

II dilepton and f+track cross section measurements. Variables are described in the text. 



cut 


CDF dilepton 


CDF £+track 


D0 dilepton 


D0 £+track 




2 


1 (1) 


2 


1 (1) 


Pt 2 (GeV) 


> 20 


> 20 


> 15 


> 15 


l»fcl 


< 2.0 


< 2.0 


< 2.5 


< 2.5 


\Vfi,track | 


< 1.0 


< 1.0 


< 2.0 


< 2.0 


Nj e t s (Ntags) 


>2(0) 


>2(0) 


>2(0) 


>1(1) 


P T (GeV) 


> 15 


> 20 


> 20 


> 20 


M 


< 2.5 


< 2.0 


< 2.5 


< 2.5 


$ T (GeV) 


> 25 


> 25 


> 25 


> 15 (e) 

> 25 (/i) 




7^0 


7^0 






A(f>($ T ,e/track) > 


20 degw/l 


5 deg to / track 


0.2 w/ii(en) 
not 180 deg(/i/i) 




Ht (GeV) 


> 200 




> 140 (en) 




Z rejection: 










M u (GeV) 


76 - 106 


76 - 106 


80 - 100(ee) 


70 - 100(e) 
70 - lW(p) 








> 2(MM) 




near Mz 


CO 


> 40 


CO 


> 20 (e) 

> 35 (p) 



fication cuts on one of the lepton candidates. Diboson backgrounds were estimated 
using Alpgen with Herwig or Pythia. The tt signal was simulated using Pythia. 

The data yield 19 events for the £+track selection with expected signal of 11.5 
events and an expected background of 6.9 events. The implicit and explicit channels 
share a significant number of events. The combined dilepton cross section measure- 
ment is given in Tabled The cross section measurement was repeated for b— tagged, 
and two-tight lepton subsamples and found to be consistent with these results. Pre- 
liminary analyses in 360 pb" 1 and 1 fb _1 of data were performed for the £+track 
channel. The former yielded a combined measurement with dilepton channelJ^ as 
shown in Table [5] The lfb -1 analysis p rovi ded a much more precise measurement 
of 8.3 ± 1.3(stat) ± 0.7 (sys) ± 0.5(lum ) ptM Fi gure[7]shows the agreement between 
the jet multiplicity and $t distributions for the lfb -1 £+track data and the signal 
and background expectations. 

The D0 experiment has performed an £+track analysis in 425 pb~ x of datsP^. 
At trigger level, single lepton plus jet triggers are employed. To offset the higher 
background resulting from the omission of the lepton identification cuts, b— tagging 
of jets using a secondary vertex algorithm is employed to produce a reasonable S/B. 
This allows the analysis to have significantly looser kinematic selection than D0s 
explicit dilepton channels. Signal and Z/-f* samples are modeled using Alpgen fed 
into Pythia. Diboson samples are modeled with Pythia. The primary background 
is instrumental from Z — > ee, fifi with fake ]£t- The data and Monte Carlo $t dis- 
tributions were observed to be in agreement, and so this background is extracted 
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Fig. 7. Distributions of jet multiplicity in CDF £+track events in the 1 fb^ 1 sample (left). The 
distribution is shown for two jet events (right). The data is indicated by points with errors, 
and signal and background expecations are overlaid. 1221 

from the simulation after normalizing the event yield to that observed in the data 
for low ^t- A significant background from fake leptons or fake tracks exists for 
this analysis. This background is estimated in the untagged sample by construct- 
ing four samples in data with loose and tight selections on the lepton and track. 
From these samples, 6— tag rates appropriate to ly+jets events, and knowledge of 
lepton efficiencies and fake rates, the number of lepton instrumental backgrounds is 
inferred. Figure [5] shows the Ht distribution of ^+track events, with ti signal and 
background estimates overlaid. There is good agreement of these expectations with 
the dat a. The £+track channel was also analyzed in a preliminary measurement in 1 
cross sections determined in combination with the explicit dilepton 
channels are shown in Table El 
4.4. Single lepton channels 

Since top quark analyses have been statistically limited until recently, the extra rate 
has caused the £+jets channel to be key in studying the entire top quark sector. 
Generally, backgrounds are easier to control if the lepton is an e or fi, although r- 
based analyses have been implemented (e.g. Ref . [84)) . This section will concentrate 
on analyses using leptons of the first two generations. 

With four quarks from ti decay, top quark events are fairly crowded in the 
central region of the detector. Initial state gluons can land near these quarks. Final 
state gluon radiation can take momentum away from the quarks. The resulting jets 
may overlap and get merged, or may split into extra jets. As a result, a substantial 
fraction of top quark events will exhibit only three jets, and a significant number will 
have more than four. The physics processes which can produce the <?+jets signature 
are W+jets production in association with jets and, at a much lower level, diboson 
{WW, WZ) production. Single top quark production is also a background to the 
£+jets ft cross section measurement. The largest instrumental background comes 
from multijet production where the lepton and $t are fake. 
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Strategies for measuring top quark pair production in the £+jets channel have 
employed three general techniques. At the selection level, the balance of efficiency 
and background level is crucial. A purely kinematic, or 'topological', approach uses 
the unique kinematic signature of the top quark to isolate it from background via a 
multiparameter discriminant. Semileptonic decays of 6-quarks occur at a significant 
rate. The backgrounds do not have nearly as high a &-jet content as top quark 
events. So a second strategy involves use of the excellent capabilities of the D0 and 
CDF detectors to tag jets with leptons from semileptonic decay of b— quarks. The 
third strategy, which is the primary approach for the most precise measurements of 
the tt cross section, involves lifetime tagging b— jets to suppress backgrounds. 



Table 3. Selection cuts for CDF 
and D0 Run II €+jets cross section measurements. 



cut 


CDF I + jets 


D0 I + jets 


Trigger: 


1 e, fi 


1 e, fi + 1 jet 


b— tag selection: 






P l T (GeV) 


> 20 


> 20 


\Ve\ 


< 1.1 


< 1.1 


M 


< 1.0 


< 2.0 


Njets 


> 3 


> 3 


P J T (GeV) 


> 15 


> 15 


M 


< 2.0 


< 2.5 


0t (GeV) 


> 20 


> 20 


H T (GeV) 


> 200 


> 


topological: 






Njets 


> 3 


> 4 


H T (GeV) 


> 


> 


A0(#r,jl) 


0, TT 





4.4.1. £+jets Channels using kinematic selection 

The kinematics of the final state products from top quark decay are quite striking 
compared to those produced by other standard model processes. It is therefore 
natural that CDF and D0 have devised selections relying solely on these properties. 
In the £+jets channel, the backgrounds are much higher than the dilepton channel. 
As a result, significantly more complex approaches have been required to isolate a 
clean signal. A value of such a selection is that it does not rely on the assumption 
that BR(t — > Wb) — 1.0. Sensitivity is retained for models which produce different 
final state jet flavor content (e.g. Ref . l85|) . 

The initial untagged ^+jets analyses were performed by D0 in Run I and con- 
tributed the primary significance to their top quark dis covery analysis^. This laid 
the basis for a final measurement at yfs = 1.8 TeV^, as well as a measurement 
with 226 pb _1 of data at ^/s = 1.96 TeV^3. To control backgrounds, four jets were 



Review of Top Quark Physics Results 27 




Fig. 8. Distribution of Ht for combined £+jets events in 230 pb" 1 of D0 Run II data (left). The 
sum of ft and backgrounds are overlaid. At right is the multi-parameter likelihood discriminant 
for data and also overlaid with signal and background expect atiorl^l 



required. The specific event selection used in the Run II analysis can be obtained 
from Table [31 The ti acceptance and efficiencies were estimated from the Monte 
Carlo. 

Theoretical calculations yield substantial uncertainties in the production of back- 
ground processes. As a result, methods are needed by which the normalization of 
the background can be calculated from data. The multijet background kinematic 
shapes were taken from the data by requiring the lepton to fail the tight identi- 
fication requirement. To estimate the normalization of the background and signal 
estimates, two steps were employed. First, loose and tight cuts for the lepton identi- 
fication were defined to permit a variation of the level of QCD multijet background 
relative to the W processes (VF+jets and ti). Efficiences for real leptons and the 
fake leptons in the QCD sample were measured in data. Then one can solve the re- 
sulting two event yield equations in two unknowns: N s for top quark plus W events, 
and Nb for the multijet yield. This is termed the 'matrix method'. In Run I, the 
apparent scaling of jet multiplicity (sometimes called 'Berends scaling' ^23) was U sed 
to anchor the W+ four jet multiplicity normalization with the yields observed for 
background-rich low jet multiplicities. The Run II measurement, however, extracts 
the level of multijet background from the matrix method and obtains the W nor- 
malization as described in the next paragraph. The W+jets kinematic shapes are 
taken from the Monte Carlo. 

A multiparameter discriminant was constructed to identify the top quark signal 
in the data. For the Run II analysis, this was formed from six observables: Ht, 
A(p(J£ T J), C, S, A and K Tmm = ARJj in p™ m / . The latter is calculated by 
determining the jet pair with minimum separation in r\ — (f> (Ai?™™), the pt of the 
second leading jet of that pair, and the scalar sum of the lepton and E T , E™ . The 
discriminant function is: 

D = [IliSiix^/biixiMniSiixJ/biixi) + 1] (2) 

where Sj and bi are the normalized distribution for each variable, i, for signal and 
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background, respectively. An assumption was made of uncorrelated variables. Figure 
[5] shows the event Ht and the likelihood discriminant for data with tt signal and 
background superimposed. The discriminant function in data was fit to extract 
a measurement of <j+t and the number of W backgrounds. The measurement for 
y/s = 1.96 TeV was a t t = §.7±\^(stat)± l {X(sys) ± 0A(lum) pb 87 . This analysis 
has been updated in 425 pb -1 to yield a t i = 6At\'%(stab) ± Q.l(sys) ± 0A(lum) 
p^[88] preliminary result using 900 pb -1 has also been obtained ^1 Statistical 
uncertainties of approximately 15% are being achieved, and the measurement in the 
latter sample was um) pb. 
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H^GeV) ANN output 

Fig. 9. Distribution of Ht for combined £+jets events in CDF data (left). The sum of tt and 
backgrounds are overlaid. At right is the likelihood for data and also overlaid with signal and 
background expectationSSJ 



CDF pursued this 'topological' approach with 194 pb -1 of collisions at s/s = 1.96 
TeV^Sl Only three jets were required. Events were rejected if A(f>(^T, ji) ~ or 
7r. The VF+jets background shape is extracted from the Monte Carlo, and the 
normalization is obtained by maximizing a binned likelihood: 

L(n th n w , n q ) = U^r (e"^)/*! (3) 

where the n's are the means for top quark in+f), W (nw) and QCD (n q ) yields, 
and n, expresses the expected number of events in an ith bin given the probabilities 
for the signal and background contributions to populate that bin. di is the number 
of events observed in the ith bin. The value of n q is fixed based on a 'sideband' 
method where the data are divided into quadrants in the lepton isolation versus 
$r plane. The number of events with non-isolated leptons and high Ifix is scaled by 
the ratio of isolated-to-non-isolated events in a W-poor, low $t region. The value 
of n t i extracted from the fit is used to determine <r t j = n t i/ e t iC. The cross section 
is calculated using two different kinematic discriminants: Ht and a variable from a 
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seven-parameter artificial neural network technique. The parameters were: Ht, A, 
1/C, the minimum dijet mass from the three leading jets M™ n , the pseudorapidity 
of the leading jet if , the minimum r\ — <p separation of two jets among the leading 
three, and the Ht calculated from the third leading and lower pt jets. The values 
of both of the discriminants for the data sample, with signal and background su- 
perimposed, are shown in Fig. [5J The neural network approach included additional 
information besides Ht to provide approximately 30% better statistical uncertainty 
than Ht alone. These variables were each cross checked in the VF+jets samples. The 
cross section was actually estimated using these and other kinematic variables indi- 
vidually (see Fig. [T0|) . and a full range of systematic uncertainties were obtained for 
the Ht case in the three-jet and four-jet inclusive bins. The full analysis provided 
a measured cross section of 6.6 ± l.l(stat) ± 1.5(sys) pb. 
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£+jets Channels using soft lepton tagging 

Hadrons containing 6— quarks undergo semilcptonic decay at a rate of approximately 
17% for each charged lepton type. When decays through c— quarks are included, 
approximately 40% of tt events have a soft non- isolated /j,, for instance. In contrast, 
typical background processes produce primarily 'light flavor' partons (u/d/s quarks 
and gluons) and the heavy flavor contribution is small. So requiring a soft lepton in a 
jet can provide a strong background suppression. While the unbiased flavor selection 
is important in its generality, the ability to demonstrate the presence of b— quarks 
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in the top quark candidate sample has also been crucial to validate whether the 
observed signal adhered to standard model expectations. 
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Fig. 11. Jet multiplicity distribution for CDF 194 pb -1 soft-lepton tagged data sample (left) 
1221 Estimated backgrounds are overlaid as a solid histogram. At right is the D0 distribution in a 
preliminary analysis of 425 pb - 1 of dataEl. 



Both Tevatron experiments have pursued strategies of selecting events where at 
least one jet is 'tagged' with a soft, non-isolated lepton. CDF has produced estimates 
of cj££ using both soft muon and soft electron semileptonic decay modes^M In Run 
II, the soft muon variant was employed in 194 pb -1 of data^l In order to optimize 
their analysis, CDF chose a cut on Ht > 200 GeV that maximizes the significance 
of signal = S/VS + B. 

Signal efficiencies were estimated in Monte Carlo with tag efficiencies scaled 
from data. In the Run I analysis, efficiencies for leptons from heavy flavor, as well 
as backgrounds, were estimated using inclusive jet events from data. Cross checks 
were performed in a sample of Z — > It + jets events with soft lepton and secondary 
vertex tags which have a much smaller fraction of top quark signal. The calculated 
yield from background in these channels agreed well with observed numbers. In 
Run II, the W and QCD backgrounds are estimated from data in two steps. The 
multijet background is estimated with the sideband method mentioned above. The 
product of this yield with the correct tag rate gives the background estimate. The 
applicability of the measured tag rate for the W component of the background is 
verified by using it to estimate the tagged yield in several different jet samples. 
By comparing this estimate with the observed tagged event statistics, the rates are 
validated. The final W+jets yield is estimated using untagged samples corrected 
for the expected QCD contribution and multiplied by the mistag rate from 7+jet 
events. Other backgrounds are much smaller than the W and QCD backgrounds. 

Dominant systematic uncertainties arise from the modeling of top quark accep- 
tance as well as the signal and background tag rate estimates. The jet multiplicity 
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distribution for the data, as well as the estimated backgrounds are shown in Fig. QT] 
The final cross section estimate for y 7 ! = 1.96 TeV is 5.3 ± 3.3(stat)~^\^(sys) pb. 
A preliminary measurement of the ct+t in 760 pb -1 of /i— tag events has also been 

executed by CDF, yielding a value of 7.8 ± 1.7(stat)±\^(sys) pb^. 

Soft muon tagging has also been utilized by D0 to identify b— jets in £+jets data 
from both Tevatron runs. The Run II analysis^! constitutes a preliminary result 
in 425 pb -1 of collider data. These analyses benefited from the large coverage of 
the muon system (I77I < 2.0 in Run II). They also benefited from the depth in 
interaction lengths of the material in front of the muon chambers, particularly the 
layers outside of the muon iron toroid, which reduced the instrumental fake muon 
background to exceptionally low levels. Only three jets were required. The Run II 
analysis improves over the earlier one by no longer requiring an Ht cut or a selection 
for the fi + jets/ fi — tag channel of a cut on the x 2 of a kinematic fit to a Z — > 
hypothesis. 

Signal and physics background (VF+jets, diboson and single top quark events) 
efficiencies and yields were estimated using Alpgen plus CTEQ5L and Pythia. 
The primary backgrounds are: Z — > where one [i is lost or overlays a jet and mim- 
ics a b— tag, QCD multijet background, and VF+jets. The former was obtained from 
Monte Carlo but the yield was normalized using the measured Z boson cross section. 
The multijet background in the tagged sample was isolated from the W— like events, 
which include IF+jets and tt components, using the 'matrix method'. VF+jets events 
were estimated by extracting the yield from the same method in the untagged sam- 
ple. The Monte Carlo was used to ascertain the flavor composition of this sample. 
This is an improvement over the Run I approach in which the event kinematic 
shapes were also extracted from the Vecbos Monte Carlo. The appropriate tag- 
ging efficiencies or mistag rates for each flavor component were corrected to reflect 
performance in data and applied to the untagged yield. The jet multiplicity distri- 
bution of the selected sample is shown in Fig. [TT1 

The cross section was determined from a maximum likelihood fit to the observed 
number of events in three and four jet samples of the e+jets and /z+jets channels. 
In each subsample, the backgrounds were constrained by the values from the matrix 
method. A gaussian term was used for each systematic uncertainty such that the 
mean cross section could shift if the expected signal and backgrounds were varied 
over the ranges allowed by the uncertainties. The cross section was measured to be 
7.3tf!£(stat + sys) ± QA(lum) pb. 

£+jets Channels using secondary vertex tagging 

Displaced vertex 6— tagging as applied to top quark candidate samples was first 
performed by CDF in the analysis of 1.8 TeV collisions^. An initial cross section 
measurement from CDF using Run II dat a 1 95 1 96 1 employed a £+jets selection with at 
least one jet tagged with a secondary vertex using 162± 10 pb -1 . The leading jet pr 
was used as a discriminating variable in a likelihood to extract a top quark fraction 
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in data. The shape of the W boson px spectrum was obtained by selecting a VF+jets 
sample in data where no jets were b— tagged to reduce top quark contamination. The 
instrumental background was obtained in the high $t sample by reversing the iso- 
lation cut on the lepton. The kinematic distributions were observed to be insensitive 
to the heavy flavor composition of the sample. To validate the background model, 
untagged W+ one or two jets events were scaled by the tag rate. Agreement with 
the 6-tagged sample in data was observed in the jet px distribution. An unbinned 
likelihood fit was performed to the observed jet pt distribution given the expected 
signal and background shapes. A signal fraction of Rfu = 0.68^°'^ was obtained. 
The top quark cross section was obtained from a(tt) = N b s Rfi t /A^e^TL where A 
is the signal acceptance, e+£ is the estimated efficiency, and N b s is the observed 
yield of events in data. The systematic uncertainty for this analysis is dominated by 
that from the jet energy scale calibration. Other significant contributors come from 
the b— tagging efficiency and the luminosity estimate. The estimated cross section 
and uncertainty are 6.0 ± 1.6(stat) ± 1.2(sys) pb and are given in Table 

In order to employ £+jets events in measurements beyond the cross section, a 
more general approach is needed to determine backgrounds. CDF has used a method 
in 162 pb -1 of data which is more similar with the lifetime tagging approach in Run I 
| 97 | 98 | -phe background estimation method was designed to account for the different 
flavor compositions of multijet and H^+jets samples. The fraction of the VF+jets 
background that contains heavy flavor was estimated from Monte Carlo. The mistag 
rate was applied to the QCD background and to the fraction of VF+jets background 
not arising from bb and cc associated production. The measured heavy flavor tag 
rate was applied to the estimated W+ heavy flavor jet fraction. In Run I, this 
approach provided good agreement of background estimation and observed yield in 
W+ one and two jet samples. In combination with the soft-lepton tag selection, the 
fe-tagged Z+jets sample was also observed to provide agreement between observed 
and expectation at all jet multiplicities. 

Signal efficiencies were calculated after correcting for the observed ratio of the 
tag rate for jets tagged with soft electrons in data and a parallel Herwig sample. 
Systematic uncertainties include the residual flavor composition of the data sam- 
ple, as well as potential differences between soft electron tagged b— jets vs. b— jets 
generally. The multijet background was estimated in tagged and untagged samples 
using the sideband method in the $t vs. lepton isolation plane. The heavy flavor 
fraction of the Vy+jets background was estimated from Monte Carlo. This fraction 
was calibrated using an MC-to-data correction factor based on data and MC jet 
samples. The yield of W+ heavy flavor events was estimated by multiplying the 
heavy flavor fraction by the number of events in the untagged sample, and then ap- 
plying the b— tag efficiencies to that. Mistagged M^+jets background was estimated 
by weighting each jet in the untagged sample by the mistag rate. A correction 
was applied to remove the estimated QCD multijet background in this sample. 
Other small backgrounds were estimated from the Monte Carlo. The background 
estimation method for the W^+jets sample was cross-checked by applying similar 
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techniques to a Z+jets sample where there would be little top quark contribution. 
Agreement of data with the expected background contributions was observed. The 
cross section was estimated in the single tag sample to be 5.6t.\'^(stat)~!iQ^(sys) pb- 
A measurement carried out in the double-tag sample was consistent. 

Using a larger data sample of 318 pb -1 , CDF has produced an extensive suite of 
measurements with different tagging treatments. The secondary vertexing scheme 
wa s up dated using a newer tagger that has higher efficiency, particularly at high 
Pt HSU. Event selection remained the same. The signal efficiencies and background 
estimates were also performed similarly to the earlier analysis. The primary ex- 
ception was that the backgrounds were adjusted with an iterative algorithm to 
correct for the top quark contribution in their control samples until the cross 
section measurement varied by less than a percent. The cross section result was 
8.7±0.9(siat)l J"g(sys) pb. A complete measurement was also performed in the same 
data sample with a jet probability taggeJI^. Two different selections (< 1%, < 5%) 
were used to vary the purity and check the consistency of their results. Event se- 
lection was the same as in the secondary vertex analysis with the addition of a 
cut on the transverse mass of the Iv pair M^f > 20 GeV. Signal efficiency and 
background levels were estimated as in the secondary vertex tag analyses. For the 
multijet background, the definitions of the sidebands were adjusted to provide a 50% 
systematic uncertainty in both the untagged and tagged samples. Control samples 
of events with one or two jets validated the background modeling approach. The 
background estimate was corrected iteratively to account for top quark contamina- 
tion until the cross section measurement was stable to within 0.1%. For the most 
sensitive jet probability tag selections, the cross section is 8.9 ± 1.0(stat)~t 1 ' (sys) 
pb. Cross sections measured in the separated electron and muon channels, and in 
the double-tagged subsample, are all consistent with these results. The results in 
318 pb -1 are somewhat higher than the theoretical value. Cross checks in the data 
using the earlier secondary vertex algorithm, as well as looking in the double tagged 
sample provided consistent values of the cross section. In 1.1 fb _1 of £+jets events, 
CDF has produced a preliminary measurement with substantially better statistical 
precision: 8.2 ± 0.5(stat) ± 0.8{sys) ± 0.5 pb ^21. 

D0 has published two determinations of the cross section in collisions at 
y/s = 1.96 TeV using single lepton c hann els with secondary vertex tagging. The 
first utilized 230 pb" 1 of collider datEpS 

and the second incorporat ed 4 25 pb 
a preliminary result based on 900 pb -1 has also been pursued HSH. Trigger 
requirements of one lepton and one jet are the same as for the £+jets analysis with 
topological selection. The data were separated into events with single tags, and 
double tags. 

In the 230 pb^ 1 analysis, signal efficiencies were established by applying a tag 
rate from soft muon tagged dijet events corrected from a Monte Carlo study so that 
it applied to inclusive b— jets. The simulation also provided the c— jet tagging rate, 
and that was corrected for the semileptonic tag efficiency in data vs. Monte Carlo. 
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For the 425 pb _1 measurement, separate multijet samples were denned which had 
different levels of 6 — jet content. One required a non-isolated /x in one jet of a dijet 
pair, the other was the subset of this sample which also had a secondary vertex 
tag of the jet opposite this /i. Event yields were extracted in subsets of these two 
samples with and without cither lifetime or soft-lepton tags. Eight equations were 
formed from these yields which allow the extraction of the 6-tagging efficiency for 
semileptonically decaying 6-quarks. The tagging efficiency for inclusive b— jets was 
obtained from bb Monte Carlo and scaled to a factor from the same eight equation 
separation scheme. 

To estimate the multijet background, the matrix method was employed in tagged 
and untagged events. All other backgrounds, including Z+jets and diboson compo- 
nents, were estimated from the simulation. NLO calculations of cross sections were 
used to normalize the rates of these backgrounds. The VF+jets background was 
found by subtracting the above contributions from the untagged sample and then 
multiplying the remainder by the tag rate appropriate to the flavor mix in VF+jets 
events. This composition is evaluated from the Alpgen W+jets simulation. In the 
425 pb" 1 analysis, an alternative parton matching scheme was used to estimate a 
20% systematic uncertainty on the flavor fractions of the W+jets sample. 

For the 230 pb _1 analysis, Fig. [T2l illustrates the distributions for Ht and a mul- 
tiparameter kinematic discriminant for data, signal and background expectations. 



• D0data 




H T (GeV) Likelihood 



Fig. 12. Distribution of event Ht (left) for the D0 £+jets 230 pb _1 sample with secondary vertex 
tag, and likelihood (right Data are indicated with points and error bars, while the expectation 
for signal plus background are shown via solid histograms. 

Figure [13] shows the jet multiplicity distribution for the single tag and double tagged 
samples in the 425 pb _1 analysis. The one-jet and two-jet bins exhibit good agree- 
ment with the expectation, thus validating the background estimation method. An 
excess is observed in the three and four jet bins which is attributed to top quark 
production. 
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Fig. 13. The jet multiplicity distribution for the D0 425 pb — 1 l +iets sample with secondary 
vertex tag for single tagged (left) and double tagged (right) eventaiSSJ. Data are indicated by 
points and error bars, while signal and background expectations are shown via solid histograms. 



The measured cross section from the earlier data sample is &.&t\'%{stat + sys) ± 
0.6(lum) pb. An alternative approach using an impact parameter tag was also used 
to calculate the cross section. This method produces very compatible results with 
the secondary vertex approach and is actually fairly uncorrelated with it in terms of 
event sample. As with the soft-lepton tag analysis in 425 pb -1 , the cross section was 
extracted by maximizing a likelihood given the observed number of events. Each 
systematic uncertainty was implemented as a Gaussian which could alter the fitted 
cross section. This provided a cross section measurement of 6.6 ± 0.9(stat + sys) ± 
0A(lum) pb. The measurement in 900 pb -1 uses much the same techniques and 
gives 8.3t Q oi{stat)tii{sys) ± 0.5{lum) pb. 

Substantial uncertainties for all D0 and CDF single lepton mode analyses arise 
from tagging efficiency uncertainty, and from the jet energy scale uncertainty. In 
addition, the secondary vertex analyses have significant uncertainty from the back- 
ground modeling. Both CDF results have approximately a five percent lepton iden- 
tification uncertainty. The 425 pb -1 D0 measurement has significant uncertainties 
from pdf's, factorization scale and heavy quark (b and c) mass. 

4.5. All-jets channel 

The all-jets channel is the most copious final state with a branching fraction of 
~46%. Without any energetic neutrinos in the final state, the all-jets mode is also the 
most kinematically constrained, and this allows a full reconstruction of the tt signal. 
However, the signal suffers from an overwhelming background from QCD multijet 
production, with a cross section many orders of magnitude larger, making the ex- 
traction of tt events extremely challenging. To improve the signal-to-background 
ratio, S/B, a set of requirements exploiting the kinematic and topological char- 
acteristics of standard model tt events is applied to the data. Unlike the ^+jets 
topological selection described in Section f4. 4. II kinematic properties alone are not 
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sufficient to isolate the top quark. They are coupled with the use of b— tagging to 
separate the heavy-flavor poor background from the tt signal. The CDF and D0 
experiments have previously measured the tt production cross section (c^) in the 
all hadronic ch annel in Run jl05 | 106 | rp ne s ^ ra ^ e gj es employed by their measure- 
ments in Run II ^^^^Q ^rely heavily on the techniques developed during their Run I 
analyses. Table[4]gives an overview of the various selection requirements applied by 
the CDF and D0 analysis discussed in this section. 



Table 4. Selection cuts for Run II all-jets cross section measurements by CDlC23 
and D0 Ell 



cuts 


CDF 


D0 


preselection 


primary event vertex (Z ver t <) 
lcpton veto 
pjf~ ts > (GeV) 

Njets > 

H° T ets > (GeV) 


60 cm 
yes 

15 

2 

6 
125 


50 cm 
yes 
15 
2.5 
6 

90 


topological selection 


variables 


H't 3 , E^t C,A 
(kinematic selection) 


H T , A, Ey 6 

<V 2 >,M 3 J n ,M 
(neural net selection) 


6- tagging 


yes 


yes 



Discriminating event characteristics 

In addition to the discriminating variables already discussed in section 4.1, the 
all-jets channels also employ some more jet-based variables. The QCD multijet 
production is dominated by a 2 — > 2 parton process producing two hard leading 
jets with extra jets produced through gluon radiation. Therefore, the additional 
jets are expected to be softer in QCD background than in tt signal. The parameters 
used are , the scalar sum of all jet px's except the two leading jets, and E^ 6 , 
the geometric mean of the transverse energies of the 5th and 6th leading jets. The 
properties which are typical for the top quark event structure, owing to the presence 
of IV-bosons and 6-quarks, are also used. The variables used are the mass-likelihood 
M. and the second-smallest dijet mass in the event M^™^. M. is defined as M. — 
(M Wl - M w ) 2 /g 2 Mw + [M W2 - M w fa\ lw + (m tl - m t2 ) 2 /cjf nt where M Wl (M Wa ) 
is the mass of the two jets corresponding to the W boson from the first (second) 
top quark, of mass m tl (m t2 ). The parameters M\y = 79 GeV and <tm w = H 
GeV are the central value and resolution of W boson mass peak, obtained from tt 
all-hadronic Monte Carlo along with the resolution of the top quark mass, cr mt = 
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21 GeV. M. is calculated for each possible assignment of jets to the Ws and b- 
quarks, while only the permutation with the smallest M. is used as the discriminator. 
M. provides good discrimination between signal and background by requiring two 
jet pairs that are consistent with the W boson mass, and two VF+jet pairs that 
are consistent with a single top quark mass of any value. The presence of two 
W bosons in tt events provides significant rejection against the QCD background. 
A further requirement that the two reconstructed top quarks have equal masses 
provides some additional discrimination. The discriminating performance of the 
two chosen kinematic variables, Et and C can be seen in Fig. [14J 
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Fig. 14. Kinematic variable distributions ( Et (left) used by CDF analysiJlEZl and 
(right) used by D0 analysi a^^^ j) in multijet data and tt Monte Carlo simulation. 



Measurement by CDF: 
In Run II, CDfES 

has performed the measurement of the tt cross section (o>j) in 
this channel utilizing approximately 311 pb -1 of data sample selected by dedicated 
multijet triggers. Events with at least four jets with pr >10 GeV and large Ht >125 
GeV are selected at the trigger level with S/B ~ 1/3500. The offline preliminary 
signal selection requires events with a large jet multiplicity, 6 < Nj ets < 8, with 
Pt >15 GeV and \r]\ <2. Application of a veto on isolated high-p^ leptons allow 
to keep the analysis orthogonal to those involving leptonic channels. After these 
requirements 364,006 events are selected for further analysis. 

The tt events are modeled with Pythia and Herwig using m t =178 GeV. In 
order to improve the S/B, the analysis employs a kinematic and topological event 
selection based on Ht, H t \ C, and A optimized to achieve the maximum signal 
significance for tt events, defined as the ratio between the expected signal and the 
statistical uncertainty on the sum of the signal and background. The values for 
the cuts after optimization are: A + 0.005H T j >0.96, C >0.78, and H T >280 GeV. 
Application of the kinematic selection yields 3315 events in data with an efficiency 
of 6.7±1.4% for the tt signal and a S/B ~l/25. In order to further enrich the sample 
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with top-like events, events are required to have at least one b— tagged jet which 
leaves 695 events containing 816 b— tags, thus reaching a S/B of about 1/5. 

The background sources for this final state are due mainly to QCD production of 
heavy quark pairs (bb and cc) and false tags from light-quark jets. The mistag rate 
is evaluated in the exclusive four jet data sample before the kinematic selection and 
is parameterized in terms of jet Et, track multiplicity N tr k, and number of primary 
vertices in the event N vert . To estimate the background, each multijet event in the 
signal sample is weighted by its b— tag probability. The sum of these weights for the 
multijet events gives the expected number of tags from the non-signal processes. 
Before the kinematic selection, when the multijet sample is still predominantly 
composed of background events only, the predicted and observed number of tagged 
jets for different jet multiplicities agree very well as can be seen in Fig. 1151 giving 
confidence in the constructed mistag rate. 

The cross section measurement is performed using the total number of tagged 
jets (not events) in order to avoid the explicit calculation of the background for dou- 
ble tagged events. The average number of tags in a tt event passing the kinematic 
selection is n™ 9 e = 0.846± 0.073, measured in tt simulation by taking into account 
the different tagging efficiencies for b—, c— and light- flavored-jets and correcting for 
the difference in efficiencies in data and Monte Carlo events. After application of 
the kinematic selection, a total of N i, s —816 candidate tags are observed, whereas 
the expected number of tags from background sources in the signal region amounts 
to iVb/cg =684±38 tags, after correcting for the presence of tt events in the pretag 
sample. The resulting excess in data of tagged jets in the signal region is ascribed to 
tt production. The tt production cross section is determined by = ^"^n^'xc 
and is found to be a t i — 7.5 ± 2A(stat.)^\(sys)^ \{lurri) pb for m top =178 GeV. 
The dominant systematic uncertainty of ~20% arises from the dependence of selec- 
tion efficiency on the jet energy scale. In Fig. 1151 the distribution of the number of 
observed tags and background is compared to the tt signal expectation assuming 
the measured cross-section 7.5 pb. 



Measurement by D0 : 

In Run II, D0| 1O8 | has measured the tt cross section (c^) in this channel based on 
~405 pb -1 of data collected using specific multije t tri gger. The analysis proceeds 
on a methodology similar to that used in Run ff^S. In addition to the single- 
tagged events, the measurement has been extended to the double-tagged events. 
The analysis proceeds with the preselected data sample composed of events with 
>6 reconstructed jets with p T >15 GeV, | ^7 1 <2.5 and H T >90 GeV. The bulk of 
the background is rejected by requiring the presence of secondary vertex tagged 
jets In order to maximize acceptance and sensitivity, separate samples with one 
tagged jet, or two tagged jets were analyzed. In the former, a tight requirement 
{Lxy/vL^y > 7) was placed on the tag to control backgrounds. In the latter case, a 
looser requirement (L xy /ar l:i:y > 5) was used to maximize efficiency. 
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Fig. 15. Left: Average number of tags per event observed in the multijet sample before kine matic 
selection compared with the estimate from the tag rate parameterization in the CDF analysiJSE] 
Right: Number of tags observed in multij et d ata after kinematic selection compared with the 
expected background in the CDF analysi rP I The tt expectation refers to the measured cross- 
section of 7.5 pb. 



Further suppression of the multijet background is achieved by applying a neural 
network (NN) selection based on six discriminating kinematic variables: Ht, -4, 
E T ' 6 , <rf >, M%£ n , M, taken from the equivalent Run I analysis^!]. 

The NN is trained to force its output near 1 for tt events and near -1 for 
QCD multij et ev ents, using the multilayer perceptron network in the Root anal- 
ysis prograrrP^. The very small S/B in the untagged data sample provided the 
background input for the training of NN, while tt Monte Carlo simulation using 
Alpgen and Pythia with m t =175 GeV was used for the signal. Fig. [TBI shows the 
performance of the NN in discriminating tt signal from the multijet background for 
both single- and double-tag samples. The NN distributions vary for the two sam- 
ples due to the variation of their flavor content. Overall, the NN displays significant 
discriminating power for the single- and double-tag events, although it exhibits a 
better discrimination for single-tag samples. 

The top quark cross section is extracted from the output of the NN. The back- 
ground dominated preselected sample is used to estimate the background. For the 
loose and tight tags, the mistag rate is derived from the data with N tag s _!■ It is 
parameterized in terms of the pr, rapidity y, azimuthal angle ef> of the jet and the 
location of primary vertex along the beam direction zpi/, in four different Ht bins. 
Fig. [17] show the distributions of the NN output for the data, the Monte Carlo 
simulation prediction for a t i =6.5 pb, the predicted background, and the predicted 
signal+background. It can be seen that background keeps dominating over the sig- 
nal even at large values of NN. The output of the NN is used to select the sample 
enriched in tt signal by applying the cut NN >0.81 (0.78) for the single (double)-tag 
analysis, optimized to minimize the fractional statistical error on the cross section 
measurement. 
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Shown are the data (points), background (the line histogram), signal (filled histogram) and sig- 
nal+background (dashed histogram). The vertical line represents the used cut of NN>0.81 (0.78). 



After the NN requirement, 495 (439) events are observed while the predicted 
number of background events is 464±4.6 (400 jig'^) in the single (double)-tag analy- 
ses. The kinematic selection efficiency times branching ratio for the tt signal in the 
two analyses are = 0.0242± u ;^|(0.0254± u ;^ u ), respectively. Based on the clear 
excess of observed events over the predicted background, D0 has measured a cross- 
section of a t i — 4A^Q(stat.)t.Q'g(sys)±0.3(lum) pb for the single-tag analysis and 
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a t i — A.Ttz^istat^ti'Xisys) ±0.3(lum) pb for the double-tag analysis. A combined 
cross-section measurement of a t i = &.5t_i'%(stat.)t\[i(sys)±0.3(lum) pb is obtained 
for a m t =175 GeV. The major systematic uncertainty in the analysis arises from 
the dependence of selection efficiency on the jet calibration and identification. 



4.6. Combined cross section 

The most sensitive of each type of cross section measurement discussed in this sec- 
tion are given in TableEl For collisions with y/s = 1.8 TeV, the quoted measurements 
are final results. Work is still progressing on the s/s = 1.96 TeV measurements. The 
input m t is 175 (178) GeV for many of the D0 (CDF) measurements. Generally, the 
measured cross section will be higher if m t is lower because the efficiency to select 
top decreases somewhat with mt- For the dilepton and i'+jets channels, the slope 
of measured cross section with difference in mass from the input value (Act/ Am) 
is ~ 0.1. The all-jcts channels exhibit significantly more variation. The range over 
which the Act is quoted at least includes 170 — 180 GeV, and several measurements 
test from 160 to 190 GeV. Currently, the total uncertainty is reaching the 15% level. 
This should still be reduced as higher statistics are accumulated. More events will 
also translate into somewhat lower systematic uncertainties, largely because of the 
increase in samples dedicated to parametrizing detector performance such as the 
jet energy scale. 



5. Single Top Quark Searches 

Twelve years after the discovery of the top quark via strong pair production, the 

first evidence of e lectr oweak production of single top quarks has been reported by the 

D0 collaboration^. This search is much more challenging than tt production due 

to smaller cross sections and larger backgrounds. A general overview of single top 

quark production is given in Sec l2.2l Here we discuss the signatures and backgrounds 

associated with the analyses and review the results by both the CDF and D0 

experiments. We review these results with a historical perspective and begin with 

early single top quark search analyses which resulted in an upper bound on the 

production cross section (section [5]). Later, in section [572]. we review the analyses 

which JciL, to the cyidcricc for-the single top quark production. ... . 
Signal Kinematics ligure 1181 shows some oi the leatures of the single top quark 

events produced in the s- and f-channels. The plots of transverse momentum and the 

pseudorapidities of objects in the single top quark events shown in this figure reveal 

that for the t-channel process the 6-quark has very low transverse momentum and is 

peaked at large pseudorapidity (forward region). Thus the 6-quark in the ^-channel 

is difficult to detect. 
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Table 5. Measured cross sections for Tevatron experiments, including the function for different top 
quark masses. The channel, experiment, collider energy and integrated luminosity are listed in the first 
four columns. The measured cross section is in column 5. Preliminary results are indicated by * in 
column 1. 



Channel 



Expt. y/s (TeV) Lum. (pb" 1 ) 



(Pb) 



all-jets //i-tag EE] 

... II OKI 



all-i 



"to. . 

l-jets/fc-tagEH 
U,l+)ets, all-jets EH 
11, l+jets, all-jetsEH 
U, £+track W 
«,£+trackE3 
gg* 1821 

^+track* E2 
«,£+track* 78 
£+jcts/topo ^ 
£+jets/topoEHl 
£+jcts/topo*IB2l 
( I jt'ts//i— tag^-i 
^+jets/^-tag*IS31 
£+jcts/^-tag*Ei] 

£+jcts/6-tag-kin ^ 

£+jets/jprobli22| 

f+jets/b-tagEffl 

^+jets/6-tag 1221 

£+jets/fe-tag*IS2l 

l+jets/6-tag*DE] 

all-jets/fe-tagLUiTJ 
all-jets/fe-tagEH 



D0 
CDF 

D0 
CDF 



1.8 
1.8 
1.8 
1.8 



110 
110 
125 
110 



7.1 ± 2.8(stat) ± l.5(sys) 
10.1 ± 1.9(stat)±l\(sys) 
5.69 ± l.2l(stat) ± 1.04(sys) 
6-5t^4(stat+sys) 



CDF 
D0 
CDF 
CDF 
D0 



1.96 
1.96 
1.96 
1.96 
1.96 



197 
425 
1200 
1000 
1000 



7.ot 2 2 \t( stat )-i.i( s y s ) ± °- 6 ( lum ) 

7.4 ± lA(stat) ± l.O(sys) 

6.2 ± l.l(stat) ± 0.7(sys) ± 0A(lum) 

8.3 ± l.3(stat) ± 0.7(sys) ± 0.5(lum) 
6.2±l^(stat)±^(sys) ± 0A(lum) 



CDF 
D0 
D0 



1.96 
1.96 
1.96 



191 
425 
900 



6.6 ± l.l(stat) ± l.5(sys) 



6.1 



+ 1.3 

-1.2 
,+0.9 



(stat) ± 0.7(sys) ± 0A(lum) 



6.3Zo S (stat) ± 0.7(sys) ± OA(lum) 



CDF 
CDF 

D0 



1.96 
1.96 
1.96 



191 
760 
425 



5.3 ± 3.3(stat)f^J{sys) 
£ 1.7(siat)^;J(sys) 
(stat + sys) ± 0.4(lum) 



7.3 



+2.0 
1.8 



CDF 1.96 162 6.0 ±1.6( stat. ) ±1.2(sys) 

CDF 1.96 318 8.9 ± l.Oistatf+Hisys) 

CDF 1.96 318 8.7±0.9(stat)t 1 1 9 (sys) 

D0 1.96 426 6.6 ± 0.9(stat + sys) ± 0.4(lum) 

CDF 1.96 1100 8.2±0.5(stat)±0.8(sys)±0.5(lum) 

D0 1.96 900 8.3to[t(stat)+ 1 ^ ) (sys)±0.5(lum) 



CDF 
D0 



1.96 
1.96 



311 
405 



7.5 ± 2.l(stat)±ll(sys)^l(lum) 
4.5±_fl(stat)±\-j(sys)±0.3(lum) 



Background processes Associated production of a W boson with jets (W+jets) 
and pair production of top quarks (tt) are the two most dominant processes which 
mimic the single top quark signatures. Events with heavy flavor produced in asso- 
ciation with a W boson give rise to irreducible backgrounds. These are Wbb and 
Wbbg for the s- and t-channel respectively. In general, the tt events have larger jet 
activity, but they can result in signatures similar to the single top quark in case the 
jets are merged, mis-reconstructed or are outside the fiducial volume. Leptons can 
also be lost in the detector cracks or lie outside the fiducial volume considered, tt 
decays can look very similar to single top quark events when one W boson decays 
to ev or \xv and the other decays to tv. Additional multijet background comes from 
events containing a mis-identified isolated lepton (e or /i) associated with hadronic 
jets. Z/7*+jets contribute to the backgrounds if one of the two leptons from Z 
decaying to e + e~ or is lost. Diboson events (WW, WZ) also contribute to 

the background, though at a much smaller level. 
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Fig. 18. Distributions of transverse momenta (a) and pscudorapidity (b) for the final state partons 
in s-channel single top quark events (top row) and for t-channel single top quark events (bottom 
row). The histograms only include the final state of t, not t 



Modeling single top quark and backgrounds The two main sources of back- 
ground are examined: ti production and non-top quark backgrounds. The non-top 
quark backgrounds include VF+jets and mis-identified multijet events. W+jets (in- 
cluding heavy flavor process qq' — > Wg with g — > bb or g — > cc, and gq — > Wc) are 
included in their relative proportions estimated using AlpgeiJ^ and normalized 
to the data. 

Both the CDF and D0 collaborations have performed searches for sing le top 
quark production in the data samples gathered in Run I (1992-1996 fi 12 | 11 fo Data 
collected during the first few years of Runll have led to updated analyses by CDF 
and D0. All of these measurements have led to upper limits on the production cross 
section of the single top quark as listed in Table [S] These are reviewed below. 



Table 6. 


95% C.L. upper 


limit 


on single 


top quark production cross 


section (in 


pb). 






Theory 


Run I 
CDF D0 




CDF 


Run II 


D0 




Luminosity (pb ) 

s-channel 

t-channel 


0.88 (±8%) 
1.98 (±11%) 


< 18 

< 13 


160 

< 17 

< 22 


Q xi) 
695 

< 14 

< 10 


NN 
955 
< 3.2 
<3.1 


LD 
230 

i +0 - 7 

2 +0 - 9 


Cut 
230 

< 10.6 

< 11.3 


NN 
230 

< 6.4 

< 5.0 


LD 

< 5.0 

< 4.4 
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5.1. Search strategy 

In this section we review the searches for the single top quark by the CDF and D0 
experiments in Run II. The CDF experiment published a search for the single top 
quark using a data sample corresponding to an integrated luminosity of 160 pb _1 UJ4l 
and had provided updated results with 600-955 pb~ 2 of data in the year 2006. The 
D0 experim ent performed a search for the single top quark using 230 pb _1 of 

The strategy of the analyses pursued during Run II is as follows: First a set 
of events with the signal topology is selected. Next, they are separated into inde- 
pendent sets based on the flavor of the lepton (e or /x) and the number of 6-quark 
jets (exactly one 6-jet or two or more 6-jets). The different sets are analyzed inde- 
pendently by statistical techniques to enhance the signal. They are then combined 
in the final stage. Finally, in the absence of any significant observed signal, binned 
likelihood fits are performed on the outputs of the statistical analysis methods to 
obtain limits on the production cross section for single top quarks. 

The single top quark search starts with the selection of events with at least one 
high pt lepton (electron or muon), significant $t and between two to four jets. At 
least one of the jets is b— tagged to separate signal events with heavy flavor from 
the VF+light jet background event. Separation of samples into s-channel events and 
f-channel searches is accomplished by requiring at least one untagged jet in the 
f-channel analysis. For both the channels, selected samples are further subdivided 
into four orthogonal sets based on the flavor of the lepton (e or /i) and the number 
of b quark jets: "single tagged" (with exactly one 6-tagged jet) or "double tagged" 
(two or more 6-tagged jets) samples. At this initial event selection the background 
is still expected to be about 90% of the sample. 

5.1.1. Collection of input variables 

To further reduce the backgrounds, both kinematic (object and global event) and 
topological information is used. Some of the many variables which are used for dis- 
criminating against the backgrounds are: jet px for different jets; b— tag information 
of the jet; H (total energy); Ht (total transverse energy); M (invariant mass); Mt 
(transverse mass); summing over various objects in the event; jet-jet separation; jet 
pseudorapidity (i-channel); top quark spin; etc. They are selected base d on exten- 
sive analysis of the Feynman diagrams of signals and background s^ The list of 
variables used by D0 analyses is shown in Table[7](CDF analyses use a subset). The 
best jet is defined as the jet which together with the W boson leads to an invari- 
ant mass closest to the top quark mass of 175 GeV. In the s— channel (t— channel) 
analysis, the top quark is reconstructed from the W boson and the 'best jet' (the 
leading b— tagged jet). The W boson is reconstructed from the isolated lepton and 
the The z-component of the neutrino momentum is calculated by constrain- 
ing the lepton and neutrino to the W boson mass, and the solution with smaller 
z-component of the neutrino momentum is chosen from the two possible solutions. 
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Table 7. List of discriminating variables. A check mark in the final four columns indicates in which signal-background 

Mil I 

pair of the neural net analysis the variable is used I I 



Variable 



Description 



Signal-Background Pairs 
s-channel t-channel 
Wbb it Wbb ft 



PT (jetl taggod ) 

Pr(jetl untaggod ) 

p T (jet2 untagged ) 

PT(jetl non _ bcst ) 

PT(jet2 non _ bcst ) 

Pt (jetl,jet2) 
M T (jetl,jet2) 
M(alljets) 
_H" T (alljets) 

(alljets — jetl t 



aggcd 



M(alljets- jetl taggcd ) 
ff(alljets- jetl taggod ) 

H T (alljets-jetl taggcd ) 

M(W,jetl taggcd ) 

M(alljets - jet best ) 
H(alljets- jet bcst ) 
ff T (alljets-jet bcst ) 

M(W,}et hcst ) 



^(jet! untagged) x Qe 
ATC(jetl,jet2) 

cos(£, jetl untagged )top t 



cos(i, <2f xz)to Pbcat 

cos(alljets,jetl taggcd ) alljcts 

cos(alljets,jet non _ bcst ) alljcts 



V 
V 

V 
V 
V 

V 



Individual object kinematics 

Transverse momentum of the leading tagged jet y/ y/ 

Transverse momentum of the leading untagged jet 
Transverse momentum of the second untagged jet 
Transverse momentum of the leading non-best jet 
Transverse momentum of the second non-best jet 

Global event kinematics 
Invariant mass of all final state objects 
Transverse momentum of the two leading jets 
Transverse mass of the two leading jets 
Invariant mass of all jets 
Sum of the transverse energies of all jets 
Transverse momentum of all jets excluding the leading 
tagged jet 

Invariant mass of all jets excluding the leading tagged jet 
Sum of the energies of all jets excluding the leading tagged 
jet 

Sum of the transverse energies of all jets excluding the lead- 
ing tagged jet 

Invariant mass of the reconstructed top quark using the 
leading tagged jet 

Invariant mass of all jets excluding the best jet 
Sum of the energies of all jets excluding the best jet 
Sum of the transverse energies of all jets excluding the best 
jet 

Invariant mass of the reconstructed top quark using the 
best jet 

Angular variables 

Pseudorapidity of the leading untagged jet X lepton charge — — 
Angular separation between the leading two jets *J — 

Top quark spin correlation in the optimal basis for the t- — — 
channeJUZl, reconstructing the top quark with the leading 
tagged jet 

Top quark spin correlation in the optimal basis for the s- y/ — 
channefn^, reconstructing the top quark with the best jet 
Cosine of the angle between the leading tagged jet and the — 
alljets system in the alljets rest frame 

Cosine of the angle between the leading non-best jet and *J 
the alljets system in the alljets rest frame 



V 



V 
V 

V 
V 



V 



V 
V 



V 



V V V 

- V 

- v - 

V 

V 



V 



V 

V 

V 

V 
V 

V 

V 



v 7 



5.1.2. Analysis methods 

Several different analysis techniques are employed to separate the signal from the 
background: Cut based, Neural Networks (NN), Likelihood and Matrix Element 
based analyses. Generally the analyses are optimized separately for s-channel and 
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f-channel analyses. For example, the D0 Neural Network analysis focuses on re- 
jecting dominant backgrounds (VF+jets and tt) by training the neural networks or 
optimizing the cuts separately for each background and for each lepton type. This 
leads to eight separate sets of cuts or networks (two leptons (e, fi) x number of 
tagged b jets (= 1, > 2) times signal type (s, i-channel)). 

• Cut Based analysis: This is the traditional method of applying simple 
selection requirements on a set of variables. In the D0 analysis, each of 
the input variables is assigned a performance metric which is evaluated by 
computing the best expected limit for a given cut on each of the variables. 
The variables with the best performance are then combined by a simple 
AND in order of an assigned rating based on their relative sensitivity. An 
optimal cut point is also evaluated for each intermediate set of combined 
variables. Again, the set with the best expected limit is chosen for the final 
analysis. This optimization is performed for each of the channels separately. 
For each channel, the optimal sets of variable and cuts are listed in Ref. 11111 

• Template Fit Analysis: Another technique used by CDF involves a max- 
imum likelihood fit to the distribution of Q ■ t], where Q is the lepton charge 
and r\ the pseudorapidity eta of the light quark jet of the events, in order to 
extract the s-and i-channel signal from the data. The variable Q-rj is chosen 
as it shows the largest difference in kinematics between s and t— channel 
events (see two left plots of Fig. fT9]) . Templates of Q ■ r\ distributions from 
expected s and t— channel signals, and the two expected background sources 
are used for the likelihood fit. To extract inclusive s and t— channel content 
a template fit is carried out using the Ht distribution, as this distribution 
is similar for both production channels. This is shown in the right plot of 
Fig. HU The likelihood fit takes into account the systematic uncertainties 
from the jet energy scale, the top quark mass, 6-tagging efficiencies and the 
luminosity. Smaller systematic uncertainties from initial state radiation, fi- 
nal state radiation, parton distribution functions, the choice of signal Monte 
Carlo generator, trigger, identification are also included. Some systematic 
effects can change the shape of the Q ■ r\ distribution, which is taken into 
consideration as well. 

The fit to the Q ■ rj distributions using a Bayesian method, leads to a 
95% C.L. upper limit on the production cross sections both in the s and 
t— channels independently. The fits using the Ht distribution result in the 
combined s and t— channel upper limits. These limits are listed in Table [6] 



package is used by D to analyze (testing and training) the events. The 
"NeuroBayes" package^A^ is used by CDF for their analysis. The networks 
were used with three layers of nodes: input, hidden, and output. Monte 
Carlo simulated events are used for the training and testing of the net- 
works. Optimization studies based on the expected upper limits on the 



• Neural Network Analysis: The Neural Network 
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single top quark production cross sections indicate that each of the chan- 
nels require different networks for treatment of the dominant backgrounds. 
D0 separates backgrounds into two categories, Wbb and it events. CDF 
has an additional category for Wcc. This leads to multiple networks (be- 
tween three to eight depending on the experiment) for the two s-channel 
or i-channel analyses. As an example, the set of input variables used for 
the D0 analysis are listed in Table [7] and denoted by a check mark. For 
each network, the set of input variables could be different (see Table [7|) and 
the selected combination of variables is chosen to give the largest signal to 
background separation while keeping the testing error at its minimum. 
The outputs of the NN for the D0 data and the expected signal and back- 
grounds are shown in Fig. [20] One notes that the NN output peaks near 
one for signal events and at low values (near zero) for background events. 
A good separation between the signal and it backgrounds is seen for the it 
networks. The Wbb networks are not as efficient in distinguishing between 
signal and backgrounds since the VF+jets topology is closer to the signal 
events than the it events. Similar behavior and separation is found in the 
CDF analysis as well. 
• Likelihood Discriminant Another popular technique for extracting small 
signals from large datasets is by constructing a likelihood discriminant C(x) 
from a vector of measurements x: 

£fg) — Psignal (x) ^ 

Psignal (■£) ~t~ Pbackground (^) 

where, P S tgnai{x) and Pbackground(x) are the probability density for the sig- 
nal and background events respectively. The strength of this method relies 
on the use of the difference in the shapes of the distributions for the signal 
and background events for a given variable. Both CDF and D0 analyses 
separate the signal into s-and t-channels while constructing the likelihoods. 
Monte Carlo events are used to construct the one dimensional probability 
density functions P S i gn ai{x) and Pbackground(x) for each of the input vari- 
ables. The final probability density functions of the signal and backgrounds 
are the products of the functions constructed for the individual variables. 
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Fig. 20. Comparison of D0 background, signal, and data for the neural network outputs, for 
the electron and muon channels combined, requiring at least one tag, for the tb-tt filter (a), the 
tqb-tt filter (b), the tb-Wbb filter (c), and the tqb-Wbb filter (d). Signals are multiplied by fO for 
readability ESI 

No correlations between variables are used. One expects the value of the 
likelihood discriminant C(x) to be peaked near one for signal events and 
near zero for background events. Many different input variables and combi- 
nations were considered during the optimization process to obtain the best 
expected limit on the cross section measurement. The list of input variables 
used for the D0 and CDF analysis are similar to those used by NN. The 
Likelihood Discriminant for the observed single tagged data (electron and 
muon data combined) together with the expected signal and backgrounds 
are shown in Fig. [21] 

5.1.3. 95% C.L. Upper Limit on production cross section 

The cross section for single top quark production is computed from the 
observed data using a Bayesian approach. The probability to observe the 
vector of event yields d assuming the single top quark production cross sec- 
tion a is given by the single top quark acceptance a, integrated luminosity 
C, and the number of events expected from background b, is given by 



p(d\cr,a,b) = P{d\a(T + b) 
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D0 Run II Preliminary, 370 pb 1 

— t-channel (x10) 
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W+jets, WW, WZ J_ 

■ Multijet _ 
• Data 
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+ 
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i d tab / 

Discriminants for single 6-tagged events from the D0 experiment -^ -U 

where a — a£ and P(x\y) is the Poisson probability to observe x event 
yield when y are expected. Using Bayes' theorem the probability for the 
single top quark production cross section to have the value a is then 

p(a\d) = -jj. J J p(d\a,a,b)TT(a,b)dadb, 

where J\f is a normalization constant, a flat prior for a was assumed and 
the prior for the other parameters n(a, b) is a product of Gaussians with 
widths given by the experimental uncertainties in these parameters. The 
dependence on all parameters except a was eliminated by integrating over 
these so-called nuisance parameters. 

The measured cross section is then given by the value of a for which p(<j|d) 
is maximized. If this occurs for a = an upper limit a ma x can be set at 
confidence level /? using the condition 



I*™* p(<j\d)da = [3. 
Jo 

he errors (5+ and 5— a: 
/ %(cr|d)dcr = 0.6827 



For a finite cross section the errors <5+ and 5— are defined by 

„+<54 



for which J_ + S + is minimized. 

Since the observed data are consistent with the background predictions for 
all analysis techniques used by D0 and CDF, following the prescription 
for cross section computation described above, the 95% C.L. upper lim- 
its on the single top quark production cross sections are computed. These 
measurements are listed in Table |6] These upper li mits represented signifi- 
cant improvements over previously published resultJU^l, mainly due to the 
larger data sets analyzed, and the use of multivariate analysis techniques. 
It is interesting to note that they approach the region of sensitivity for 
models incorporating fourth quark generation scenario or flavor-changing 
neutral-currents . 
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5.2. Evidence for single Top quark production 

The evidence for single t op q uark production was reported by the D0 Col- 
laboration in December 2008 ^^ 1 This analysis is based on a large data sample 
corresponding to a luminosity of 0.9 fb _1 collected between 2002 and 2005. Event 
sample selection procedures are similar to the earlier searches as described in the 
previous section. Events are classified into twelve subsamples based on the lepton 
type (e or n), the number of jets (two, three or four jets), and the number of tagged 
6-jets (one or two tags) in the event. The event yields are tabulated in Table [8] 



Table 8. Numbers of expected and observed events in 
0.9 fb" 1 for e and /i, one b tag and two b tag channels 
combined from the D0 experiment. The total background 
uncertainties are smaller than the component uncertainties 
added in quadrature because of anticorrelation between the 
IV+jets and multijet backgrounds resulting from the back- 
ground normalization procedure 1112] 



Source 


2 jets 


3 jets 


4 jets 


s-channel 


16 


3 


8 


2 


2 


1 


t-channel 


20 


4 


12 


3 


4 


1 




39 


9 


32 


7 


11 


3 


tt^£+jcts 


20 


5 


103 


25 


143 


33 


Wbb 


261 


55 


120 


24 


35 


7 


Wcc 


151 


31 


85 


17 


23 


5 


Wjj 


119 


25 


43 


9 


12 


2 


Multijets 


95 


19 


77 


15 


29 


6 


Total background 


686 


41 


460 


39 


253 


38 


Data 


697 




455 




246 





A multivariate technique known as "Decision Tree"' I is used for discriminating 
between signal and background events in the sample. Decision Trees, DT, are a 
machine learning technique which are used to compute the probability of the event 
to be a signal event. Decision trees provide a way to represent rules underlying 
data with hierarchical, sequential structures that recursively partition the data. It 
is a binary tree with a selection cut implemented at each node such that each event 
eventually ends up into a well-separated class called "leaves" . A purity value defined 
as the ratio of signal to background events and predetermined from the training 
samples is associated with each of the leaves. The output of the decision tree is 
the distribution of the final purity values for the sample and is discrete. A major 
improvement in this techniq ue, c omes from the implementation of the adaptive 
boosting algorithm AdaBoostP^^ll. This boosting algorithm modifies the weights of 
the misclassified events and rebuilds the tree iteratively. The final decision tree with 
improved performance is an average over many trees produced during the boosting 
process. 
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Forty nine variables (see reference f or details) are used as inputs to the 
Decision Tree. These variables based on individual object kinematics, global event 
kinematics, and angular correlations are constructed in order to discriminate be- 
tween signal and background events. For each of the three searches s— channel, 
t— channel and s + t channel, decision trees were trained on the twelve event sub- 
samples, leading to a total of thirty six decision trees. Thus each one of the 36 
decision trees is trained for a given signal signature against the sum of VF+jets and 
ti backgrounds. The Monte Carlo samples used for training the trees and testing 
their performance are independent of each other. 

Tests on the performance of the decision trees on the collider data are done by 
studying samples which are rich in background. Two such samples were created: 
a W+jets sample and a ti sample. The W+jets sample was selected by requiring 
Ht < 175 GeV for the events with two jets and one tagged 6-jet. The ti sample 
comprises of events with four jets, one tagged 6-jet and Ht > 300 GeV. In Fig. [221(a) 
and (b), the decision tree discriminants, Out, are shown for both the background 
enriched data samples. The prediction of the background models are overlayed on 
these data distributions. Good agreement is visible between the prediction and the 
observed events in these background enriched samples. 

Output of the decision tree analysis in the high end of the discriminant spectrum 
for combined e and \x data samples is displayed in Fig. [22] (c). The expected signal 
and backgrounds are overlaid and summed over the twelve decision trees for s + 1 
channels. A single top quark signal is well accommodated with the data sample. The 
characteristics of the events in the high decision tree discriminant region have been 
examined, by making a selection Odt > 0.65. For these selected events, in Fig. [221(d) 
the distribution of invariant mass from the best b— jet and reconstructed W boson is 
displayed. Good agreement is observed between the data and a prediction including 
background and a single top quark signal. 

Two other supporting analyses were carried out. One is based on using Bayesian 
Neural Networks and t he ot her is known as the Matrix Element technique. The 
Bayesian Neural networlsP^Sl at the first stage is similar to a simple Neural Network 
as described in section [531 Instead of choosing one set of weights to characterize the 
network with Bayesian Neural Networks, a posterior probability density is computed 
using all possible weights. The final network is obtained by computing a weighted 
average of a large number of networks, obtained after many sets of training cycles, 
where the weights are the probability of each network given the training data. The 
BNN uses a subset of the input variables used in the DT analysis. The training 
samples are the combined s and t channel Monte Carlo signals and the combined 
background set. The performance of this technique is shown in Fig. [23] (left panel), 
for events with electron and two jets, with one of them tagged as a 6-jet. 

The Matrix Element method encodes all kinematic information of the event and 
contains all the properties of the interaction. Hence this technique uses maximal 
information and is firmly anchored with an understanding of the underlying physics 
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Fig. 22. Boosted decision tree output distributions for (a) a W+jets-dominated control sample, 
(b) a if-dominated control sample, and (c) the high-discriminant region of the sum of all 12 inclusive 
tb+tqb decision trees. For (a) and (b), Ht = E^, + Et + Y)E T i eiB . Plot (d) shows the invariant 
mass of the reconstructed W boson and highest-p^ b-tagged jet for events with Odt > 0.65. The 
hashed bands show the ±1 standard deviation uncertainty on the background. The expected signal 
is shown using the measured cross section (D0 experiment) HE] 




tb+tqb ME Discriminant 



Fig. 23. Bayesian Neural Network discriminant (left) and the Matrix Element discriminant (right), 
in the high-discriminant region, for the combined event sample (D0 experiment) 11221 

processes. Matrix elements of the main signal and background diagrams are used 
to compute an event probability density for signal and background hypotheses. The 
signal and background Feynman diagrams used for the subsamples with two jets 
and three jets are shown in Fig. [2H A discriminant Ome is calculated from the 
probabilities for the event to be compatible with the signal hypothesis, P S i gn ai and 



Review of Top Quark Physics Results 53 



background hypotheses Pbkg- 

rt Psignal 

Ome = p (5) 

■T signal+P bkg 

where P S i gn ai is the properly normalized differential cross section for an event 
which contains objects with reconstructed four-momenta x: 

1 da 
a dx 

This technique is similar to those employed in the top quark mass measurement and 
shares some of the same tools. Figure [23] (right panel) shows the high end region 
for Ome obtained from the combined event sample and it accommodates the single 
top quark signal reasonably well. 



Psignal — ~ x ~7^Z (6) 




Fig. 24. Representative Feynman diagrams corresponding to the leading-order matrix elements 
used for event probability calculation. The upper and middle row are for events with exactly two 
jets. Upper row, signals: ud— >tb, ub^td; middle row, backgrounds: ud^Wbb, sg^Wcg, ud-^W gg. 
The bottom row shows diagrams for events with exactly three jets. Left two plots: signals, ud-^tbg, 
ug^tbd; right plot: background, uc(— > Wbbg 11221 



A detailed study of the systematic uncertainty was carried out. Dominant sources 
of systematic uncertainties which are accounted for in the analysis are: normaliza- 
tion of the three major sources of the backgrounds (it, W + jets and multijet 
backgrounds) which includes a component arising from the heavy flavor fraction, 
the 6-tagging rate functions for the signal and backgrounds, the jet energy scale 
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uncertainty, and the uncertainty on the integrated luminosity. Some of these uncer- 
tainties are dependent on the shape of the underlying spectrum. Uncertainties were 
assigned for each of the backgrounds, as a function of jet multiplicities, lepton type, 
and number of tags. To derive the affect of each of the uncertainties, the inputs 
were shifted by ±1<t, and the analysis was redone. Systematic uncertainties on the 
signal acceptance were also computed in a similar fashion. 

Given that an excess of events compared to the background estimates is ob- 
served, the cross section is computed using the method described in Sec. 15. II Before 
the final cross section results were computed, a verification of the cross section com- 
putation procedure was performed by generating many sets of pseudo-experiments 
or 'ensembles'. These ensembles were subjected to the full analysis chain, including 
systematic uncertainties. 

Generated ensembles include: 

(1) Ensembles at a few different values of total s+t channel cross sections. 

(2) SM ensemble with inclusive a(s+t channel) = 2.9 pb. 

(3) Ensembles at the experimentally measured inclusive s+t channel cross section. 

(4) zero-signal ensemble 

A pool of 1.7 M weighted signal + background events are used to generate 
the zero signal ensemble. In the ensemble generation process, the relative and total 
yields of each type of background is fluctuated in proportion to systematic errors. To 
generate one of the ensembles, a random sample is drawn from a Poisson distribution 
about the total yield. 

All generated ensembles are used to evaluate the linearity of the cross section 
measurement by comparing the observed cross sections with the input cross sec- 
tions. The three techniques were tested and display a good agreement with a linear 
response function. 

The SM ensemble is used to estimate the compatibility of the measured value 
with the SM expectations. The ensemble generated at the measured cross sections 
were used to verify the estimates of the systematic and statistical uncertainties. The 
zero-signal ensemble is used to determine the sensitivity of each measurement. The 
expected p— value is defined as the fraction of zero-signal ensembles in which a SM 
cross section of at least 2.9 pb is measured. The observed p— value is computed as a 
fraction of zero-signal pseudo-datasets in which at least the observed cross section 
is measured. 

Final results of the measured cross sections and the significance are given in 
Table [9] for the primary DT analysis, together with the two supporting analyses: 
BNN and ME. 

Figure [25] shows the first measurements of the single top quark cross sections at 
the Tevatron. The plots are from the DT analyses. 

The measured cross section for single top quark production a(pp~^s + t — 
channel) = 4.9 ± 1.4 pb, which is consistent with expectations from SM. The 
observed p— value, or the probability that the background fluctuates to give the 
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Table 9. Results from the different analysis techniques for measure- 
ment of the cross sections and significance of the analyses from the D0 
experiment. SD below denotes the number of standard deviations. 11121 





DT 


ME 


BNN 


o~(pp—*s + t — channel) 


4.9 ± 1.4 pb 


4-6±i;l Pb 


5.0 ± 1.9 


expected p-value 


1.9% 


3.7% 


9.7% 


observed p-value 


0.035% 


0.21% 


0.81% 


observed significance 


3.4 SD 


2.9 SD 


2.4 SD 




tb+tqb Cross Section [pb] 



Fig. 25. Expected SM and measured Bayesian posterior probability densities for the tb+tqb cross 
section. The shaded regions indicate one standard deviation above and below the peak positions 
(reproduced from Ref. IllOl l. 

measured cross section value of 4.9 pb or greater is 0.035%, which corresponds to a 
3.4cr evidence for the single top quark production at the Tevatron. 

The CDF collaboration confirmed the D0 evidence for single top quark produc- 
tion in August 2007 by analyzing a data sample corresponding to an integrated lumi- 
nosity of 1.5 fb _1 . CDF performed analyses using two different techniques: a mul- 
tivariate likelihood technique and the matrix element discriminant techniqu ci^^H 
These are extensions of the techniques already described earlier and applied to 
smaller datasets. The single top production cross section measured by the matrix 
element technique is a(pp^s + t — channel) = 3.0^4 pb, which is consistent with 
expectations from SM and those measured by the D0 experiment. This result corre- 
sponds to a 3.1(7 excess over SM background. The likelihood discriminant exhibits 
a 2.7cr excess over SM backgrounds and leads to a compatible single top quark 
production cross section measurement. 
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6. Measurements of the Top Quark Mass 
6.1. Role in electroweak model 

Mass of the fundamental particles is currently understood to be generated from the 
spontaneous breaking of ele ctrow eak symmetry via the Higgs mechanism. The high 
value of the top quark mas provides a means to establish the Higgs mass scale 
because the two can be related via radiative corrections to the W mass, Myy. One 
loop corrections to Mw give M w — si n?i^(\-i±r) wnere includes a correction 

j j- I 3GFmf \ j xt> / \\GfM% cos 2 Ow i M% \ 

depending on m t (_ 7 _- 3 _jl_) and on the Higgs mass ( 

The fact that m t is 35 times higher than the next heaviest fundamental fermion 
gives it a place of importance for another reason. Given the Higgs vacuum expecta- 
tion value, v ~ 247 GeV, the current world average m t points to a top quark Yukawa 
coupling of Y t = 2*1^1 ~ l. 

Since the electroweak theory makes no prediction of 
these couplings and they range over many orders of magnitude, it seems surprising 
that one of them would be unity. This value raises interesting questions about the 
precise nature of the Higgs mechanism, and some believe that the measured value 
of mt may i ndicate extra Higgs doublets, as in Minimal Supersymmetric Models 
For these reasons, the measurement of the top quark mass has rapidly 
become one of the most important efforts in the whole top quark sector. It is in fact 
one of the key parameters to come out of the current Tevatron program. 



6.2. General techniques for mass extraction 

The top quark mass can be extracted from each of the broad channels in the ti 
final state: dileptons, £+jets and all-jets. Across these channels, there are three very 
general approaches to experimental measurement of m t . In order of increasingly 
sophisticated usage of the information in each event, they are: 

• Fitting Kinematic and Decay Parameters: The value of mt is reflected in 
several individual observables of top quark events. For instance, the direct t — > b 
decay means that in the top quark rest frame, the b— quark has momentum 
~ m t /2. In the lab frame, the motion of the top quark and jet specific effects 
reduce the strength of this correlation. However, the observed px spectrum 
of 6-jets is still anticipated to strongly follow m t . In addition, the decay length 
distribution associated with 6-jets also correlates with m t . By fitting one or more 
of these parameters one can discriminate between top quarks of different masses. 
Such a fit in a sample of data can be performed to pull out a measurement of 
m t . 

• Kinematic Reconstruction: ti final states leave as many as 18 kinematic 
observables (17 for €+jets, 14 for dilepton). Because of the specific decay chain 
from tt, there are correlations among the kinematics of these objects. For in- 
stance, two light-quark jets from a W boson should be consistent with a dijet 
mass ~ Mw- Constraint equations describe these relations and can be used to 
calculate solutions that provide a mass estimator for each candidate event. By 
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fitting these mass estimators for a set of events, the mass of the top quark is 
extracted. 

• Matrix Element Fitting: A special case of kinematic recontruction involves 
the full use of the information about top quark production and decay. By us- 
ing the leading-order matrix elements, in conjunction with a full knowledge of 
the experimental resolutions of the final state object momenta, a fit can be 
performed to the data. This is used to provide a probability that an observed 
event configuration is consistent with a top quark of a certain mass. 

Almost all analyses are based on the latter two strategies. These approaches 
generally yield substantially better uncertainty than the first because they use many, 
if not all, of the measurements made for each event. The matrix element approach 
in particular tends to extract the maximum information and so has slightly better 
performance. Because each method relies on different information from the complex 
tt event, both CDF and D0 pursue multiple strategies for each class of final state. 
By ascertaining the correlations between these approaches, they provide the best 
overall estimate of the top quark mass. 

All analyses that measure the top quark mass generate pseudo-experiments to 
test the performance and calibration of their methods and to determine the effects 
of systematic uncertainties. For these tests, simulated event samples {ensembles) 
are generated that correspond to the observed event sample in number of events 
and signal and background composition. These ensembles are then processed in the 
same way as the collider data. In this way the experiment at hand can be simulated 
many times. The mean (m) = J2iLi rn i/^ °f tnc measured top quark masses rm 
from N ensembles, their rms cr(m) = \J (( m ) — n*i) 2 /N , and the mean and 

rms of their pulls d = (m; — m t )/5mi can be determined. Here m t is the top 
quark mass assumed in the simulation and Srrii is the statistical uncertainty in the 
measurement from the i th ensemble. A well calibrated analysis method gives (m) = 
m t , (d) = 0, and a(d) = 1. Ensemble tests can be used to calibrate the methods, i.e., 
determine corrections to be applied to the measurements so that these relations are 
satisfied. Various sources of systematic uncertainties on the mass measurement are 
considered. Some of the most prominent sources which are carefully evaluated are: 
1) jet energy scale, 2) initial state gluon radiation, 3) final state gluon radiation, 4) 
parton distribution functions, 5) Monte Carlo generators, 6) background model, 7) 
^-tagging, 8) Monte Carlo statistics. 

6.3. The dilepton final state 

It is important to study m t in all of the final states in which the top quark can 
be identified for two main reasons. Non-standard decays can impact specific final 
states differently The kinematic analysis involved in a mass measurement therefore 
provides an additional test of the hypothesis that the signal events conform to the 
t — > Wb decay chain. If this chain is other than expected, sufficient statistical preci- 
sion would reveal discrepancies between the measured m t as estimated in different 
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channels. Dilepton channels provide both an independent statistical sample with 
which to measure m t . As the integrated luminosity at the Tevatron increases or 
the LHC turns on, the statistical limitations of this channel become less relevant. 
Systematic uncertainties are similar in magnitude to those for the single lepton chan- 
nels, and so the dilepton-based measurement has an important role in improving 
the world-average top quark mass. 

6.3.1. Fitting with kinematic parameters 

The effort to extract a measure of mt in dilepton events has involved each of the 
techniques mentioned in Section 16.21 The most basic of these involves fitting to 
kinematic distributions in candidate events. Because the top quark is so much more 
massive than any of the fundamental fermions in its decay chain, the magnitude 
of mt is directly manifested in the momenta of these particles. For instance, the 
transverse momenta of the 6-quarks are of order rrit/2. 

This approach has been tried by CDF using eight events in their explicit dilepton 
sample. Two different approaches were combined. If one ignores 6-quark and lepton 
masses, then the top quark mass can be expressed approximately as raj = + 
i- CO sV ib • O ne can rewrite this expression in terms of well-measured quantities, 

m\ = (Mf b ) + yfM*, + 4M*,(M? b ) + (M? b ) 2 (7) 

Here the b refers to the jets which are presumed to arise from the 6-quarks. There 
is some confusion arising from the combinatorics of mapping the jets and leptons 
that should come from the same top quark. Considering the allowed jet-lepton 
configurations, each event produces two Mu, measurements. These values are used 
to provide a value for (Mjfy for the sample using a prescription tuned with ti events 
generated with Herwig. A top mass measurement is extracted from this parameter 
in the data sampkil^l This is augmented with another approach which uses the 
average pr of the two b— jets in the event. This average is very correlated with mt/2. 
A maximum likelihood fit is applied to templates generated from this parameter for 
signal and background samples. A combined mas s estimate of 161 ± 17(stat) ± 
10(sys) GeV was obtained from these analyses^3 

6.3.2. Fitting using kinematic reconstruction 

While the technique of fitting to kinematic parameters provides useful, general esti- 
mates of m t , it does not use all of the information in an observed event that might 
maximize the measurement's precision. For instance, momentum correlations be- 
tween final state particles which can be used to test whether the event is consistent 
with two top quarks of equal mass, or two W bosons with mass equal to Mw- These 
capabilities are gained with an explicit kinematic reconstruction of an event. Per- 
forming such a reconstruction with dilepton events has the challenge that there are 
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two neutrinos. There are 18 unknowns for the six initial fermion three-momenta, 
but only four full three-momenta measured plus two components for the $t- Three 
additional constraints arise from requiring each Iv pair to give Myy, and both top 
quarks must have equal mass. As a result, a -1C underconstrained fit results and 
some information must generally be provided as input to solve these events. Gener- 
ally, approaches assume an input value of nit in steps over the allowed range (e.g. 
150 GeV to 200 GeV). 

The specifics of how the events are solved varies with method (see below). A two- 
fold ambiguity results for the solutions for the momentum of each neutrino. This 
ambiguity in kinematic reconstruction is compounded by the different assignments 
of lepton and jet to one of the t — > Wb decay chains. A particular configuration 
refers to a particular jet assignment in tandem with one of the neutrino solutions 
given that assignment. Generally, a weight is calculated for each configuration which 
reflects the relative consistency of that configuration with some additional observed 
property in the event. This weight is summed for all configurations which are solved 
for a particular input nit- 

This integration is performed for a range of assumed top quark masses to provide 
a probability or 'weight' vs. nrit for each event. Several parameters of the weight dis- 
tribution carry information about the actual m t . These 'mass estimator' parameters 
are used to compare events in data with expectations from top quark signal events 
of different masses, and background. Often, analyses will choose just the mass for 
which the weight distribution is maximum ('peak') as the mass estimator. The peak 
does indeed carry most of the sensitivity to the actual nit for the two kinematic 
reconstruction strategies described below. However, more information is available 
in the rest of the weight distribution. Other approaches extract this information by 
forming a coarsely binned template from the weight distribution, or taking instead 
the first couple moments of this distribution. In any of these cases, the distribution 
of mass estimators is generated for data, and this must be compared via maximum 
likelihood fit to templates of distributions for signal and background samples. 

6.3.3. Neutrino weighting 

One method for solving tt events utilizes the expected neutrino rapidity distribu- 
tions. This metho d is termed 'neutrino weighting' (uWT) and was first developed by 
D0 in Run fEU. A kinematic fit is performed by omitting two measured variables, 
Jfi x and $ y . To compensate for this loss of information, one assumes rapidities for the 
two neutrinos. An integration over the neutrino rapidity distributions is performed 
such that for each choice of neutrino rapidity pair and jet combination, the neutrino 
momentum solutions are calculated. For each configuration, a weight is calculated 
by comparing the measured event $t with the kinematically reconstructed $t 

W = T, «p { J «p [ 2a%y ) (8) 
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This weight is summed for all configurations for each value of m t . Templates using 
the parameters from the full integrated weight vs. mt distribution can be used to 
extract a top quark mass. 

D0 has performed this analysis with Run iUHl an d R un j J135 | d a t a by coarsely 
binning the weight distribution. The Run II analysis differs primarily in that the 
£+track channels were included. A maximum likelihood fit of these templates is 
performed to analogous templates for simulated signal and background, and instru- 
mental backgrounds from data. A probability density estimator approach is taken 
| 139 | xhe signal probability for a given m f , f s (W\mt), and background probability, 
fb(W) are established by comparison of data templates (W) with those for signal 
and background. The value of m t can be extracted by maximizing the likelihood 

L(m t ,n b ,n) = G(n b - n b ,a)P(n s + n b ,n) x II" =1 [(n s / s +n b f b )/(n s + n b )} (9) 

where a Gaussian and Poisson constraint are placed on the estimate of n b and 
n s + n b given the observed event yield, respectively. As is typically done with these 
analyses, the method is tested by performing the analysis on independent pseudo- 
experiments of simulated signal and background events. In common with most of the 
mass analyses, the relation of average fitted top quark mass vs. input nit was verified 
in each channel to have slope near unity and with small offset. This helps to keep 
systematic uncertainties from the method to a minimum. Any residual difference 
of the slope and offset from unity and zero, respectively, is used to correct the 
fitted mass. The statistical uncertainty from a given ensemble is extracted from 
limits in mt when —ln(L) varies by 0.5 from its minimum value. It is cross-checked 
by verifying that pull widths are near 1.0. Statistical uncertainties for the Run II 
analysis were corrected for a small deviation of pull widths from unity. The plot 
of the —ln(L) vs. m t resulting from the maximum likelihood fit to the Run II 
data is given in Fig. [27] The top mass was measured to be m t = 179.5 ± 7 A(stat) 
GeV. A 5.6 GeV systematic uncertainty is dominated by the jet energy scale. This 
measurement was combined with the matrix-weighti ng d etermination (see Section 
[673T4jl to yield m t = 178.1 ± 6.7(stat)± 4.8(sys) Ge\fl2ll The Run I measurement 
was mt = 170.0 ± U.8(stat) GeV 0^1 

The dilepton analys is with y WT has been pursued in preliminary analyses in 
the 1 fb _1 data sample &37 | 136 | Q ne q ues y on f or these analyses has been what is 
the most optimal way to use the output of the weight calculation (i.e. the weight 
distribution) to get the smallest expected uncertainty in the top mass. Ref. 11371 
presents an analysis in eu, events of three different approaches to extract the value of 
mt from the weight distribution. These analyses attempted to determine an optimal 
use of the most important parts of a weight distribution. One of these three methods 
constituted an optimized version of the binned template approach described above. 
Another approach used just the value of the top mass for which the weight was 
maximized, and then performed a 2-dimensional fit to this peak value and the 
input mt- A third approach considered various properties of the weight distribution 
(e.g. moments, integral weight in high and low mass bins, etc.) and settled on 
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the first two moments as the minimal number of parameters that provided the 
20% improvement in sensitivity usually associated with the use of the full binned 
template. The 2-dimensional fit and moments approaches each were found to match 
the apriori expected sensitivity of the binned-template method, and were combined 
into a measurement using dilepton events (ee, e/z, Hfj.) in 1 fb" 1 Ei A value of 
m t = 172.5 ± 5.8(stat) ± 3.5(sys) GeV was obtained. 
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Fig. 26. ^VKT results from the 359 pb" 1 CDF Run II ^+track sample. The l inear ity of the fit 
(left) and data overlaid with signal and background templates (right) are showrffSSl 



CDF has pursued the vWT approach in 359 pb 1 of data in Run II. The basic 
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kinematic reconstruction is performed similarly to that used by D0 . In Run II, they 
applied the technique to their £+track sample of 46 events. CDF chose the mass 
estimator as the peak of the weight distribution. Thus the data m pea k distribu- 
tion is compared to the expectation from signal and background (see Fig. 1261) . The 
maximum likelihood fit on data gives m t = 170.7^g'g(stai) ±4.6(sys) GeVl^OI. Sys- 
tematic uncertainties were estimated for jet energy scale, gluon radiation (ISR and 
FSR), background template shape, pdfs, and showering Monte Carlo generator. A 
further improvement for CDF is obtained by combining this analysis with two others 
of lesser statistical power. These use an assumption of the longitudinal momentum 
of the top quarks, or the azimuthal angle of the neutrinos, to solve each candidate 
eve nt. T he final result from this combination is nit — 170.1 ± 6.0(stat) ± A.l(sys) 

In Run I, which used an explicit dilepton selection, the mass estimator for 
each event was taken as the weighted mean in a window around the mass with high- 
est weigh t. This approach generated a value of m t — 167 .4 ± 10.3(stat) ± i.8(sys) 

CDF performed a consistency check with the fitting procedure garnered in 
single lepton double-tagged events by treating the two light-quark jets from the W 
has a lepton and a neutrino. By mimicing the dilepton signature in this way they 
obtain a mass estimate in these events which was within the expected resolution. 



6.3.4. Matrix weighting 

A nother ap proach to kinematic reconstruction of top quark events has been pro- 
posed E2H3H jjj this method, the measured $t is not omitted from the kinematic 
reconstruction, thereby permitting a tt event to be solved to an eightfold ambiguity 
if m t is assumed. Unlike the uWT approach, both the t and i are solved simultane- 
ously A weight is calculated for each solved configuration with an assumed mass, 
m t : 

w = f{x)f(x)P(E n \m t )P(E l2 \m t ) (10) 

where fix) is the pdf taken at the Feynman x obtained for the solution. The 
P(E[\m t ) are the probabilities for a lepton to have energy Ei in the solved top 
quark rest frame. This approach is termed 'matrix weighting' (MWT). 

The D0 Collaboration has pursued this approach in both Run I and Run II. 
Using the 3 70 p b" 1 data sample from Run II, they considered the explicit dilep- 
ton channel The candidate sample was broken down into the low background 
6-tagged subsample, and the higher background sample with strictly topological 
selection. This topological selection loosens kinematic and electron identification 
requirements on the e/i channel relative to what was done for the vWT analysis 
described above. This separation by 6— tagging was done to enhance the sensitivity 
to the signal events in data. The maximum likelihood is performed by taking the 
peak of the weight distribution: 



L{m t ) = ((n s Si(m t ) + n b bi)/(n s + n b )) 



(11) 
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Fig. 27. The m pea k distribution from the MWT technique applied to 370 pb of dilcpton 
events from D0 (left). The vWT (open circles) and MW T (cl osed circles) fits are shown (right) 
for the data. The former also includes the £+track channeii^lU 



where the product is taken over the i bins of the m pea k distribution from data (rij) as 
compared to a signal template (sj) and a background template This constitutes 
a different approach than used for Run I since there are single templates from the 
Wpeafe distributions of the data, signal and background rather than templates from 
full weight distributions for each event in data, signal and background samples. 
From 26 dilepton events, a value of m t = 176.2 ± 9.2(stat) ± 3.9(sys) GeV was 
determined. The plot of the —ln(L) vs. m t is given in Fig. 1271 For comparison, 
the Run 1 measurement was m t = 168.2 ± Y2A(stat) Ge W^t For the Run 11 
measurement, the dominant systematic uncertainty was from the uncertainty in jet 
energy calibration. The matrix-weighting approach has been applied to the 1 fb" 1 
sample in dilepton events and yielded m t = 175.2 ± Q.l(stat) ± 3A(sys) GeV^^l 
6.3.5. Matrix element fit 

Use of the full matrix-element to encompass the expected correlations of particle 
kinematics in a mass analysis was first applied by D0 in complete form with sin- 
gle lepton eve nts^^ l The CDF collaboration has performed such an analysis with 
dilepton events^22l. The analysis uses the explicit dilepton event selection from the 
cross section measurement described above. This gives 33 events in 340 pb -1 . 

This analysis differs from those described above which integrate over just one 
class of parameters to perform the fit, such as neutrino rapidity or Feynman x. 
Instead the leading order matrix element for qq — > ti production is used. Some initial 
assumptions are implemented which permit a more tractable calculation. Lepton 
momenta are assumed to be perfectly measured. Quark directions are assumed to 
be perfectly indicated by observed jet angles. An integration is then performed over 
the neutrino and quark energies as calculated from the LO matrix element. This 
integration is constrained by a transfer function which expresses the mapping of 
quark energies to observed jet energies obtained from a sample created with Herwig 
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fed into GEANT. The probability density can then be expressed as 



P s (x\m t ) 



1 



/ 



d$ \M t t{qu Pi;m t )\ 2 x J| W(p l ,j l )fp DF {q 1 )f PDF (q 2 ) (12) 



tr(m t ) 



jets 



where the matrix element is performed as in Ref . ' I When performing this integra- 
tion, the contribution from the leading order gluon-gluon production is omitted. A 
similar probability density calculation is performed for those background processes 
arising from Z and WW production, as well as W + jets in which one jet fakes a 
lepton. One can then obtain a joint likelihood for an event sample by multiplying 
the likelihoods for each event from the background and signal probability densities. 

In order to test the validity of the result, ensembles are constructed from simu- 
lated signal and background samples which were allowed to fluctuate with Poisson 
fluctuations. The instrumental background was obtained from data. The resulting 
slope vs. input (pole) nit, shown in Fig. 1281 is not unity primarily because of the sim- 
plification of the background components in the probability density calculation. The 
measured top quark mass and estimated statistical uncertainty are corrected for this 
slope. Additionally, the resulting pulls are 1.51, so the statistical uncertainty is cor- 
rect ed fo r this. The resulting top mass measurement in data is m t — 165.2±6.1(si<rf) 
Figure. [28l shows the probability density vs. m t , as well as the —ln(L) vs. 
mt (inset). Systematic uncertainties are estimated for several sources. The chief one 
is the jet energy scale (2.6 GeV) which is obtained by varying this calibration within 
its allowed range. The effect of uncertainties in pdfs is quantified in several ways, for 
instance comparing CTEQ5L and MRST72, and by exploring the range given in the 
CTEQ6M set of functions. Combining all sources gives a 1.0 GeV pdf uncertainty. 
Initial and final state radiation are varied giving a 0.5 GeV and 0.7 GeV uncertainty, 
respectively. Pythia and Herwig generated events are used to provide a Monte 
Carlo generator uncertainty (0.8 GeV). Other uncertainties are calculated for back- 
ground modeling and the fitted mass slope correction. The total uncertainty tallies 
to 3.4 GeV. When combined with the template-based methods, the value of m t was 
determined to be 167 .9 ±5.2(stat) ±3.7 (sys) GeV. The matrix element metho d was 
also applied in 1.03 ftr 1 of data and yielded 164.5 ± 3.9{stat) ± 3.9{sys) GeV^^. 

6.4. Single lepton channels 

£+jets events the top quark mass is reconstructed from the momenta of the charged 
lepton (e or //), the transverse momentum of the neutrino which is inferred from 
the $t in the event and the momenta of the four jets in the detector, two from 
the fragmentation of the b quarks and two from the hadronic W decay. There are 
two unknown parameters in the kinematic description of lepton+jets events: j> z and 
m t . On the other hand the event kinematics have to satisfy four constraints. The 
charged lepton and neutrino momenta must add up to a system whose mass equals 
the mass of the W boson. Adding the momentum of the 6-jet from the semi-leptonic 
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Input M, [GeV] M, [GeV] 

Fig. 28. CDF matrix element results in 340 pb — 1 dilepton events. The scale between fitted and 
input top mass is shown (left) and the probability density as a function of ni( is shown at righJESl 

top quark decay must give a system with the mass of the top quark. The invariant 
mass of the di-jet system from the hadronic W decay must equal the W boson 
mass and the three-jet system from the hadronic top quark decay must have the 
mass of the top quark. We can thus perform a 2-C kinematic fit of the £+jets decay 
hypothesis to these events to determine the most likely value of the top quark mass. 
This simple picture is complicated by the presence of jets from gluons radiated from 
the initial state and by the splitting of jets from the top quark decay. Moreover it is 
not possible to uniquely assign the jets in the event to the partons from tt decay. If 
only the four jets with the largest transverse momentum values are considered, there 
are 12 different ways to assign jets to partons. If one jet is identified as originating 
from the fragmentation of a 6-quark and is only assigned one of the 6-quarks this 
reduces to eight permutations. If two jets are identified as 6-jets, this further reduces 
to two permutations. 

The experiments have applied several techniques to extract the top quark mass 
from ^+jets events. Almost all are based on the kinematic reconstruction of the top 
quark mass. An exception is a measurement by the CDF collaboration that used the 
decay length distribution of b— jets to determine the top quark mass. The simplest 
kinematic reconstruction technique is the template method which has been used by 
CDF and D0. The matrix element method and the ideogram method were both 
developed at D0 and are now used by both experiments. We will describe these 
methods in the following sections. 

The kinematic reconstruction of the top quark mass requires that all energy and 
momentum measurements are well calibrated. The electron energy scale and the 
muon momentum scale can be calibrated precisely enough with Z — > l + t~ decays 
that the residual uncertainty is negligible. The calibration of the jet energy scale, 
however, is more difficult and gives rise to the dominant systematic uncertainty in 
the top quark mass measurement. In addition to the apriori calibration of the jet 
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energy scale, most recent measurements of the top quark mass in the ^+jets channel 
also make use of the hadronic W boson decays in these events to gain an additional 
constraint on the jet energy scale. These analyses use the known W mass as an 
input parameter and then determine simultaneously the top quark mass mt and an 
overall jet energy scale parameter aj es that multiplies all jet energies. 

6.4.1. Template method 

For the template method one chooses an estimator for the top quark mass, typically 
the best fit mass mfn from a kinematic fit of the event. A simulation is used to 
compute the expected distribution of this estimator based on the expected number 
of signal and background events as a function of the top quark mass. These distri- 
butions are referred to as templates. One then performs a fit of the templates to 
the distribution of the estimator observed in the data to determine the value of the 
top quark mass that best predicts the data distribution. 

The early measurements of the top quark mass by CDF and D0 in 199fP^used 
this method. D0 has used this method for its first t op qu ark mass measurement 
from Run II using of approximately 230 pb _1 of dat a ^1 , CDF also performed a 
measu reme nt of the top quark mass using this method on 318 pb" 1 of data from 
Run iiEH. 

The CDF analysis divides the sample of events with a high-p-r electron or muon, 
$t and at least four high-p-r jets into four subsamples of different combinatorics 
and background contamination levels defined by the number of 6-tagged jets: at 
least two 6-tagged jets, at least one 'tight' 6-tagged jet, at least one 'loose' 6-tagged 
jet, and no fe-tagged jets. For every event with x 2 < 9 for the kinematic fit the 
hypothesized top quark mass nifit that minimizes x 2 is entered into a histogram. 
The masses of all pairs of jets that are not tagged as &-jets are entered in another 
histogram. Both quantities are histogrammed separately for the four subsamples. 

The templates are histograms of mfu and rrijj from Monte Carlo simulations, 
parameterized in terms of the hypothesized top quark mass m< and the jet energy 
scale parameters ctj es and separately for each of the subsamples. Sample templates 
are shown in Fig. [29j 

A simultaneous fit of the data to the parameterized templates gives best agree- 
ment for m t — 173.5^3 g(stat © jes) ± 1.3(sys) GeV. The apriori calibration of the 
jet energy scale a ± 8a enters into this fit as a Gaussian constraint on Aa = ± 1, 
where Aa = (oto — aj es )/6a. Best agreement is achieved for Aa = — O-lOj^'g 8 . 
Figure [30] shows a contour plot of the likelihood in the m t — Aa plane. 

An identical analysis is carried out on a fac tor of two larger data sample 
corresponding to a luminosity of 680 pb _1 H3S1, leading to a measurement of 
m t = 173.4 ± 2.5(stat © jes) ± 1.3(sys) GeV. 

A variation of this analysis which uses a 2-dimensional templates of top quark 
mass and hadronic W boson mass constructed using a Kernel Density Estimate 
(KDE) has been carried out by the CDF collaboration. With this technique and a 
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Fig. 29. Sample templates for mf it and rrijj for the CDF ^+jets dataEU. 
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data sample with integrated luminosity of 1.7 fb 1 the top quark mass is measured 
to be m t = 171.6 ± 2.1(stat ® jes) ± l.l(sys) GevQUl 

6.4.2. Matrix element method 

The matrix element method was developed by D0 between Run I and Run II. It was 
used first on the Run I data s et to obtain a top quark mass measurement of mt = 
180.1±3.6(stat)±3.9(sys) Ge\fl31. A significantly improved statistical precision was 
achi eved compared to the previous measurement mt = 173.3 ± 5.6(stat) ± 5.5(sys) 
on the same data set which used the template technique. 
The idea is to compute for each event the probability density to observe the 
event as a function of the top quark mass, based on the full kinematic information 
from the event. This probability density is given by 

Pevt(x\m t , a, /) = fptop(x\m t , a) + (1 - f)p bkg (x\a). (13) 

Here x stands for all measured quantities in the event, i.e. the momenta of the 
charged lepton, the jets and the $t- a is a scale parameter for the jet energies, and 
/ represents the fraction of top quark decay events in the data sample. The event- 
by-event probability densities are combined into a joint likelihood for all events in 
the sample, 

n 

-\nL(x 1 ...x n \m t ,a, f) = -^lnp evt (x\m t , a, /). (14) 

t=l 

The measurement of the top quark mass is then the value of mt that maximizes L 
for any value of the parameters a and /. 

The signal probability density is given by the differential cross section normalized 
to the total cross section a t j(m t ) for all events accepted by the analysis cuts. It is 
given by 

p sig (x\mt,ct) = — j — - f dzdzf(z)f(z)do- t j(y,m t )W(x\y 7 a), (15) 

where z and ~z are the fractions of the proton and antiproton momenta carried 
by the initial partons, f(z) is the parton distribution function for the proton, y 
represents the momenta of all partons taking part in the hard scatter event. The 
transfer function W(x\y,a) gives the probability to measure the observables x for 
the parton momenta y and the jet scale parameter a. The parton cross section 
da t j is calculated based on the leading-order matrix element M. for the process 
qq — ► tt —>■ tvbqq'b: 

\M\ 2 

da t - t (y,m t ) = l -^d<S> e , (16) 

xxs 

where d$6 is the Lorentz-invariant six-particle phase space element and s is the 
parton center of mass energy. In order to make p S i g calculable with reasonable 
computing power, a number of simplifying assumptions must be made. Typically 
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these are assuming that pr(it) = 0, that all angles are well measured and that the 
transfer function factorizes into contributions from each final state object. After 
these assumptions the remaining six integrations can be carried out using Monte 
Carlo integration techniques. The calculation is carried out for all possible jet-parton 
assignments and for both solutions of p z of the neutrino. These contributions are 
then added together for the final value of p S i g . A similar calculation is carried out 
for pbfcg except that this quantity does not depend on the top quark mass. 

The power of this method originates from the use of all kinematic information 
from the events and the use of signal and background probabilities for each event. 
This effectively weights events in the event likelihood L that are more likely to be 
from top quark decay stronger than events that are more likely to be background. 
The downside is that the phase space integration is very computationally intensive. 

The D0 Collaboration published a measurement of the top quar k ma ss using 
the matrix element method based on 370 pb _1 of data from Run In this 

measurement D0 makes use of 6-tagging to give higher weight to jet permutations 
that assign 6-tagged jets to 6-quarks. This leads to a reduction in the combinatoric 
background and improves the statistical precision of the measurement. The result of 
this measurement is m t = 170.3^ 4 g(stat® jes)^'g(sys) GeV. The measurement was 
updated using a data sample corresponding to an integrated luminosity of 1 fb~ 1 H4SI. 
Figure [3T1 shows a contour plot of the event likelihood L defined in equation [T4l as a 
function of the assumed top quark mass nit and the jet energy scale factor a. Plots in 
Fig. [32] show the projections onto the two axes. The updated measurement applied 
to fe-tagged events yield a top quark mass of m t — 170.5 ± 2.4(stat jes) ± 1.2(sys) 
GeV. The first uncertainty is derived from the width of the likelihood and accounts 
for both the statistical uncertainty and the uncertainty from the jet energy scale 
calibration. Ensemble tests show that the matrix element algorithm results in a 
measured value of the top quark mass that tracks the input top quark mass with a 
small offset. Figure [33] shows a plot of measured top quark mass versus input top 
quark mass for this analysis. The result is corrected for the observed offset. 

An analysis usin g the same technique is carried out by CDF using a sample 
of 940 pb _1 of data 1 149 1 phe events are required to have a lepton, $t and only 
four jets, out of which one is expected to be tagged as a 6-jet. This sample yields a 
measurement of mt = 170.9 ± 2.2(stat ©jes) ± 1.3(sys) GeV. 

By combining the signal probability computed using the matrix element tech- 
nique with a Neural network discriminant to reject backgrounds, CDF experiment 
has measured mt = 170.9 ± 1.3(stat) ± 1.2(jes) ± 1.2(sys) GeV This measurement 
is based on a data sample with integrated luminosity of 1.7 fb _1 H211. 

6.4.3. Quantized dynamical likelihood method 

CDF has developed another method, based on Ref. 0313 called the quantized dy- 
namical likelihood method. In this analysis only events with exactly four jets are 
used. The four jets are assigned to partons in every possible way with the restriction 
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Fig. 31. Contour plot of the likelihood for the matrix element analysis of the DO e+jets data (left 
panel) and /i+jets data (right panel) LHH1. 
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Fig. 32. Likelihood versus top quark mass (left), and jet energy scale (right) f or th e matrix 
element analysis of the D0 e+jets data (top plots) and /j+jets data (bottom plots) H3S1 



that 6-tagged jets are only assigned to one of the b quarks. 
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Fig. 33. Measured top quark mass from ensemble tests versus inpu t top quark mass for matrix 
element analysis of the D0 e+jets data (left) and /x+jets data (right j ^-^H 



For each permutation a likelihood is computed as 

da 

Liymt) = — , 
d$ 6 



(17) 



where d$6 is the Lorentz-invariant six particle phase space element and da is the 
differential tt production cross section integrated over the momentum fractions z 
and z of the initial state partons and the transverse momentum pr(tt) of the tt 
system: 



da = 



dzdzf(z)f{z)d 2 p T (tt) f (p T (tt))da tJ (x,m t ) 



(18) 



where f(pT{tt)) is the expected px distribution of the tt system, da t j(x,mt) is the 
parton cross section defined in section 16.4.21 and all other parameters are also as 
defined in the same section. 

The likelihood L\ is typically averaged over 50,000 random generations of the 
parton momenta x from the observed jet momenta using Monte Carlo derived trans- 
fer functions. These functions are taken over all possible jet-parton assignment and 
over the two possible values of the z-component of the neutrino momentum to ob- 
tain the event likelihood L(m t ). This joint likelihood for all events is computed 
by multiplying together all the event likelihoods. The top quark mass is measured 
as the hypothesized top quark mass for which the joint likelihood is maximized. 
Finally, this result is corrected for the effects of backgrounds in the sample. This 
correction is determined from ensemble tests using Monte Carlo data. 

Using a data sample of 63 tt candidates in the lepton+jets channel with at least 
one 6-tagged jet from 318 pb _1 of data, CDF measures m t — 173.2^2 4(stat) ± 
3.2(sys) GeV. 



6.4.4. Ideogram method 

The ideogram methooP^] j s based o n ide as that were used by the Delphi Collabo- 
ration to measure the W boson mas d^*^ . It was first adapted for the measurement 
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of the top quark mass by the D0 Collaboration. It is based on the same general 
idea that was developed for the matrix element method in equations [IS] and [H] 
The difference lies in the definitions of p S i g and Pbkg- In order to reduce comput- 
ing requirements, this analysis method makes use of the same kinematic fit as the 
template method. 
The signal probability is defined as 



p S ig(o\m t ,a) =p sig (D) 



24 

E 



exp 



/ / G(mi,m',ai)B(m',mt)dm' + (1 - f)S(mi,m t ) 

(19) 



Here D is a discriminant based on the kinematic observables in the event that is 
constructed such that most top quark like events have D « 1 and most background 
like events have D ss and p S i g (D) is the probability density for signal events with 
discriminant value D. The sum over i runs over all 12 jet permutations and the two 
solutions for p z of the neutrino, xl 1S the minimum value of the x 2 goodness of fit 
variable for the kinematic fit described in section [6 . 4 . 1 1 and m, is the corresponding 
value of the hypothesized top quark mass. G(x, xo, a) is a Gaussian bell curve with 
mean xq and width a and B{x,xq) is a Breit-Wigner function of mean x$. The 
integral is a convolution of the Breit-Wigner line shape of the top quark with a 
Gaussian resolution function and represents the contribution from jet permutations 
with the correct jet-parton assignments and / indicates the probability that a jet 
permutation corresponds to the correct jet-parton assignments. S(nii,mt) repre- 
sents the contribution of jet permutations with the wrong jet-parton assignments. 
Figure [34] shows the event likelihood curves for simulated events with zero, one, 
or two tags. The background probability density pbkg is determined using a Monte 
Carlo simulation. 
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Fig. 34. Event likelihoods for simulated D0 £+jets events with 0, 1, or 2 tags322 



The D0 collaboration has applied this technique to the same data set as was used 
by the matrix element analysis described in section 16.4.21 except that the ideogram 
analysis uses events with four or more jets. Figure[35] (left plot) shows a contour plot 
of the event likelihood L defined in equation [14] as a function of the assumed top 
quark mass m t and the jet energy scale factor a. Figure [35] (right plot) shows the 
likelihood versus top quark mass. The result of this measurement is m t = 173.7 ± 
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4.4(stat jes)l2'o( s y s ) GeV. The first uncertainty is derived from the width of the 
likelihood and accounts for both the statistical uncertainty and the uncertainty from 
the jet energy scale calibration. Ensemble tests show that the ideogram algorithm 
results in a measured value of the top quark mass that tracks the input top quark 
mass with a small offset. The final quoted result includes this calibration. 
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Fig. 35. Contour plot of the likelihoo d (le ft - ! and Likelihood versus top quark mass for the 
ideogram analysis of the D0 £+jets dat aj [ 



6.4.5. Decay length distribution 

In addition to the measurement techniques described in the past few sections, the 
CDF collaboration has measured the top quark mass from the decay length distri- 
butions of 6-hadrons from top quarks decays. The mean momentum of the 6-quarks 
from the decays of top quarks increases with the top quark mass. A harder mo- 
mentum spectrum leads to a higher Lorentz-factor and a longer decay length in the 
lab frame. Thus the decay length spectrum of 6-hadrons from top quark decays can 
yield a measurement of the top quark mass, albeit with significantly poorer statis- 
tical sensitivity than the methods that are based on the kinematic reconstruction 
of the top quark decay. The advantage is that this method has different system- 
atic uncertainties. Most importantly, it does not depend strongly on the jet energy 
scale calibration. The result of this measurement is m t = 180.7±^(stat)±8.6(sys) 
GeV US! 



6.5. All-jets channels 

Events in which both W bosons decay to quark pairs provide another avenue with 
which to measure the top quark mass. The lack of a neutrino simplifies the kine- 
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matic reconstruction by removing quadratic ambiguities and providing for an over- 
constrained fit from the momenta of six jets. The large all-jets branching fraction 
provides a potentially substantial source of data that is essentially uncorrelated 
with the leptonic modes. On the other hand, a very high background persists even 
though fe— tagging is used. This substantially impacts the statistical power of this 
channel. The Tevatron experiments have pursued several analyses to measure m t in 
the all-jets channel: three template-based measurements, an ideogram analysis and 
a matrix element analysis. 

6.5.1. Template methods 

Templated app roach es for the all-jets mass extraction have been tr ied b y CDF in 
Run II 155 and £GEH, and by D0 using Run I data published in 2005^21 The most 
sensitive of these is the Run II CDF analysis which is a preliminary result that used 
1.02 fb _1 . All three use event kinematic selections very similar to those described 
in Section [4~5l At least six jets are required in all cases, and CDF relaxed its initial 
Y>Et cut to minimize the bias to the mass measurement. The D0 analysis used soft 
//-tagged all-jets events, while CDF used secondary vertex 6-tagged events. 

To model signal, top quark events were modeled by both experiments using 
Herwig. CDF considered masses from 150 GeV to 200 GeV in their Run II analysis. 
Backgrounds were estimated similarly by both experiments. An untagged multijct 
sample was selected just as the signal sample. Known mistag rates were then applied 
to properly weight events when constructing templates. Estimatation of signal and 
background levels gave 334 events and 573 events, respectively, for CDF's Run II 
sample. The Run I analyses had 65 and 136 events total for D0 and CDF, with 48 
and 108 expected background events, respectively. 




Fig. 36. The CDF measurement of the top mass in the all-jets final state using Run I data 
(leftJEEl. Expected background and top q uark signal are shown in solid, and data points provide 
data. At right is the D0 distribution 1 156 1 j nse t s at upper right of both plots provide the fit used 
to determine the mass. 
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Fig. 37. The result of the template analysis of 1 fb -1 of CDF Run II data is shown with signal and 
expected background components (left). The inset at upper right provid es the resulting —ln(L) 
vs. rat- The distribution of expected statistical uncertainties (right) ISSI j s a j so provided showing 
the result from data (see arrow). 



Kinematic reconstruction of the candidate events was performed. For this re- 
construction, each of the six highest Et jets was associated with one of the quarks 
from a top quark or W boson decay. Tagged jets were specifically associated with 
b— quarks. There is an ambiguity in how the jets are assigned in pairs to W bosons 
and how these are associated with b— jets to form each top quark. Resolving this 
situation involved considering all possible configurations and calculating a x 2 with 
respect to a correct tt configuration. Each configuration presents two jet triplets 
purported to originate from each top quark. For example, D0 defined a \ 2 f° r each 
configuration, 



2 _ / m tl - m t2 \ ( M Wl - M Wo \ f M W2 - M W( 
\ 2 x a„ H J \ a Mw J \ <?m w 

where Mw {= 77.5) GeV is the mass reconstructed in tt simulation for the two jets 
from a W boson. The values of My^ and m ti correspond to the masses calculated 
from the jet pair or triplet for the zth W boson or top quark, respectively, given the 
configuration being considered. All analyses resolved the combinatoric ambiguity 
by selecting the one with the lowest value of the \ 2 - Using this combination, the 
invariant masses of three jet groupings were calculated. They are shown in Fig. [36] 
for the Run I measurements, and Fig. [37] for the analysis in 1 fb -1 . 

This distribution is fit using expectations for backgrounds and for a top quark of 
varying mass and background. CDF in Run II has used a convolution of one gaussian 
and two gamma functions to represent the signal shape, and two gaussians plus a 
gamma function for background. The background shape was adjusted for a small 
estimated top quark contamination. CDF obtained measurements of mt = 174.0 ± 
2.2(stat) ± 4.8(sys) GeV and 186 ± 10 GeV in Run II and Run I, respectively. The 
primary systematic uncertainties in the former arose from the jet energy scalibration 
and the Monte Carlo modeling of the top quark signal. When applying the maximum 
likelihood fit, D0 allowed both signal and background levels to float. The fitted mass 
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is 178.5 ± 13.7 (stat) GeV. The systematic uncertainties for jet energy calibration, 
6-tag rate, and background statistics totalled 7.7 GeV. 

6.5.2. Matrix element analyses 

As with other channels, matrix ele ment tech niques have now been exploited in 
the all-jets channel, to date by CDF | 157 | 158[ rp ne f ormer j g ac t U ally a variant of 
the ideogram method used in ^+jets samples but also incorporating the matrix 
element approach to perform the m t determination. This preliminary result was 
obtained from 310 pb" 1 of 6-tagged all-jets data. Each event presents 90 potential 
configurations of grouping jets to form W's and top quarks. All combinations are 
considered when performing a kinematic fit. Goodness of fit and the probability 
that two jets are 6-jets were used to pick the right configuration. A two dimensional 
likelihood was performed in terms of the masses of the two top quarks, m t \ and rriti- 
The shape of the signal was obtained from Pythia and Herwig using Gaussian 
resolution functions for jets obtained from data. The QCD background was modeled 
with Alpgen. The likelihood yielded a measurement of mt = 177.1 ± 4.9(stat) ± 
4.7(sys) GeV. The leading systematic uncertainties were the jet energy calibration 
and the background shape. 
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Fig. 38. Contour plot of the favore d reg ion in the jet energy scale vs. rat for the CDF matrix- 
element analysis of 943 pb — 1 of data LU^left) . The distribution of fitted top mass for the double 
tagged sample is shown at right. Expected signal and background are shown. 

U sing a pure matrix-element approach, CDF has also generated a preliminary 
resuh V^ in 943 pb _1 which is an indicator of things to come. The event selection 
requires six jets with a scalar sum of their Et of more than 280 GeV. Spherically 
distributed events which tend to have central jets are kept if they have at least one 
6-tagged jet. Extraction of the top quark mass proceeds by first estimating a prob- 
ability vs. top quark mass using a matrix-element calculation. Transfer functions 
were used to relate reconstructed jet spectra to parton level momenta. This weight 
distribution is then used as a template and the data is compared to templates gen- 
erated from signal and background Monte Carlo. By relating the fitted dijet masses 
from W bosons according to each jet-quark combination chosen to the expectation 
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from W — > jj, an estimate can be made of the residual jet energy scale as well as the 
measured correction to the top quark mass. This relation is shown in Fig. 1381 This 
serves to validate the method and provide a systematic uncertainty for the residual 
jet energy scale. The fitted mass distribution for double-tagged candidate events is 
also shown in Fig. [35] The CDF analysis yielded a result of nit = 171.1 ± 4.3 GeV 
where 2.1 GeV of the uncertainty is due to systematic effects. 



6.6. Combined fits and electroweak constraints 

The Tevatron Electroweak Working Group^^l has performed combinations of the 
top quark mass measurements by D0 and CDF, properly taking into account corre- 
lations between the measurements in different channels and from both experiments. 
As o f Ma rch 2007, the most precise value for the top quark mass is 170.9 ±1.8 
GeV More details on the combined results are given in Table [TDJ Their impli- 
cations are discussed in Section 19.21 

The measured value of the top quark mass can be compared to the value of the 
top quark mass from the global electroweak fit performed by the LEP Electroweak 
Working Group^H. Based on the Z lineshape measurements a value of 172.6^10 2 
GeV is most consistent with the standard model. If one includes in addition the 
measured W boson mass and width into the fit, the best top quark mass goes up 
to 178. 9^g 1 6 7 GeV. Either value is completely consistent with the value from the 
direct measurement. 



Table 10. Summary of Top Quark Mass Measurements. 



Channel 


Experiment 


Measured Mass (GeV/c 2 ) 


It 


CDF 


164.5± 5.6 


ti 


D0 


172.5± 8.0 


£+jets 


CDF 


170.9± 2.5 


^+jets 


D0 


170.5± 2.7 


All-jets 


CDF 


171.1±4.3 


Lcpton+jets (L xy ) 


CDF 


183. 9± 15.8 



7. Properties of the Top Quark 

As with the other fundamental fermions, the ability to study the top quark 
presents new opportunities to test our current understanding. Measuring the charge 
of the top quark helps to confirm if the particle discovered in 1995 is in fact the 
weak isospin partner to the 6-quark. Determining V t b will either confirm or refute 
whether we live in a world confined to three quark doublets. Several properties of 
the top quark have been measured so far and the increasing sizes of data samples 
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are placing real constraints on fundamental physics. 
7.1. W boson helicity measurements 

The high mass of the top quark raises the question whether there are new interac- 
tions near the electroweak symmetry breaking energy scale. Measuring the helicity 
of the W boson from top quark decay probes the tWb vertex and provides a strin- 
gent test of the standard mod el. A general form of the Lagrangian for the tWb 
interaction is given in RefP^U It contains four form factors: /f, f± that param- 
eterize the V — A and V + A interactions, and f% which reflect the strength 
of an anomalous weak magnetic moment. In the SM, the coupling of W bosons to 
fcrmions is purely V — A, and therefore the only non zero form factor is . If the 
tWb couplings are standard, then top quarks decay to left-handed W bosons ( W_ ) 
or to longitudinal W bosons (Wo)- In the presence of non-standard couplings such 
as V + A, some admixture of right-handed W bosons (W+) is expected. 

7.1.1. Sensitive Variables 

Measurement of the W boson helicity can be performed using any leptonic top 
quark decay. CDF and D0 have analyzed both their dilepton and £+jets samples 
to extract limits on the W boson helicity using three sensitive variables. 

• the helicity angle 9*, defined as the angle between the charged lepton and the 
top quark directions in the W boson rest frame, 

• the transverse momentum of the lepton in the laboratory frame, 

• the invariant mass, Ma, of the charged lepton and b— jet thought to come from 
the same top quark. 

If we define 8* as the angle of the decay positron in the W boson rest frame, 
with the polarization axis defined by the direction of the W boson in the top rest 
frame, then the angular distribution of the lepton with respect to the polarization 
of the W boson is given by 

w(cos(9*)) = /-f (1 - cos(9*)) 2 + /o|(l - cos 2 (0* )) + + cos(6* )f 

where /_ , fo and /+ are the fractions of left-handed, longitudinal, and right-handed 
W bosons, respectively. In the SM /_, f and /+ are expected to be 0.7, 0.3 and 0, 
respectively. 

Computation of the helicity angle requires the reconstruction of the top quark 
and the pairing of the lepton with a jet in the event. Reconstruction of the top 
quark using semileptonic decays is performed using a constraint kinematic fit to the 
it hypothesis (see Section HT4|) . In the fit, Mw and n%t are fixed to their respective 
measured values. The fit is performed using all possible jet assignments, 12 in gen- 
eral, six if there is a single 6-tagged jet, or two if there are two b— tagged jets. The 
permutation which gives the lowest x 2 is chosen. Once the event is reconstructed, 
the helicity angle is computed after a boost to the rest frame of the reconstructed 
W boson. 
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When dilepton events are used, each event has two leptons and hence contributes 
twice to the measurement. Kinematic reconstruction of dilepton events is described 
in Sect ion |6~51 Here, the presence of the two neutrinos kinematically underconstrains 
the system. For an assumed top quark mass, the kinematics are solved algebraically 
with a four-fold ambiguity in addition to the two-fold ambiguity which arises from 
the pairing of the lepton with the jet. Once the detector resolution effects are folded 
into the measurement, an average value of the helicity angle for each lepton is 
available. 

Because the Ma and lepton px are directly determinable from the lab frame, 
using them greatly simplifies the analysis. They can also be readily applied in dilep- 
ton events as well as £+jets events where one of the four final state jets might not 
have been reconstructed. The Ma approach relies on t he fac t that the helicity angle 
cosine, cos(9*), can be approximated by the expressior l^^ 

M 2 

Each lepton-jet pairing proves a measurement of the angle. The lepton px method 
relies on the correlations in W boson momentum and that of its decay lepton. Left- 
handed W bosons tend to emit the charged lepton in the direction opposite to their 
direction of flight, longitudinal W bosons perpendicularly to their direction of flight, 
and right-handed W bosons emit the leptons preferentially along their direction of 
flight, leading to increasingly harder lepton pt spectra. However, this method is less 
powerful than measuring the helicity angle. 



7.1.2. CDF Results 

The CDF experiment has pursued measurements of the polarization fractions of 
the W boson via both the Ma and lepton pt measurements. These measurements 
have used three channels: dilepton, £+ jets with one 6— tagged jet, and €+jets with 
two tagged jets. In each case, templates are generated for tt decays with V + A and 
V — A couplings as well as background events. A likelihood calculation is used to 
extract constraints on the non-standard couplings. 

In an analysis published in 2005, the £+jets sample was used in y/s — 1.8 TeV 
colli sions to provide a measurement of the V + A decay rate at the tWb vertex 
HSZ1. They utilized the Ma method to place a limit on / + . To strengthen the limits, 
they combined these results with earlier measurements o f the W boson polarization 
using the lepton p T , f = 0.91 ± 0.37(stat) ± 0.13(systj^l A limit of /+ < 0.18 
was obtained at 95% C.L. 

This was followed with a measurement in y/s = 1.96 TeV collisions A 162 
pb _1 sample of £+jets events was used in both single and double tagged chan- 
nels. Both the Mil, and lepton pt measurements were employed for these channels. 
Dilepton events from 193 pb _1 were also used to obtain /o and / + via the lepton px 
technique. Each dilepton event provides two measurements for each event. These 
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samples are described in Section 14.31 and 14.41 The primary difference in £+ jets 
event selection involved an additional requirement of a fourth jet with px > 8 GeV 
and 1 77 1 < 2. This allowed the events to be kinematically reconstructed with the ti 
hypothesis. m t taken to be 175 GeV and the lepton-jet matching was performed 
based on the result of this fit. 

In order to extract a measurement, templates were generated for the different 
samples. In all cases, one W boson was simulated with SM expectations and the 
other was forced to a chosen helicity. Leptons from the decay of this second W 
boson were used to generate Ma, and p T templates for ft decays with V — A and 
V + A couplings and for background events. 

These templates were used in a fit to data. To extract fo, f + was fixed to 
zero, and to extract /+, fo was fixed to the standard model expectation of 0.7. 
Interestingly, the value obtained in the dilepton sample (fo = —0.54^° 35 ± 0.16) 
was unphysical because of the soft pt spectrum of the observed events. This result 
is consistent with fo = 0.7 only at the 0.5% level. It differs by 2 a from the £+ 
jets measurement of fo = 0.95^q'42 ^ 0.17. As indicated in section [6j the dilepton 
sample was analyzed and determined to be marginally consistent with the standard 
model. All results were combined to obtain fo = 0.74+Q 34 (stat + sys) and /+ = 
0.00i^(stat + sys). 

This gives a limit of /+ < 0.27 at 95% c.l. To account for correlations among 
uncertainties, Monte Carlo experiments were conducted to obtain the necessary cor- 
relation coefficient. Significant sources of systematic uncertainty came from back- 
ground modeling, the uncertainty in the value of mt, and the calibration of jet 
energies. 

In 700 pb _1 of Run II data, CDF has used t he M g technique in dilepton, £+ 
jets with one b— tag and £+ jets with two tags 023 The event selection follows 
Section 14.31 and 14.41 A fourth jet was not explicitly required in the £+jets sample. 
In constructing templates for the £+jets samples, the following observations can be 
made. In the single tag case, the tagged jet comes from the same top quark as the 
identified lepton in about 50% of the cases. For the double tag sample, there are 
two possible combinations of b~ tagged jet with charged lepton. This led to the use 
of a one-dimensional Ma template for the former, and a two-dimensional template 
for the latter. Templates for background were obtained from Alpgen Wbb events 
and, for the single tag CcLSG, el 15% admixture of multijet instrumental background 
from data. The dilepton analysis considers each charged lepton as an independent 
measurement. In each case, there are two possible matches to the leading two jets 
in an event, which are assumed to be the b— jets. Like the double-tagged i'+jets 
analysis, a two-dimensional template of Ma is generated for signal. The background 
was modeled by a 50%:30%:20% mix of Z — * l£, PF+jets, and diboson processes. 
The VF+jets instrumental background was obtained from data according to the 
prescriptions described in Section l4~3l The expectations from backgrounds for f+jets 
(dilepton) channels were validated in background-dominated one and two (one) jet 
samples in data. Figure [39] indicates the distribution of Ma for the £+jets data 
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sample. Both V + A and V — A expectations are shown for comparison. 
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Fig. 39. Distribution of in 700 pb 1 of CDF £+ jets events, with expected tt and background 
overlai d (le ft). Likelihood for dilepton, £+jets and combined analyses indicating estimated fv+A 



To obtain an estimate of fv+A, a binned log likelihood fit was employed to 
the data using templates from tt decays with V + A and V — A couplings and 
from background events. The tt and background cross sections were treated as 
nuisance parameters in the fit. The fitted values of these parameters were near 
the input values (e.g. — 7.3 ± 0.9 pb). The validity of the fit method was 
explored by using ensemble testing where the tt and background composition was 
allowed to vary by the estimated uncertainty in these estimates. Ensembles were 
also constructed with a varying V + A component. The measurement from the 
f+jets analysis yields the fraction of V + A current fv+A = 0.21 ± 0.28. Dilepton 
events yield fv+A — —0.64 ±0.37. These results are approximately 1.8<7 apart. The 
combined result gives fv+A — —0.06 ± 0.25(stat + sys). The likelihood vs. fv+A 
for each channel, and the combined result, are shown in Figure The primary 
systematic uncertainty came from the uncertainty in the jet energy scale. By using 
a Bayesian approach, an upper limit of fv+A < 0.29 was obtained at 95% c.l. 
These results translate into a value for the fraction of right-handed W bosons of 
/+ = -0.02 ± 0.07(stat) ± 0.04(sys) and a limit of /+ < 0.09 at 95% C.L. 

The most recent analysis of the VT-boson helicity fractions from CDFH^ was 
performed using a dataset corresponding to an integrated luminosity of 1.7fb _1 , The 
distribution of the helicity angle cosine, cos(9*), was used to obtain the fractions of 
longitudinal fo or right-handed /+ W^-bosons. Fully reconstructed £+jets tt events 
were used for this analysis. While all combinations of jet-parton assignments were 
processed via the kinematic fitter, the one with the smallest fit \ 2 was chosen as the 
final event configuration and was used for reconstructing the variable cos{9*). An 
unbinned likelihood fit to the observed cos(9*) distribution for the data was per- 
formed using templates for the longitudinal, right- and left-handed signal events and 
the background events. From this fit a measurement of the uncorrected longitudinal 
and right-handed fractions of VF-boson events in the data was obtained. To convert 
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this measurement to the true longitudinal and right-handed fractions of events, a 
correction for the selection and reconstruction efficiencies was applied. Three differ- 
ent types of fits were employed. In all three fits the constraint /_ + fo + /+ = 1 was 
imposed. Two 1-D fits correspond to the measurements of the fraction / + (fo) with 
fo (/+) fixed to the value expected from SM. If one constrains the right-handed 
fraction / + = 0, then the measurement of fo — 0.57 ± O.ll(stat) ± 0.04(syst). Con- 
straining fo = 0.7 lead to a measurement / + = —0.04 ± 0.04(stat) ± 0.03(syst) or 
/+ < 0.07 at 95% C.L.. The third fit was a simultaneous 2D fit for both parameters 
/+ and fo and lead to a measurement of fo = 0.61 ± 0.20(stat) ± 0.03(syst) and 
/+ = -0.02 ± 0.08(stat) ± 0.03(syst). 



7.1.3. D0 Results 

During Run II, D0 has performed measurements of t he W boson helicity using 
the cos{0*) distribution. In the most recent publicatio n^^ , both dilepton and £+ 
jets events are used for the measurement. The event selection in both channels is 
similar to that used for the nit measurements. The backgrounds for the f+jets events 
are mainly VF+jets production and the instrumental background from misidentificd 
multijet events. Discrimination between it events and background events is obtained 
by constructing a discriminant which is close to one for signal it events and near zero 
for background events. This discriminant is computed from variables that exploit 
the differences in the event kinematics and jet flavor of the signal and background 
events. 

The kinematic variables used for building the discriminant include Ht, 
A.(/>(£,$1t), aplanarity (A), sphericity (S) and the minimum dijet mass mjj m i n . 
These are described in Section I4TT1 The y 2 from the kinematic fit was also used. In 
addition a variable which distinguishes the flavors of the jets is obtained by using 
the impact parameters of the tracks in the jet with respect to the primary vertex. 
The impact parameters of tracks are used to form a probability for the jet to origi- 
nate from the primary vertex. The average of the two smallest probabilities are used 
to characterize the event. This variable has smaller values for top quark events than 
for background events. For the dilepton channel, backgrounds arise mainly from 
WW+jets or Z+jets processes. After the selection of the events using kinematic 
quantities, a fairly good signal to background ratio is obtained and no further cut 
on an event discriminant is required. 

Measurement of the W boson helicity is performed by fitting the cos(6*) dis- 
tributions observed in both the dilepton and £+ jets data samples. The observed 
distribution is fit to two components: the expected backgrounds and a it signal sam- 
ple which is generated at specific values of / + . During the fit, fo is fixed at 0.70 and 
/+ is varied. The cos(9*) templates for signal samples are constructed by generating 
it events with two different values of f+=0 and / + =0.30 and using a linear interpo- 
lation to generate a value of /+. The backgrounds with real leptons are generated 
using ALPGEN and pythia, and the multijet background templates are extracted 
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from the background control samples (see Section WM for details). Figure |4TH shows 
the cos{6*) distribution observed in £+ jets and dilepton events for a luminosity 
corresponding to 370 pb _1 . Overlaid on the figure are the expected distributions 
from the standard model prediction and a model with a pure V + A interaction. 
A binned Poisson likelihood £(/+) is computed for the data to be consistent with 
the sum of signal and background templates at each of the chosen /+ values in the 
range < /+ < 0.3. 

A parabola is fit to the —ln[L(f+)] points to extract the likelihood as a function 
of f + . Systematic uncertainties are folded in prior to the determination of the final 
result. The dominant contributions to the systematic uncertainties arise from jet 
energy scale measurements and the input top quark mass. They are assumed to be 
fully correlated between the dilepton and £+ jets events. With the assumption that 
f is fixed at 0.7, D0 measures /+ = 0.109 ± 0.094(stat) ± 0.063(syst) using £+ jets 
events, and / + = 0.089 ± 0.154(stat) ± 0.059(syst)using dilepton events. Combining 
the two yields /+ = 0.056±0.080(stat)±0.057(syst). A Bayesian confidence interval 
is computed using a flat prior distribution which is non-zero only in the physically 
allowed region of < /+ < 0.3 leading to /+ < 0.23 at 95% C.L. 

An earlier analysis with 230 pb _1 data was also carried out by D0 E3. The 
analysis strategies and techniques were similar to the ones described above. The 
main difference involved the further subdivision of the £+ sample into two categories 
depending on whether a tagged & — jet was found in the event. A likelihood analysis 
of the angular distribution of the leptons, leads to the result / + = 0.0±0.13(stat)± 
0.017 (syst) and /+ < 0.25 at 95% C.L. All measurements are consistent with the 
SM prediction / =0.7, /+=0. 




Fig. 40. Distribution observed in D0 ^+jets events (left) and dilepton events (right; I The 
standard model predictions are shown as the solid lines. A model with a pure V + A interaction 
would result in the distribution given by the dashed lines. 

D0 also has carried out a first model independent measurement of the W boson 
helici ty fra ctions on a larger dataset corresponding to an integrated luminosity of 
1 fb -1 !-^-. As described above, all earlier measurements of helicity fractions /+, fo 
have been performed by keeping one of them fixed to the value predicted by the 
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SM and measuring the other by a fit to one of the kinematic variables sensitive 
to the VF-boson helicity. The model independent study is based on a simultaneous 
measurement of /+, /o, with the condition that /+ + /o + /- = 1- This analysis also 
benefits from better event selection efficiency. The event selection efficiency and the 
background rejections compared to analyses on the smaller datasets are improved 
by using a likelihood based event discriminant. The likelihood discriminant is based 
on the following variables: Ht, centrality C, kximin), the sum of all jet and charged 
lepton energies h, the minimum dijet mass of the jet pairs m(jj) m i n , aplanarity A, 
sphericity S, $t and the dilepton invariant mass (for definitions of these variables 
see Sectiorj4]). The likelihood discriminants T>, are computed independently for each 
of the different channels considered in the analysis (ee, e^i, ee, e+jet, and /z+jet) 
and an appropriate selection on T> is applied to maximize the signal to background 
for each channel. In addition, the neural network 6-tagging algorithm was used to 
identify the 6-jets. 

The simultaneous measurement of /o and /+ is performed using the distribution 
of cos{9*). A further enhancement in this analysis compared to previous ones is 
the use of the hadronicaly decaying W boson in the cos(6*) distribution in the 
£ + jets events. The distributions of cos(9*) for the leptonic and hadronic W bosons 
obtained in the I + jets events is shown in Fig l41l 




cose* |cos9*| 



Fig. 41. Distribution of cos ( 8*) fo r the leptonic W bosons (left) and hadronic W bosons (right) 
obtained in the l+jets sample! The data are represented by dots, the expected signal is shown 
as the dashed histogram and the shaded histogram corresponds to the expected background. 

The simultaneous fit for /o and /+ is performed using a binned Poisson likelihood 
£(/o;/+). The consistency of the data with the sum of signal and background 
templates as a function of the helicity fractions is computed, keeping the background 
normalization constrained to the expected value within uncertainties. With this 
procedure D0 obtains 

/o = 0.425 ± 0.166(stat.) ± 0.102(syst.)and./+ = 0.II9 ± 0.090(stat.) ± 0.053(syst.). 

In Figl42l the 68%, and 95% C.L. contours from the fit are shown compared to 
the prediction from SM. While the measurement may suggest a smaller fraction 
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of longitudinal W bosons and a larger fraction of right-handed VK-bosons as com- 
pared to SM predictions, they are statistically consistent with the SM. Given the 
uncertainties in the measurement, D0 computes a 27% chance for observing the 
discrepancy. 




Fig. 42. Result of the model-independent simultaneous fit to the W boson helicity fractions fo 
and /_)_. The star denotes the SM prediction, the ellipses represe nt th e measured 68% and 95% 
C.L. contours. The triangle indicates the physically allowed region.USE] 



Measurements of fo and /_ can also be performed using the earlier prescription 
of fixing one of the them to the SM value. Fixing / + to the SM value gives fo — 
0.619 ± 0.090(stat.) ± 0.052(syst.) and fixing f to the SM value leads to a value of 
/+ = -0.002 ± 0.047(stat.) ± 0.047(syst.). 

7.2. Measurement of B(t —> Wb)/B(t — s- Wq) and \V tb \ 

The top quark decays predominantly via t — > Wq, where q can be a d, s, or a 

fo-quark. Flavor changing neutral current (FCNC) decays of the type t — >Vou where 

V = g, Z, 7 and q = u,c are of the order of 10~ 10 or smaller within the SM^SHl. In the 

SM, the ratio of top quark branching fractions, R = BR(t — > Wb)/BR(t — > Wq), 

I v h l 2 

can be expressed in terms of the CKM matrix elements, R — iv^^+ivd^+lVtdi 2 • 
Based on the assumption of three quark generations and unitarity of the CKM 
matrix, and using the very small experimentally mea sured values of \V u b\ and \ V c b\, 
one obtains 0.9990< \V tb \ <0.9992 at 90% ClM This also constrains R to be 
in the interval 0.9980-0.9984 at the 90% C.L. 48 . Thus, in the SM, the top quark 
is expected to decay nearly 100% of the time into a W boson and a &-quark. A 
significant deviation from the SM prediction of BR(t — > Wb) ~ unity could arise in 
the presence of additional quark generations or non-SM processes in the production 
or decay of the top quark. The identification of tt events at the Tevatron allows us 
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to check for such a possible deviation which would imply the appearance of new 
physics. In Run I, CDF measured R — 0.94 ^'^(stat+sys) and set a lower limit 
of R >0.61 (0-56) at 90% ( 95%) C.L., in agreement with SM predictions. In Run II, 
both CDeQ2H and D0 221] have reported direct measurements of R by examining 
the candidate ti events in data samples of ~160 pb _1 and ~230 pb _1 , respectively. 
These analyses assume that t — > Xq, where X ^ W is negligible. 

Analysis strategy: In the SM case with R w 1, the ti event signature contains 
two b— jets from the two top quark decays. However, with R of an unknown value, a 
ti event might have zero, one or two b— jets since the top quark can also decay into 
a light quark (d or s) and a W boson. The ratio R = BR(t -> Wb)/BR(t Wq) 
determines the fraction of ti events with zero, one and two b— jets and therefore how 
ti events are distributed among samples with zero, one or two-tags. Thus R can be 
extracted from the relative rates of identified b— jets in ti events. If each of the two 
top quarks in the event has a probability R to decay to a fe-quark, and there is an 
efficiency e b to tag the b— jet, the efficiencies to detect zero, one or two 6— jets in 
the event are 



e = (1 - Re b ) 2 ,e 1 = 2Re B R(l - Re b ), e 2 = {Re b ) 2 (22) 

2 1 1 
Re b = — = - — = - 23 

ei/e 2 + 2 2e / ei + l ^J7 2 + 1 V ' 

2 1 1 

Rtb = ; = ; - , (24) 

N 1 /N 2 + 2 2AV7V1 + 1 TAV^+l V ' 

where Ni are the number of ti events observed with i tags, where i can be zero, one 
or two. Any two tagging rates determine Re b uniquely, while all three tagging rates 
jointly overdetermine Re b . The ratios of tag rates can only determine the product 
Re b , as it is not possible to distinguish missing tagged jets due to tagging inefficiency 
(eb less than unity) from missing tagged jets due to the absence of 6-quarks (R less 
than unity). R can be extracted by measuring e b independently, from ti simulation 
calibrated with complimentary data samples. As R depends only on relative tag 
rates, this measurement is independent of any assumptions of the overall ti cross 
section (a^ ). 

The measurement thus follows three basic steps. After identifying samples en- 
riched in ti events, the background level in those events is estimated for different 
b— jet multiplicity. Then the expected tag rates (and, implicitly, their ratios) in ti 
events are predicted as a function of Re b . Finally, the observed tag rates are com- 
pared to the expectations, yielding the most likely value of Re b , and allo wing to set 
a lower limit on R. This strategy was pioneered by CDF's Run I analysi a^^ i 

Measurement by CDF: Based on ^160 pb _1 of Run II data, CDF followed the 
same technique and determined R from the relative numbers of ti events with dif- 
ferent multiplicity of b— tagged jets in the ^+jets and dilepton channelJ^H The 
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event selection, b— jet identification and the associated background estimation are 
essentially equivalent to the ones in the cross section analyses discussed in earlier 
sections. In particular, in the one-tag and two-tag subsamples of the £+jets sample, 
apriori estimates of the background are made with a collection of data-driven and 
simulation techniques. In the zero-tag €+jets sample, where the W+jets production 
rate dominates over the ti pairs, the event rate is determined by exploiting the 
unique kinematics of ti events with an artificial neural network (ANN) . A binned 
maximum likelihood fit of the ANN output distribution is performed for the ti frac- 
tion in the subsample. In the dilepton channel, most of the jets in the background 
events arise from generic QCD radiation. To determine the background distribu- 
tion across the subsamples with zero, one and two b-tags, a parameterization of 
the probability to tag a generic QCD jet, derived from jet-triggered data samples, 
is applied to the jets in the dilepton sample, correcting for the enriched ti content 
of the sample. Table ITT1 shows the observed number of events (N obs ) in the £+jets 
and dilepton samples, and the corresponding estimates of background levels (N bk9 ) 
and expected event yields (N exp ). The ti event tagging efficiency ej depends on the 
fiducial acceptances for jets that can potentially be b— tagged, and the efficiency to 
tag those jets. These efficiencies in turn depend on the species of the jet's progenitor 
quark. The efficiency thus depends strongly on i?, as R not equal to unity implies 
fewer b— jets available for b— tagging, and more light-quark jets instead. The jet ac- 
ceptances and tagging efficiencies are used to parameterize ei(R) using ti simulated 
events. The nominal values of ti for R equal to unity are also given in Tabic ITT1 



Table 11. Summary of observed and expected event yields as a function of b— jet multiplicity in the 
Z+jets and dilepton samples in the CDfllZl] and D0 HZZT analyses. 





O-tag 


1-tag 


2-tag 


0-tag 


1-tag 


2-tag 


CDF (Z+jets) 


CDF (dilepton) 


6i(R=l) 


0.45+0.03 


0.43+0.02 


0.12+0.02 


0.47+0.03 


0.43+0.02 


0.10+0.02 


jybkg 


62.4+9.0 


4.2+0.7 


0.2+0.1 


2.0+0.6 


0.2+0.1 


<0.01 


N l exp 
1 . 


80.4+5.2 


21.5+4.1 


5.0+1.4 


6.1+0.4 


4.0+0.2 


0.9+0.2 


Afobs 

i 


79 


23 


5 


5 


4 


2 


D0 (1+3 jets) 


D0 (l+>4 Jets) 


N tt 


32.4+1.6 


32.3+1.6 


8.2+0.5 


35.6+2.8 


41.5+3.3 


13.5+1.4 


N exp 
i 


1275+44 


79+5 


11.4+0.8 


297+19 


56+4 


14.4+1.4 


Kfobs 
i 


1277 


79 


9 


291 


62 


14 



The expected event yield in each of the three tagged subsets of ^+jets and 
dilepton samples is 

Nr = N? n MR) + N» k9 (25) 
where A/|* c = J2i{N° bs — N bk9 ) is an estimate of the inclusive number of ti events 
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in the sample. The analysis constructs a likelihood function using this information 
to get the best estimate of Rtb- The full likelihood is a product of independent 
likelihoods for the £+jets and dilepton samples. Each likelihood function is built 
as the product of the Poisson functions comparing N° bs to N^ xp for each value of 
i, multiplied by Gaussian functions which incorporate systematic uncertainties in 
the event-tagging efficiencies and backgrounds, taking into account the correlations 
across different sub-samples. The resulting likelihood as a function of R is shown in 
Fig. [43l along with the negative logarithm of the likelihood, which gives a central 
value of R — lA2^'ll(stat)^\l(sys) most consistent with the observation. 

Measurement by D0 : The D0 experiment measures simultaneously the ratio 
of branching fractions R together with the tt production cross section, using 230 
pb _1 of data and sele cting a sample enriched in tt events which are consistent 
with f+jcts final statefl^l The tt enriched sample is divided into I + 3jets and 
£+ > 4jets, which are further categorized on the basis of zero, one and two or more 
b— jets. The analysis fits simultaneously R and the total number of tt events (N^) 
in the zero, one and two-tag samples to the number of observed one-tag and two-tag 
samples, and, in the zero-tag events, to the shape of a discriminant variable (D) that 
exploits kinematic differences between the background and the tt signal. The event 
selection, background determination and secondary vertex b— tagging algorithm are 
essentially equivalent to that of the corresponding cross section analyses described 
in the section 4.4. 

In an analysis based on the SM, with R unity, the event tagging probabilities 
are computed assuming that each of the signal events contain two b— jets. In the 
present analysis with R not equal to unity, the event tagging probability becomes a 
function of R since a tt event might have zero, one or two b— jets. To derive the event 
tagging probability as a function of R, the event tagging probability is determined 
for the three following scenarios (i) tt — » W + bW~b, (ii) tt — > W + bW~q~i, and (iii) 
tt — > W + qiW~qi, where qi denotes either a d— or s— quark. The probabilities 
Pi-tag to observe i — tag — zero, one or > two b— jets are computed separately for 
the three types of tt events, using the probabilities for each type of jet (b, c or light 
flavor jet) to be b— tagged. The probabilities Pi-tags in the three scenarios are then 
combined in the following way to obtain the tt tagging probability as a function of 
R: 

Pi-tag(ti ) = R 2 Pi-tag(i) + 2JJ(l - R)Pi-tag{H) + (1 - Rf Pi-tag (Hi) (26) 

where the subscript i — tag runs over zero, one and at least two tags. Table [IT] 
compares the observed number of events with the sum of the predicted backgrounds 
and the fitted number of tt events as a function of the number of 6— tags. To measure 
R and , a binned maximum likelihood fit is performed to the data. The values of 
R and o>r that maximize the total likelihood function are R = 1.03^g '^(stat+sys) 
and — 7.9^'5(stat+sys)±0.5 (lum) pb, respectively, and in good agreement 
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with the SM expectation. The result of the 2-dimensional fit is shown in Fig. 03] 
in the plane (R, <r^ ), along with the 68% and 95% confidence level contours. A 
preliminary simultaneous measurement of and R was performed in 900 pb _1 of 
data^^2l. This gives a measured value of R = 0.97^Q;gg(stat + sys). 




Fig. 43. Left: The upper plot shows the likelihood as a function of R (inset) and its negative log- 
arithm. The intersections of the horizontal line ln(L)=-0.5 with the likelihood define the statistical 
1it errors on R. The lower plot s hows 95% (outer), 90% (central), and 68% (inner) C.L. bands for 
Rtrue as a function of R by CD F [ Right: The 68% and 95% statistical confidence contours in 
the (R,N^ ) plane by D0 LLZ2J The point indicates the best fit to the data. 

With their measured values of R, CDF and D0 have set a lower limit on 
R. Assuming that non-VF decays of the top quark can be neglected, that only 
three-generations of fermions exist, and that the CKM matrix is unitary, implying 
R = \Vtb\ 2 , they have extracted the CKM matrix element \Vtb\- The results of this 
measurement are summarized in the Table 1121 All the measurements of R and | Vtb \ 
are consistent with SM expectations. These results are presently limited by statistics 
and will profit from the upcoming larger data sets. 

Table 12. Published measurements and lo wer 

limits of R = BRH—> Wb)/BR(t -> Wq) and |V tb | from CDfCE] and 
D0 Q22] in R un n. 



CDF (/ £dt=160 pb- 1 ) D0 ( / £dt=230 pb" 1 ) 



R 


i 19 +0.27 
± - lz -0.23 


1 n3+ 019 


R (95% C.L.) 


>0.61 


>0.61 


R (68% C.L.) 




>0.78 


\V tb \ (95% C.L.) 


>0.78 


>0.78 


\V tb \ (68% C.L.) 




>0.88 
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First direct determination of \Vtb\ by D0 : Single top quark cross section 
measurement provides a direct determination of the CKM matrix element \Vtb\ 
since the cross section is proportional to the square of this quantity (a <x \V t b\ 2 )- A 
value inconsistent with the SM expectation \Vtb\ —1 would be a signature for new 
physics such as a fourth quark family. The D0 experiment has performed the first 
direct measurement of \ Vtb\ based on the single top quark cross section measurement 
derive d from decision trees analysis (discussed in section 5.2) using about 0.9 fb _1 
of datcM This result makes no assumptions on the unitarity of CKM matrix or 
the number of families, but it does require a few assumptions. First, it is assumed 
that the single top quark production meachanism only involves interaction with a 
W boson, and not from one of the various beyond SM scenarios that include extra 
scalar and vector bosons or FCNC interactions. The second assumption made is 
that \V t d\ 2 + \V ts \ 2 <§C \V t b\ 2 , which is experimentally supported by the BR(t — > 
Wb)/ BR(t — > Wq) measurements done on ti events. This requirement implies that 
BR(t — > Wb) ~100% and that single top quark production is completely dominated 
by the tbW production. Finally, it is assumed that the tbW interaction is CP- 
conserving and of the V — A type, but it is allowed to have an anomalous strength. 
Two measurements have been performed, one for the strength of the V — A coupling 
\Vtbfi\ in the Wtb vertex, where /f is an arbitrary left-handed form factor, with 
no requirement that it be less than one; and one assuming /f =1 resulting in |Vt&| 
being restricted to between and 1. The limits calcualte d usi ng Baysian approach 
are: \V tb f[\ =1.3±0.2 and 0.68< \V tb \ <1 at the 95% C.lPI. 



7.3. Top Quark Charge 

Many of the currently measured properties of the top quark are still only poorly 
known, and the indirect constraints set by precision electroweak data leave plenty 
of room for new physics. In particular, its electric charge (qt), one of the most fun- 
damental quantities characterizing a particle, has not been directly measured so far. 
It still remains not only to confirm that the discovered quark has a charge +| as 
assigned by the SM, but also to measure the strength of its electromagnetic (EM) 
coupling to rule out anomalous contributions to its EM interactions. In the top 
quark analyses of the CDF and D0 experiments, the correlations of the 6-quarks 
and the W bosons in pp — > ti — > W + W~bb are not uniquely determined which re- 
sults in a two-fold ambiguity in the pairing of W bosons and 6-quarks, and, hence, 
in the electric charge assignment of the top quark. In addition to the SM assign- 
ment t — > W + b, it is conceivable that the 't-quark' is an exotic quark Q with charge 
q= — | which decays via Q — > W~b. This alternative interpretation is consistent 
with current precision electroweak data. It is possible to fit Z — > £ + £~ and Z — » bb 
data assuming a top quark mass of m t =270 GeV, provided that the right-handed 
6-quark mixes with the isospin +| component of an exotic doublet of charge — | and 
— | quarks, [Q\ , Qa)^^- In such a scenario, the particle we have been exploring 
may in fact be the — | charge quark and the top quark, with mt=270 GeV, would 
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have so far escaped detection at the Tevatron. 

Specific qt measurements are therefore required to rule out or confirm one of 
the hypotheses. In order to determine gt, one can study the charge of its d ecay 
products or investigate photon radiation in ti events. The latter methocP^] con- 
sists in cross-section measurements of ti pairs where a photon is radiated, either 
at the production (pp — > if 7) or the decay (pp — > ti , t — > Wiry); the radiation 
being dependent on the value of qt- Measurement of qt using radiative if processes 
is hopeless at the TevatrorlA^I - even with 20 fb^ 1 of data, a charge — | top quark 
can be excluded at «95% confidence level. However, at the LHC, with 10 fb _1 , one 
can hope to determine qt with a precision of ~10%. The first method - measuring q t 
by reconstructing the charges of its decay products (the final-state 6— jets and the 
W bosons) for semi-leptonic decays - seems feasible at the Tevatron. The W boson 
charge in leptonic decays is accurately given by the lepton coming from it and the 
&-jet charge can be determined from a measurement of charges associated with the 
tracks in the jet. It should be noted that since both t and i quarks are present in 
every event, the analysis is only sensitive to the absolute value of qt- 



The D0 experiment has performed the first determinatio n I of qt with 
^370 pb _1 of Run II data on a double-tagged semi-leptonic (pp — > ti — > l^vjj'bb) 
sample of ti candidate events. The first step in the measurement involves selecting 
a pure sample of ti events in data in the £+jets channel. Events with an isolated 
high pt (> 20 GeV) e or /i accompanied by four or more high pt (> 15 GeV) jets 
are considered for the analysis. Events with W bosons are selected by requiring $t 
> 20 GeV. To remove the multijet background, J£t is required to be acollinear with 
the lepton direction in the transverse plane. The purity of the sample is significantly 
enhanced by requiring at least two jets in the event to have a secondary vertex b- 
tag. In the resulting sample of 21 candidate events, Wbb production is the largest 
background and represents ~5% of the sample. 

Each selected ti event has two "legs", one with a leptonically decaying Wi 
(t — ► Wb — > Ivbi) and one with a hadronically decaying Wh (t — > Wb — > qq'bh)- 
The second step in the analysis consists in assigning the correct jets and leptons 
to the correct "leg" of the event, uniquely specifying which &-jet comes from the 
same top (or anti-top) quark as the lepton. A constrained kinematic fit is performed 
for this purpose. The four highest pr jets can be assigned to the set of final state 
quarks according to many permutations and there are at least two ways to assign 
the b— jets to b\ and bh- For each permutation, the measured four vectors of the jets 
and leptons are fitted to the ti event hypothesis and the bi and bh jets are identified 
by selecting the permutation with the highest probability of arising from a ti event. 
For 16 of the 21 selected events, the kinematic fit converges leading to correct as- 
signment, thus providing 32 measurements of qt since top quark charge is measured 
twice in each event. In each ti event, the observable \qt \ is computed, which is the 
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sum of the lepton charge from the W^-boson decay and the charge of 6/. The 6— jet 
charge is determined by taking a py-weighted sum of the charges of the tracks with 
Pt > 0.5 GeV contained within a cone of AR =0.5 around the b— jet's axis and 
thus one can distinguish between the b and 6 on a statistical basis. The performance 
of this jet charge algorithm is calibrated using an independent dijet data sample 
enriched with semileptonic b quark decays. In the calibration sample, the charge of 
the 6— jets is derived from the lepton charge, accounting for c— jet contamination 
and B oscillation corrections. The performance of the algorithm is also calculated 
as a function of the pr and rj of the jet, so that kinematic differences between the 
control sample and the ti are accounted for. 

The subsequent step is to use the expected shapes of jet charge for 6-jet and 
&-jet in data to derive the expected shape of |g| for the SM (q t = +2/3) and the 
exotic (q t = —4/3) scenario. The resulting SM and exotic templates, normalized 
to unity, are used as probability density functions, p sm and p cx , respectively, in 
the C.L. calculation. Figure. [44] compares the charge distribution of the top quark 
candidates reconstructed in data with +2/3 and -4/3 charge hypotheses. 




Top quark charge [e] Fraction of exotic quarks 



Fig. 44. The 32 measured values of the top quark charge (binned) compared to the Standard 
Model and exotic scenario templa tes (l eft )=-^U Likelihood fit of the fraction of exotic quark pairs 
in the selected data sample (right) 1 ^®! 



A likelihood ratio defined as 

is then computed to determine the most probable model where the numerator (de- 
nominator) measures the likelihood of the observed set of charges qi arising from a 
SM top quark (exotic quark). The subscript i runs over all 32 available measure- 
ments. The value of the A data measured in data is compared with the expected A sm 
and A ex distributions derived by performing ensemble tests using the SM and exotic 
scenarios respectively. The observed set of charges is found to agree well with those 
of a SM top quark. For the exotic heavy quark hypothesis, only 7.8% of the pseudo- 
experiments give a higher A ratio than the one measured in data. Therefore, D0 
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yields a p— value, corresponding to the probability of consistency with the exotic 
model, of 7.8%. The top quark is indeed consistent with being the SM \q t \ =2/3 
quark. Figure 02] shows the fraction of exotic quark pairs (p) in the data obtained 
by performing an unbinned maximum likelihood fit to the observed set of which 
yields < p <0.52 at the 68% C.L and < p <0.80 at the 90% C.L. 

8. Searches for Non-Standard Top Quark Production and Decay 

Several alternative models predict substantially different production and decay 
mechanisms for the top quark than were described in Section O For instance, the 
possibility to detect a heavy resonance which decays to tt is testable using elements 
of the techniques for kinematic reconstruction. The final states are sensitive to 
alternative decay chains from the top quark, such as thru a charged Higgs or W 
boson. This section will review the experimental efforts in these directions. 
8.1. tt resonance search 

It may be that the top quark is a special member of the fermion family. Its large 

mass has provoked suspicion that it may have a unique and significant role to play 

11 77I44J 

in electroweak symmetry breaking. Topcolor model d 11 suggest that there are 
new strong interactions which can generate the large top quark mass while leaving 
other quarks light. This interaction results in a massive tt condensate, sometimes 
termed Z . Technicolo r on the other hand, posits interactions at the electroweak 
scale which provide electroweak symmetry breaking without the need for a scalar 
Higgs field. In the extended models, this gives masses to all fermions. However, the 
generated value of m t is too small. The top-color assisted technicolo J ^ 77 | 44 ] moc jel 
combines these approaches to predict the correct m t . Variants of this approach 
predict a Z' which couples preferentially to third generation fermions and has no 
couplings to leptons. 
8.1.1. General Methods 

Efforts to look for a tt resonance have employed the €+jets final state. In such events, 
a full kinematic reconstruction is performed as in the mass analyses described in 
Section 16.41 with the exception that m t is either held fixed or within a window 
determined by uncertainty on the mass. All jet combinations and neutrino solutions 
are tried. The chosen combination is one that minimizes the x 2 , given the known 
errors and measured parameters of the event. The combinatorics of jet assignment 
mean that this choice is not always the correct one, although b— tagging can enhance 
the probability. Generally, a cut on the x 2 is applied to further reduce backgrounds. 
The solution for each event permits the calculation of the apparent tt invariant 
mass, M t (, from that event. 

A narrow resonance with Tz> = 0.012Mz' is the standard assumption for the 
signal. Generally, this width is substantially less than the resolution inherent in 
the experimental measurements, so natural widths up to a few percent may still 
be accomodated. However, this assumption does mean that any results are not 
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strictly applicable to some models, such as Kaluza-Klein scenarios. The primary 
backgrounds to the search are standard tt production, VF+jets and QCD multijet 
production. Single top quark and diboson events give small but non-negligible con- 
tributions. Evidence of resonant production is sought by fitting the expected 
resonance plus the expected tt and backgrounds to the observed spectrum in data. 
A limit on any anomalous production cross section results on a limit on the mass 
of the resonant state. 



8.1.2. D0 Searches 

In 2004, D0 published a search for tt resonance in Run I dat A sample of 41 
£+jets events, 4 with soft il— tag selection, were used for the analysis. The resonance 
signal was modeled with Pythia in a mass range of 400 GeV < Mz' < 850 GeV. 
The CTEQ3M pdf was used. The tt background was simulated using Herwig. The 
Vy+jets background was modeled with Vecbos, and Herwig was used for the par- 
ton shower simulation. The instrumental background from multijet production was 
derived from a data sample in which actual leptons h ave b een rejected. No evidence 
for resonant production was observed. A Bayesian fit HSH W as performed which con- 
sidered three components separately: signal (Z' — > tt), standard tt production, and 
W + jets and multijet events. The latter element had W and QCD fractions set 
from the template mass analysis. The resulting limit is Mz 1 > 560 GeV/c 2 at 95% 
C.L.. 
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Fig. 45. Distribution of tt invariant mass, M^, from D0 Run II data (left) 
observed limits, with comparison to the Z' estimation is also shown (right). 



Expected and 



Two preliminary studies have been executed in Run II using template meth- 
ods for kinematic reconstruction. A mass range 350 GeV/c 2 < Mz' <900 GeV/c 2 
was considered in both analyses. The first search used 370 pb" 1 of £+jets events 
tagged with a secondary vertex tagger. Signal Monte Carlo samples were generated 
using Pythia. Studies indicated that the choice of jet-to-parton configuration was 
correct about 65% of the time after the x 2 cu t- The tt and Vy+jets backgrounds 
were generated using Alpgen plus Pythia. The QCD background was estimated 
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via matrix method in the data. The data yielded 108 events which showed no evi- 
dence of a resonance. A Bayesian approach was used to extract a limit. Kinematic 
shape uncertainties arose from jet energy calibration and efficiency and the knowl- 
edge of the kinematic dependence of the 6— tag and mis-tag rates. Normalization 
uncertainties came from theoretical uncertainty in the ti cross section as well as 
the integrated luminosity. Correlations among uncertainties were included in the 
fit. A limit of M z > < 680 GeV/c 2 was obtained at 95% C.L. OH. A further study 
has been performed with 197 £+jets events obtained from 900 pb _1 of data. These 
events were tagged with a neural network tagger. In this sample, 154 events were 
expected from standard ti process. Much of this ana lysis is similar to the 370 pb _1 
analysis. A limit of Mz> > 680 GeV/c 2 was obtained HHS1 Results are shown in Fig. 

su 

8.1.3. CDF Searches 

An initial searc h for the ti resonance was conducted by the CDF collaboration us- 
ing Run I data One of the leading three jets in each event were required to 
be b— tagged with a secondary vertex. Vecbos was used to model W+jets events. 
Pythia was used to model the signal as Z' — ► ti. In order to improve the recon- 
structed mass resolution of the ti system, the requirement on the computed top 
quark masses was relaxed to fall in the range 150 GeV to 200 GeV. This allowed 
a more efficient acceptance of combinations with correct jet-parton assignments. 
Twenty events were removed out of the initial sample of 83 because they yielded a 
poor minimum \ 2 - No evidence of a resonance was observed in the data. A binned 
likelihood fit of the data to signal and background was performed to extract a limit 
on the resonance production rate. Masses from 400 GeV to one TeV were scanned. 
The primary sources of systematic uncertainty were from signal and background 
shapes, and jet-related systematics such as energy scale and gluon radiation. Un- 
certainties were incorporated into the limit calculation as Gaussian uncertainties. 
The 95% C.L. upper limit for T z > = 0.012M Z / is M z > > 480 GeV. 

CDF has improved upon this result by looking at 680 pb _1 of Run II data^^l. 
The signal and standard ti production were modeled with Pythia. Multijet back- 
ground came from data and VF+jets events were modeled with Alpgen plus Her- 
WIG. The top quark and boson pair processes were normalized based on the theoreti- 
cal cross sections and the estimated luminosity to be 199 and 14 events, respectively. 
The QCD and W+jets backgrounds were then taken as the remaining contribution 
to the event yield observed in data: 450 events. The QCD/W ratio was fixed at 
10% as taken from the cross section measurement ^1 Event reconstruction was 
done via matrix element fit using resolution functions for the jets. All possible jet 
configurations were summed over in this approach, and a probability distribution 
vs. m,j.T was obtained. Simulations showed that the mean of this distribution was 
best correlated with the actual resonance mass. The mass resolution was set by the 
jet energy calibration and the unknown p z of the neutrino. Incorrect jet assignments 



96 Kehoe, Narain, Kumar 





Fig. 46. Distribution of ft invariant mass, M^jr, from CDF Run II data (left) CSS. Expected and 
observed limits, with comparison to the Z' estimation is also shown (right). 

also produced a significant low mass tail. No evidence for an excess over standard 
expectations was observed. The presence of a resonance would be detected by a like- 
lihood calculation based on the distributions < > for signal of varying masses 
and the backgrounds. Systematic uncertainties were extracted by varying the jet 
energy scale, initial and final state radiation, and the Q 2 for iy+jets production by 
±1<7 and generating these mass distributions for use in the likelihood fit. A limit of 
M z > < 725 GeV/c 2 at 95% C.L. was extracted from the data (see Fig. |46j). More 
recently, CDF has analyzed 955 pb _1 of data using a template fitting method. A 
limit of M z < < 720 GeV/c 2 was obtained QUI. 

8.2. Charged Higgs decays 
8.2.1. Two Higgs Doublet Models 

In the electroweak model, the generation of mass is governed by one scalar Higgs 
doublet manifesting in one observable particle, the Higgs boson. The peculiar prob- 
lems associated with the Higgs mechanism have led to several efforts to produce 
models with extended Higgs sectors. The simplest configuration has two Higgs dou- 
blets and this is included in several scenarios, including supersymmetry (SUSY). 
Two types of model exist. Type I models have only one doublet coupling to fermions. 
Type II models, however, posit one doublet which couples to up-type fermions 
including charged lept qns, and another doublet that couples only to down- type 
fermions. In the MSSlvlESl, the Type II case holds. Five physical particles are pre- 
dicted: two neutral scalars (ft, , H°), a neutral pseudoscalar (A ), and two charged 
scalars (H ± ). 

Evidence of the existence of the charged Higgs has been sought in e + e~ 
collision d^®^ . Indirect searches from CLEO provide the constraint m# > [244 + 
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63/(tan/3) 1/3 ] GeV at the 95% C.L. 186 , where tan/3 is the ratio of vacuum ex- 
pectation values of the two Higgs doublets. Observed r — > v T K and K — ► vjjj) 
branching ratios provide a limit of 0.21 GeV -1 > tan/3/m ff at 90 % C.Lp2Zl. 
These limits are stricter than those quoted for the Tevatron searches below, but 
are more theory dependent. The direct search for H ± at LEP, however, provides a 
hard limit of run > 78.6 The Tevatron searches consider potential physics 

beyond this limit. 



8.2.2. Production and Decay 

Single charged Higgs production would have negligible cross section at the Tevatron. 
Weak pair production of charged Higgses is calculated to have a cross section of 0.1 
pb. However, a large enhancement occurs when they are generated through the 
production and decay of top quark pairs. Then the Higgs pair production rate is 
potentially as high as the 6.7 pb tt cross section. Additionally, the top quark decay 
produces associated b— jets which allows a further background reduction relative to 
direct diboson production. 




Fig. 47. At left is shown the dependence of BR(H —* TV) and BR(t — > Hb) on tan 0. Coverage in 
the mjj vs tan/3 plane is shown at right when BR(t — > Hb) > 0.5 (cross-hatched) and 0.9 (solid). 
Regions where specific final states dominate are indicated. 



If the H^ are lighter than the top quark, then the decay t — > Hb can com- 
pete with the standard t — > Wb decay. In particular, BR(t — > Hb) = m 2 cot 2 /? + 
rn 2 tan 2 (3 + 4m 2 rn 2 is significant when tan/3 is either very large or very small (see 
Fig. [47]). Figure [47|also indicates that this branching ratio declines as mn increases. 
There are then three top quark decay modes to consider: tt — ► HHbb, HWbb, and 
WWbb. The final states are dictated by the H and W decays, and the analyses all 
assume fermionic decays for the charged Higgs. 

In the Type II models, H ± decays are expected to be quite different than the 
normal W boson decays. Its decay branching ratios are dependent on m t , M H and 
tan/3. For values of tan/3 > 2, the H — > tv decay dominates and grows to 100% 
when tan/3 > 5, as shown in Fig. 1471 The r is identified explicitly only through 
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its hadronic decays to a narrow jet, as described in Section 13.21 For smaller tan (3 
values, the H — > cs decay is primary and leads to two quark jets in the final state. 
However, when Mh > 130 GeV the decay chain H — > tb — > Wbb exceeds this 
mode. The dominance of the H — > tv or cs modes have led to two different search 
strategies for low and high tan/3 regions. 

8.2.3. High tan/3 Searches 

The more direct searches attempt to identify the presence of the H ± explicitly from 
the presence of the r in a selected top quark event sample. These searches have 
the benefit that they can avoid some dependencies on theoretical calculations, since 
they measure the actual rate of r production. In general, an assumptio n is made 
that BR(H — > tv) ~ 100%. CDF conducted such a search in Run I data using 
|5t triggers. Two event topologies were explored for this analysis. One searched for 
events possessing a r + /(= e, fi, r) + $t+ two jets where one of the jets had to 
be tagged as coming from a b— quark. The other topology exhibited two r jets that 
were not azimuthally back-to-back. It was employed to account for the case where 
Mh approached m t and the b— jets were too soft to be reconstructed. Instrumental 
backgrounds from fake r jets are dominant. They are modeled from the data by 
folding a jet — ► r fake rate into the multijet +J£t sample. W and Z boson pro- 
duction was modeled using Vecbos and normalized to the measured cross sections 
in data. Standard ti backgrounds were also determined to be significant. The total 
background estimation was 7.4 ± 2.0 events and seven events were observed. To 
extract a limit on H production, the polarization of the r's should be considered 
because it affects the v angles and produces larger $t- Isajet was modified to 
account for this and expected signal efficiencies were extracted. For high tan /3 (i.e. 
> 100), the region M H < 147 GeV was excluded at 95% C.L.. 

Theoretical progress led to the realization that higher orde r radiativ e corrections 
impact branching ratio calculations in the high tan/3 regionQ^HIIHI] This means 
that limits in the tan/3 vs. Mh plane depend on m odel parameters. CDF incorpo- 

II QOI 

rated this into their analysis of I + r events in Run I — I Both tracking-centric and 
calorimeter-centric r identification schemes were used. The explicit identification of 
jets by b— tagging was not implemented in this analysis. The primary backgrounds 
were Z and W boson production in association with jets, ti also provided a sig- 
nificant contribution. This was estimated by anchoring the expected event yield to 
CDF's measurement of Uxj = 5.1 pb in the £+jets channel. The total background 
was estimated to be 3.1 ±0.5 events, and four events were observed. By considering 
the observed data and the estimated backgrounds and their uncertainties, a limit 
was determined for BR(t — > Hb). Acceptance for WW,WH and HH modes was 
calculated using Pythia with m t — 175 GeV. Tauola was used to provide the 
correct polarization of the r. Higgs masses from 60/c 2 to 160 GeV/c 2 were scanned, 
and a limit of BR(t -> Hb) < 0.5 or 0.6 in the range 60 GeV < M H < 160 GeV 
was obtained. In the MSSM at very high tan/3, this limit is not valid because the 
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relevant Yukawa coupling is non-perturbative. 

A direct search was also pursued in 2002 using 62.2±3.1 ptr 1 by DefESH. Events 
comprised of r + $t + jets and possessing a spherical distribu tion to the object 
Et's were selected. A neural network (NN) based on jetnetESI incorporated 
the $t and two eigenvalues of the normalized momentum tensor to discriminate 
signal and background. The NN was trained on signal generated with Isajet with 
forced decays of H — > tv and r — > hadrons. The Higgs mass was taken as 95 
GeV since it only weakly affected the NN performance. The top quark production 
cross section was taken as 5.5 pb. The instrumental background was modeled using 
multijet events from data. W+jets events were modeled using Vecbos plus Isajet 
for soft QCD evolution. A total of 5.2 ± 1.6 events were expected from background 
and standard top quark processes, while three events were observed in data. The 
probability for the number of expected events to fluctuate to the number observed 
was calculated for individual choices of tan and m# . Systematic uncertainties for 
jet energy scale uncertainty, signal modeling and t identification were the main 
sources considered. Values of tan/3 > 32.0 were excluded at 95% C.L. for mn = 75 
GeV. The result is shown in Fig. 051 

8.2.4. Low tan/3 Region and General Searches 

As already noted, the unique signature from Higgs — > r jets is confined to the high 
tan [3 region. However, at low tan (5, it is also possible for substantial numbers of 
top quarks to decay to charged Higgses. The Higgs decays to cs in this region are 
less distinctive than in the previous cases. A search incorporating the low tan /3 
region was first performed on 109.2 pb -1 by D0 in R un I. They employed a tech- 
nique to look for the disappearance of SM top quarks^^ and used the selections 
for the standard £+jets channels, both topological and soft fi tagged, in the ti pro- 
duction analysis^!. This strategy has the advantage of incorporating the cs decays 
into the search. However, at high Mji > 130 and tan/3 < 2, the Higgs boson will 
preferentially decay via a virtual top quark to Wbb. This gives a final state that 
is not different enough from standard top quarks, so the search is insensitive in 
this region. Since the presence of non-standard top quark decays would affect the 
measurement of the top quark production cross section, the theoretical value of 
5.5 pb at 1.8 TeV was used and assumed to be unaffected by the presence of the 
new Higgs sector. Also, studies indicated that the non-standard decays would not 
affect the existing measured top quark mass by more than five percent, so m t = 175 
GeV was used for the analysis. Isajet was used to calculate signal efficiencies and 
was modified to properly account for Wbb final states. Pythia was also modified 
to provide a cross-check. Because the selection was the same as used for the cross 
section measurements, those existing background estimates were valid. The esti- 
mated background plus t — > Wb signatures provide 30.9 ± 4.0 events compared to 
30 events observed in data. Efficiencies were calculated for different mjj's from the 
Monte Carlo to extract the sensitivity at different points in the (m#,tan/3) plane. 
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Fig. 48. D0 Run I limits jS^j in the Mu vs. tan /3 plane incorporating direct and indirect searches 
(top). CDF Run II limits USD using radiative corrections are shown at bottom. In the absence of 
these corrections, the CDF results would be the most stringent experimental limits. 



Leading order theoretical calculations of branching ratios were used to calculate the 
probability to observe the charged Higgs decays. Primary systematic uncertainties 
were jet energy scale, signal model and particle identification. At 95% C.L., most 
of the points giving a BR(t — > Hb) > 0.45 were excluded. Combining this analysis 
with the direct limit from D0 gives a branching ratio limit of BR(t — > Hb) < 0.36 
at 95% C.L. for 0.3 < tan/? < 150 and m H < 160 GeV (see Fig. 08). 

Since this analysis, it was realized that i? — > Wh° can be significant. CDF has 
incorporated this element, along with the developments described for t heir earlier 
analyses, into the first search for charged Higgs bosons in Run iJHH] Using 193 
pb _1 , the yield of events in four categories of final state were compared: dilepton, 
£+jets with b— tag, £+jets with at least two b— tags, and £ + t events. The latter was 
compared to the other channels to produce direct limits on the presence of charged 
Higgs bosons in top quark decay. This analysis was the first to use the newer Yukawa 
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Table 13. Limits by D0 CSS an d CDfES] in high tan region 
for charged Higgs branching fraction. 



Experiment 


Lumi. (pb^ 1 ) BR(t - 


■+ Hb) < @ 95%C.L. 


D0 


66.2 and 109.2 


0.36 


CDF 


193 


0.4 



coupling and radiative cor in the determination of acceptances for 

the sensitivity calculation. The analysis considered decays of the Higgs to rv, cs, t*b, 
and Wh° — > Wbb. Pythia was modified to include the t*b decay mode and used to 
estimate the efficiencies for various top quark widths, the charged Higgs mass and 
width, and m^o. In order to properly calcu late the overall acceptance for various 
parameter choices in MSSM, the CPSuperlfES 

generator was used. Much of the 
low tan [3 region was excluded, including some of the high m# region. In the high 
tan/3 region, a limit of BR(t — ► Hb) < 0.4 at 95% was obtained, as shown in Fig. 
l48l If no assumptions are made for charged Higgs decay (i.e. H — * tv), the limit 
is relaxed to < 0.91. For comparison (see Table IT3)) . these limits would be stricter 
than the previous D0 limits if the radiative corrections were ignored. 



8.3. W' in Top quark decays 

The existence of new forces in nature can be revealed through the observation of 
additional gauge bosons beyond those of the s tandard model. Various extensions 
of the SM postulate larger gauge groups ^^^ and therefore new forces associated 
with additional charged gauge bosons, which are generically called W. For instance, 
the left-right symmetric modeP^ expands the SU(2)l x U(l) electroweak group 
to SU(2)l x SU{2)r x U(l), predicting the existence of three new gauge bosons: 
two charged W bosons and one neutral Z' boson. In this model, the W boson 
appears as a heavier counterpart of the left-handed W and is responsible for right- 
handed interactions, in the same way as the SM W boson mediates only left-handed 
interactions. 

Previous indirect searches based on low energy phenomena such as fj, decay, the 
Kl — Kg mass difference, neutrinoless atomic double beta decay, and semileptonic 
branching ratio b — » Xlv have resulted in stringent model-dependent limitA^I. 
Direct searches for the production of W bosons have focused on its leptonic decays, 
W' — > lv\, by looking for anomalous prod uction o f high trans verse mass lv\ pairs. 
Searches using the decay mode W — > e^J^QSEQH (W' — > fiU/ p^^ ) from Tevatron 
Run I exclude a W boson with mass m w > < 754 GeV (660 GeV) at 95% C.L.. 
The combination of th ese leptonic channels yields the lower limit on m\y> of 786 
GeV at 95% ClP. These mass limits all assume that the new vector boson's 
couplings to leptons are those given by the SM, with the additional assumption 
that the neutrino produced is much lighter than the W boson. Searches for the 
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W boson in the quark decay channel W — > qq' based on resonant structure in the 
dijet mass spectrum have also placed strong constraints on its mass. Searching in 
the hadronic final state avoids assumptions regarding the neutrino mass m Vl but 
is background limited. Previous direct searches in t his ch annel have ruled out W' 
bosons in the mass range 1< m W ' <261 GeV by UA^^S], 300< m w > <420 GeV by 

CDF EM! and 300< mwi <800 GeV by D 0EQ5] at 95% C L 



8.3.1. W in Top quark sector 

The top quark sector offers great potential for the search of new gauge bosons. The 
single top quark final state is especially sensitive to the presence of an additional 
heavy boson, owing to the decay chain W — + tb, where the top quark decays to 
a 6-quark and a W boson which subsequently decays leptonically or hadronically. 
The leading order Fcynman diagram for W' boson production resulting in single 
top quark events is shown in Fig. l49l 




Fig. 49. Leading order Feynman diagram for single top quark production via a heavy W' boson. 
The top quark decays to a standard model W boson and a fe-quark. 



This diagram is identical to that of the SM s-channel single top quark production 
where the SM W boson appears as a virtual particle. The W boson also has a 
t— channel exchange that leads to the single top quark final state. However, the 
cross section for a t— channel W process is much smaller than the SM t— channel 
single top quark production due to the high mass of the W boson. A heavy W 
boson signal would appear as a peak in the tb invariant mass distribution in these 
processes. Although this search is only sensitive to W bosons with mass above the 
tb kinematic threshold of approximately 180 GeV, it is relatively free of background 
compared to the W' — > qq' decay mode because of the distinct signature from 
the top quark decay t — > Wb. Furthermore, the interpretation of the data is less 
sensitive to assumptions regarding the right-handed neutrino (vr) sector or the 
lepton couplings of the W boson. 
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8.3.2. Run I search for W -> tb 

CDF attempted the first search for a W boson decaying to a tb final state using 
a 106 pb _1 of Run I data sample^^. The analysis considers the ^+jets final state 
topology of single top quark production arising from leptonic decays of W boson, 
W — ► eu, \iv. The event signature consists of a high-pr lepton, significant $t from 
the neutrino, and two 6-quark jets. The candidate events are selected by requiring 
a high-pr electron or muon with large $t and accompanied by two or three jets, 
at least one of them being 6-tagged. The fully simulated Pythia Monte Carlo 
was used to determine the expected contribution from signal events in the data 
sample as a function of W mass. The W boson is required to have a right-handed 
coupling to the tb final state since negligible signal yield differences between right- 
handed and left-handed couplings are expected. The analysis considered both the 
cases when mw < m VR and mw m VR , where m VR is the mass of vr that 
couples to the W . The analysis attempted to interpret the data for mw >225 
GeV since the acceptance calculation becomes increasingly uncertain as one nears 
the tb kinematic threshold. The dominant background contributions to this search 
arises from the pair production of top quarks, single top quark production, and the 
associated production of W bosons with one or more heavy quarks (Wbb, Wcc) 
which provides the largest single background contribution. A candidate sample of 
57 events agreed reasonably with the 48±6 expected from background. Therefore, 
no significant evidence for W' boson production is seen. 

In such a scenario, the analysis sets a limit on the W production cross section 
times branching ratio, er x BR(pp — > W' — » tb), employing the invariant mass distri- 
bution of the Wbb final state. The Wbb mass distribution for the 57 candidate event 
sample is shown in Fig. 1501 and is compared with the expected mass distribution for 
a W' boson with mw' =500 GeV and the sum of the background processes. To esti- 
mate the size of a potential signal contribution, an unbinned maximum likelihood fit 
is performed to the observed mass distribution, allowing for both a signal and back- 
ground contribution for different values of mw ranging from 225 to 600 GeV. Based 
on the fit results for the fraction of events arising from W' production, the analysis 
sets 95 % C.L. upper limits using a Bayesian technique on the relative contribution 
of a W' boson for each value of mw', & x BR(W' — > tb)/a x BR(W' — > tb)sM, 
where the denominator is the expected a x B for a given W' boson mass assuming 
SM couplings. Figure E01 shows the upper limits on the W' boson production cross 
section as a function of the mw 1 for the two different assumptions on m VR . This 
analysis excludes a W' boson at 95 % C.L. with masses 225 < raw 1 <536 GeV for 
mw ^> m VR and 225 < mw> <566 GeV assuming raw < m un- 

8.3.3. Run II search for W' -> tb 

In Run II, D0 has conducted a comprehensive sea rch f or the W' boson in the top 
quark decay channel by analyzing 230 pb" 1 of datsP^l The search considers three 
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Fig. 50. CDF's Run I serach for W boson in top quark decayJSEl Left: The Wbb mass spectrum 
of the candidate events after constraining the lepton-neutrino invariant mass to the W boson mass. 
The distribution expected from the production of a W boson with a mass of 500 GeV is illustrated 
by the dashed curve. The distribution expected from the background processes is shown by the 
solid curve. Right: The upper limits on the W' boson production cross section as a function of the 
W' boson mass. Limits are shown for the case m w i S> M VR (solid) and m w i < M VR (dashed). 
The intercepts at [a.BR(W -* tb)/a.BR(W -> t5) SM ]=l correspond to the 95% C.L. limits on 
the W' boson mass with SM strength couplings. 



models of W' boson production. In each case, the CKM mixing matrix elements for 
the W boson is set to the SM values. The first model (W' L ) assumes the coupling 
of the W boson to SM fermions to be identical to that of the SM W boson. Under 
these assumptions, there is interference between the SM s-channel single top quark 
process and the W boson production process, although this interference term is 
small for large W boson masses. In the second and third model (W^), the W 
boson has only right-handed interactions. In the second model, the W' R boson is 
allowed to decay both to leptons and quarks, whereas in the third model it is only 
allowed to decay to quarks. The branching fraction for the decay W — > tb is about 
3/12 when the decays to quarks or l epton s are both included and is about 3/9 when 
the decay to leptons is not allowee r^l . COMPHEP 4.4.3 matrix element event 
generatoJ^ is used for the modeling of the W boson production process. Fig. 
[ST1 compares the invariant mass distribution for the W models with left-handed 
coupling (including interference) and right-handed coupling (no interference) with 
the SM s-channel single top quark produ ction , for a W boson mass of 600 GeV. 

Table [Ml shows the NLO cross sections^SI for a W boson x branching fraction 
to tb for three different W boson models. The systematic uncertainty on the cross 
section varies between about 12% at a mass of 600 GeV and 18% at a mass of 800 
GeV. The theoretical W' boson production cross se ction is more than 15 pb for 
masses between 200 and 400 GeV for all three model d '. The upper limits on the 
single top quark production cr oss se ction in the s-chann el we re 6.4 pb by D0 based 
on the analysis of same dataselPSl and 13.6 pb by which do not depend 
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Mftb) [GeV] 

Fig. 51. The invariant mass o f the top-bottom quark system at the parton level for different 
models of W' boson productiorJ223 Shown are the SM s-channel distribution, the W' L — > tb 
boson distribution, including the interference with the SM contribution, and the W R — > tb boson 
contribution, for a W' boson mass of 600 GeV. 

much on whether the W' boson coupling is left-handed or right-handed. Thus, the 
W' boson production decaying to tb is excluded in this mass region and the analysis 
explores the region of even higher masses for W' searches. 



Table 14. NLO production cross section for a W' boson X branching fraction to tb for 
three different W boson modelJEE] The production cross section for W' L interactions also 
include the SM s-channel contribution as well as the interference term between the two. 





Cross section X BR(W' - 


-* tb [pb] 




W mass [GeV] 


SM+W' L 


W' R (^lor q) 


W' R (-> q) only 


600 


2.17 


2.10 


2.79 


650 


1.43 


1.25 


1.65 


700 


1.03 


0.74 


0.97 


750 


0.76 


0.44 


0.57 


800 


0.65 


0.26 


0.34 



This search adopts a strategy similar to that used in CDF's Run I analysis, i.e. 
looks for events that are consistent with W' — > tb production and the W boson from 
the top quark decaying leptonically (W — > ez/, /iv; including W — > rv with r — > ev). 
It utilizes the same dataset, basic event selection, and background modeling as the 
D0 single top quark search described in Section I5TT1 and Ref. 11151 It selects signal- 
like events and separates the data into independent analysis sets based on final-state 
lepton flavor (electron or muon) with two or three jets and 6-tag multiplicity, single 
tagged or double tagged. The independent datasets are later combined in the final 
statistical analysis. 
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The large mass of the W boson is expected to set it apart from all background 
processes, hence a complete kinematic reconstruction of the invariant mass of the W' 
boson is performed by adding the four- vectors of all reconstructed final state objects: 
the jets, the lepton, and the neutrino. Figure 15121 shows the distribution of effective 
mass of £+ $t + jets when the Iv invariant mass is constrained to the W boson 
mass. The data, sum of background, and expected W boson contributions are shown 
for different couplings and masses. The observed event yield is consistent with the 
background predictions within uncertainties. The dominant sources of systematic 
uncertainty on the signal and background acceptances arises from the modeling of 
b hadrons in the simulation, the jet energy scale, object identification and trigger 
efficiencies, and modeling of jet fragmentation. 
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Fig. 52. D0 's Run II search for W' boson in top quark decayJ^QZ] The distribution of W boson 
invariant mass for several values of W boson mass and background processes for (a) left-handed 
W boson couplings, and (b) right-handed couplings when only the decay to quarks is allowed. 
Electron, muon, single-tagged and double-tagged events are combined. 

In the absence of a significant excess over the background predictions, a binned 
likelihood analysis is performed on the observed invariant mass distribution to ob- 
tain upper limits on the a x BR(pp — » W' — > tb) for discrete W mass points in 
each model. Figure [53] shows the a x BR(pp — > W' — > tb) limits together with the 
NLO cross sections and the expected limits, along with their uncertainties. At the 
95% C.L., the shaded areas above the solid lines are excluded in this analysis. The 
intersection of the solid line with the lower edge of the uncertainty band on the pre- 
dicted cross section defines the 95% C.L. lower mass limit for each model. Together 
with the limit from the SM s-channel single top quark search, the presence of a W 
boson with SM-like left-handed coupling is excluded if it has a mass between 200 
and 610 GeV. In addition, the presence of a W boson with right-handed couplings 
that is allowed to decay to leptons and quarks (only quarks) is excluded if it has 
mass between 200 and 630 GeV (670 GeV). This is the first direct search limit for 
the W boson which takes into account interference with the SM properly. It is also 
the most stringent limit on the presence of a W' boson in top quark decays. 
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Fig. 53. D0 's Run II search for W' boson in top quark decayJ223 Cross section times branching 
fraction limits at the 95% C.L. versus the mass of the W boson with (a) left-handed couplings 
and (b) right-handed couplings. Also shown are the NLO cross sections and the expected limits. 
The shaded regions above the circles are excluded by the measurement. 



More recently, CDF has also produced a preliminary result! ' on the search 
for resonant W boson which decays via W — » tb using a sample of approximately 
1 fb _1 . As in the previous analyses, the search looks for unexpected structure in 
the spectrum of the invariant mass distribution of the reconstructed W boson and 
two leading jets (Mwjj)- Expected contributions from SM processes are derived 
from selections and background studies of the single top quark analysis. Resonant 
tb production is modeled as simple W with SM-like couplings to fermions. The 
analysis finds no evidence for resonant W production and excludes a W boson for 
mw 1 <760 GeV in the case of mw' ^> m VR and raw' <790 GeV in the case of 
raw 1 < ra VR at 95% C.L.. 



9. Prospects for Future Measurements 
9.1. Tevatron synopsis 

The Tevatron program in Run II is currently undergoing rapid progress and mea- 
surements of the top quark are greatly improving. We have seen a flowering of 
techniques for the isolation and study of the top quark. Some approaches pursued 
in Run I, such as the topological selection of £+ jets events or the secondary vertex 
tagging have been adopted by both experiments as useful strategies. New selec- 
tion mechanisms such as £+ track in the dilepton case, and 6-jet tagging tagging 
based on neural network discriminant in general have produced competitive results 
and significantly expand the fraction of top quark events that are analyzable. At- 
tempts to maximize sensitivity for various measurements have generally meant that 
selections permit significant background yield. This necessitated improvements to 
background modeling that included increased ability to determine important prop- 
erties of background in data. An example involves the study of the flavor content 
of the b— tagged backgrounds. Usage of Monte Carlos has been improved by efforts 
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to more correctly describe jet production. 

Substantial increases in integrated luminosity have so far brought the measure- 
ment of the tt cross section from 30% to < 15% uncertainties. These measurements 
are in general agreement with the theoretically predicted value. Many different chan- 
nel definitions have been tried within the dilcpton, £+ jets and all-jets categories. So 
far, no serious discrepancies have been observed in the relative rates of these chan- 
nels. Extensive kinematic tests of the various channels also reveal no clear difference 
from standard model expectations: production and decay seem well-modeled. The 
Run II measurements will surpass the 10% level of precision. They are currently 
systematically dominated and of comparable precision as the theoretical uncertain- 
ties. This will hopefully motivate further theoretical work. It will also provide an 
interesting test when it comes to comparing LHC and Tevatron production. 

The major physical significance of m f has motivated both Tevatron experiments 
to devote considerable effort to its measurement. Already, Run II has surpassed 
initial expectations. Improved techniques have been a significant part of this effort. 
Both experiments have attempted a variety of fitting techniques to exploit the data 
to the fullest extent. A critical contribution in this regard has been the matrix 
element approach, which is now being used on all channels. Soon measurements in 
all channels will be dominated by systematic uncertainties. In Run I, it was clear 
from the outset that the calibration of jet energies would be a major concern. This 
still remains the case, and the measurement of top quark properties, particularly 
the mass, has driven major progress in this area. One of the results of this effort 
is that we now see the use of the top quark sample as a means to anchor this 
calibration via the use of jets from W — > jj decays in situ. This is a critical 
transition that will have significant impact on future Tevatron and LHC physics. 
Currently, the cumulative uncertainty in nit is 1%, with each experiment achieving 
1.5% in their best measurements. This result has been achieved with only a fraction 
of the expected total Run II luminosity. As energy scale and other uncertainties 
are further reduced, the measurement will likely become even more precise. Such a 
measurement will remain competitive with the LHC for some time to come. 

An exciting aspect of the field is the rising experimental capability to test elec- 
troweak physics through the properties of this particle. With the advent of a signal 
for single top quark production, a new direct probe of electroweak interactions is 
available. The measurement of \Vtt\ is an important element of this effort. Run II 
should yield a 10% measurement of this quantity. Strong constraints have been laid 
on potential V + A couplings at the Wtb vertex. The top quark charge is already 
determined to be like its c and u brethren. 

Direct searches of new physics in the production and decay of the top quark have 
yielded significant constraints. Techniques to look for charged Higgs Bosons in top 
quark decay have been refined as the calculational picture has improved. Limits are 
occupying an increasing fraction of the experimentally available parameter space. 
With the increased statistics, better limits of it resonance production are being 
set, as well as on the presence of potential extra W bosons in single top quark 
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production. 

9.2. Looking to the Future 

Many of the basic techniques described in this review can be used at the LHC, 
and there will be many new challenges. Instantaneous luminosities of 10 x relative 
to the Tevatron, and tt cross sections 100 x larger will provide a thousand- fold in- 
crease in statistical samples. This will allow, for instance, very large event samples 
with both b— jets identified. Early studies indicate that 1% precision in nit mea- 
surement is accessible (e.g. Ref. 209). To achieve or surpass this level will require 
large, well-constrained samples. The top quark pair production and decay will be 
measurable across a wide array of channels. The measurements that will improve 
most dramatically are the other top quark properties, such as Vtt or W boson he- 
licity. The searches for new physics will also gain considerable sensitivity from the 
ability to make precise analysis of the rates and kinematics of final state particles 
in top quark events. 




m, [GeV] m H [GeV] 



Fig. 54. The expected region of Higgs mass in the mw vs. mt plane (left). The combined 
Tevatron measurement of mt (and my^) are shown. Also shown is the fit to all electroweak data 
(right), in which the value of mt provides substantial weight EL25J 

Finally, the study of the top quark is inextricably linked with the question of elec- 
troweak symmetry breaking and the Higgs mechanism. By including the measured 
value of m t in the global electroweak fit one obtains the most stringent con straint 
on the mass of the Higgs boson. The preferred value for mn is 76^4 

GeV^, which 

is not favored by the ran > 114.7 GeV limit from direct searches at LEP-rP223. 
This result is shown in Fig. EU The 95% confidence level upper limit on the Higgs 
boson mass set by this analysis is m# < 144 GeV. It is interesting to consider the 
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measurements in relation to potential new physics. For example , Figure 1551 indicates 
the regions in the my/ vs. mt plane favored by MSSMs 1 1 [ and by the standard 
model. If current central values of mw and m t remain the same as more precise 
measurements are made, a signature of new physics may arise. In the standard 
model, improving the precision with which the top quark mass is measured will be 
a crucial ingredient in further tightening the constraints on the Higgs boson mass. 

The most favored region provides for a complex set of Higgs boson production 
and decay processes. In this sense, the study of the Higgs boson has much in com- 
mon with what confronted experimentalists for the top quark discovery. Final state 
channels including most of the objects identifiable at a hadron collider detector: 
e, 7, fj,,y,jets,T, and b— jets. This implies that the full range of experimental capa- 
bilities at the experiments will need to be mastered to maximize sensitivity. Multi- 
ple channels will be used to gain significance before any one channel sees a pristine 
signal. Decay chains include gauge bosons, and there are concurrent W/Z + X 
backgrounds that go with this. Some of the modeling questions and background 
estimation approaches that have been important in study of the top quark will also 
be important for the Higgs boson search. 
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Fig. 55. Regions in mw vs - m t plane favored by MSSMs and by the standard model 
Anticipated sensitivities from future measurements are also shown. 



One thing is clear, the race is on in earnest for the Higgs boson. Will the Teva- 
tron have a shot at a signal? Will the LHC effort turn on rapidly and overwhelm 
the search quickly? Is there a Higgs mechanism at all, or are there new strong dy- 
namics in its place which elevate the top quark to a more fundamental status. The 
study of the top quark lies at an important cross-roads in the pursuit of answers to 
these questions. We look forward with fascination to the chance to work with our 
colleagues in the coming years' study. 
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